Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,956:1 +G 5.2% coding (771/1269 nt) umuC → translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,231,9561.G 5.2% 396.8 / 17.3 96coding (771/1269 nt)umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base . (46/45);  new base G (2/3);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

ATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAT  >  NC_000913/1231808‑1232103
                                                                                                                                                     |                                                                                                                                                   
aTACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑a                                                                                                                                                    >  2:253658/1‑150 (MQ=255)
          ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATcc                                                                                                                                          >  1:1067459/1‑150 (MQ=255)
          ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATcc                                                                                                                                          <  2:32865/150‑1 (MQ=255)
             gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTc                                                                                                                                       <  2:1075994/150‑1 (MQ=255)
                tATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                    >  1:862321/1‑150 (MQ=255)
                  tCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTTAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGc                                                                                                                                  <  2:505382/150‑1 (MQ=255)
                   ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCg                                                                                                                                 >  2:926703/1‑150 (MQ=255)
                   ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCg                                                                                                                                 >  2:1416942/1‑150 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgtt                                                                                                                               <  2:134552/150‑3 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGc                                                                                                                               >  2:1273307/1‑150 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGAGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGc                                                                                                                               >  2:847511/1‑150 (MQ=255)
                      tacaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCa                                                                                                                              <  2:1775636/147‑1 (MQ=255)
                      tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgttc                                                                                                                              <  2:6397/150‑4 (MQ=255)
                        accATTTTAATGTCGTGCTCGACAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAAGTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGAATATCCGTCGATGCGGCAgg                                                                                                                            <  1:1596615/148‑1 (MQ=255)
                        aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAgg                                                                                                                            >  2:769767/1‑150 (MQ=255)
                         aCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGc                                                                                                                           <  2:474621/150‑1 (MQ=255)
                          cATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGcc                                                                                                                          >  1:371631/1‑150 (MQ=255)
                                aaTGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGt                                                                                                                    <  2:1158341/150‑1 (MQ=255)
                                  tGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAg                                                                                                                  <  2:1514311/150‑1 (MQ=255)
                                     cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                               <  1:563441/150‑1 (MQ=255)
                                     cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                               <  2:222437/150‑1 (MQ=255)
                                         cTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAAGTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                           >  1:1732274/1‑150 (MQ=255)
                                          tCGAAAGAACGGTGCGTGAACTGCGCGGCGACCCCTGTTTGCACCTGGAAGAGTTTGCCCCGACGAAGCCGGAAATTATCTGTTCCCGCGCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGAGGCAGGCCATGTGGAGTGACGCTg                                                                                                          >  2:85126/1‑150 (MQ=255)
                                          tCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTg                                                                                                          >  2:1361493/1‑150 (MQ=255)
                                           cGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                         >  2:1797784/1‑150 (MQ=255)
                                            gAAAGAACGGTGCGTGAACTGCGCGGCGACCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGATGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTACTTACGCTccc                                                                                                        >  1:1615201/1‑150 (MQ=255)
                                            gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGcc                                                                                                        >  1:615490/1‑150 (MQ=255)
                                              aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCg                                                                                                      >  2:330840/1‑150 (MQ=255)
                                              aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCg                                                                                                      <  2:371631/150‑1 (MQ=255)
                                                  aCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcg                                                                                                  <  1:1462737/150‑1 (MQ=255)
                                                   cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                 >  1:1298328/1‑150 (MQ=255)
                                                   cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                 <  1:727599/150‑1 (MQ=255)
                                                    ggTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcggc                                                                                                >  2:150271/1‑150 (MQ=255)
                                                     gTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGgcggcg                                                                                               <  1:519066/150‑1 (MQ=255)
                                                      tGCGGGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGcggcgg                                                                                              <  2:1805346/150‑1 (MQ=255)
                                                       gCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCGCCGACGAAGCAGGAAATTATCTGTTCCCGCACGTTTGGTGAACGCATC‑ACGGATTAGCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCGGGGCGGCGGa                                                                                             >  1:319001/1‑150 (MQ=255)
                                                        cgggAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaa                                                                                            <  1:900345/147‑1 (MQ=255)
                                                         gTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGGTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaa                                                                                           >  1:1203695/1‑150 (MQ=255)
                                                          tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaa                                                                                          <  2:966215/150‑1 (MQ=255)
                                                           gAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaaa                                                                                         >  2:712372/1‑150 (MQ=255)
                                                            aaCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAc                                                                                        <  2:597081/150‑1 (MQ=255)
                                                              cTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACtt                                                                                      >  2:557687/1‑150 (MQ=255)
                                                               tGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTc                                                                                     <  1:462322/150‑1 (MQ=255)
                                                                gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGATCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCg                                                                                    <  1:511766/150‑1 (MQ=255)
                                                                 cgcgGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGc                                                                                   >  1:1435357/1‑150 (MQ=255)
                                                                        aaCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGc                                                                            <  2:1564567/150‑1 (MQ=255)
                                                                         aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                           <  1:1014405/150‑1 (MQ=255)
                                                                         aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGACCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                           >  2:17509/1‑150 (MQ=255)
                                                                           ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATc                                                                         >  2:115425/1‑150 (MQ=255)
                                                                            cgcggcgCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCa                                                                        <  2:387431/144‑1 (MQ=255)
                                                                                 tGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGTTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                   >  1:464218/1‑150 (MQ=255)
                                                                                   cAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                 <  1:1796709/150‑1 (MQ=255)
                                                                                     aCTGGTATAGTTTGCACCGACGAAGCAGGAATTTATCTGTTCCCGCTCGGTTTGTGAACGCATC‑ACTGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGATGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCg                                                               <  1:309196/150‑1 (MQ=255)
                                                                                         gtagagTTTGCACCGACGAAGCAGGAAATTATGTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTTCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                           <  1:336781/148‑1 (MQ=255)
                                                                                         gatgagTTTGCAGCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑TCGGATTATCCGGCGATGCGGCAGGCCATTTGTAGTTACGGTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCTTCAATATTGTCGGttt                                                           <  2:1484584/147‑1 (MQ=255)
                                                                                           agagTTTGCCCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                         >  2:170262/1‑150 (MQ=255)
                                                                                               tttGCACCGACGAAGCAGGAATTTATGTGTTCCCGCTCGTTTGGTGAACGCATA‑TCGGTTTATCCGTCGATGCGGCAGGCGATTTGTAGTTACGCTGCCAGGGCGGCGGAAAATCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTcc                                                     <  1:916881/150‑1 (MQ=255)
                                                                                                   cACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGt                                                 >  2:646866/1‑150 (MQ=255)
                                                                                                    aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGtt                                                <  1:1559330/150‑1 (MQ=255)
                                                                                                    aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGtt                                                >  2:457948/1‑150 (MQ=255)
                                                                                                        acgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatt                                            <  2:1824045/150‑1 (MQ=255)
                                                                                                         cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatta                                           >  2:623070/1‑150 (MQ=255)
                                                                                                         cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGAGGAAAAACTTCGCAGAGAGCATCAATATTGTCGGTTTATCTCCACGTttatta                                           >  1:1454480/1‑150 (MQ=255)
                                                                                                         cgaAGCAGGAAAGTATCTGATCCCGCTCGTTTGGTTAACGCATC‑ACTGATTATCCGTCGATGCGTCAGGCCTTTTGTAGTTACGATTACCGGTCTGCGGAAAAACTTCGCAGCGAGCAGCAATATTGTCGGATTATCTCCACGTttatta                                           <  1:116399/150‑1 (MQ=255)
                                                                                                             gCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGAc                                       >  2:157814/1‑150 (MQ=255)
                                                                                                             gCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGAc                                       <  2:862321/150‑1 (MQ=255)
                                                                                                              cattaaaTTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCATGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACg                                      <  2:1391509/146‑1 (MQ=255)
                                                                                                              cAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACg                                      >  1:444989/1‑150 (MQ=255)
                                                                                                               aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                     >  1:1756712/1‑150 (MQ=255)
                                                                                                               aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                     >  1:401906/1‑150 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  1:1179830/150‑1 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  1:253658/150‑1 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  2:1298328/150‑1 (MQ=255)
                                                                                                                 gAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                   >  1:982237/1‑150 (MQ=255)
                                                                                                                 gAAAGTAGCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGTATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                   <  1:299475/150‑1 (MQ=255)
                                                                                                                     ttATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCCTTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGCGCATCTATATTGTCGGTTTAGCTCCACGTTTATTAAGACGTCACCAt                               >  1:98851/1‑150 (MQ=255)
                                                                                                                      tATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAtt                              >  1:1153490/1‑150 (MQ=255)
                                                                                                                       aTCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAttt                             >  2:348139/1‑150 (MQ=255)
                                                                                                                         cTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgc                           >  1:1089175/1‑150 (MQ=255)
                                                                                                                          tGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgcg                          >  1:669557/1‑150 (MQ=255)
                                                                                                                           gTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgcgc                         <  1:662950/150‑1 (MQ=255)
                                                                                                                             tCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTc                       >  2:642239/1‑150 (MQ=255)
                                                                                                                               ccGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACTCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACAATTTGCGCTCaa                     >  2:1597885/1‑150 (MQ=255)
                                                                                                                                cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                    >  2:966389/1‑150 (MQ=255)
                                                                                                                                cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                    >  2:712395/1‑150 (MQ=255)
                                                                                                                                     gTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAcc               <  1:197207/150‑1 (MQ=255)
                                                                                                                                         ggTGAACGCATC‑ACGGATTATCCGTCGATGCGGCATGCCATTTGTAGTTTCGCTGCCCGTGCGGCTGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTAACCATTTGCGCTCAATGAACCttat           <  1:847511/150‑1 (MQ=255)
                                                                                                                                          gTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          >  1:687996/1‑150 (MQ=255)
                                                                                                                                          gTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          <  1:926703/150‑1 (MQ=255)
                                                                                                                                          cctgccGCATCGACGGATAATCCGT‑GATGCGGCAGGCCATTTGTTGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          >  2:782769/6‑150 (MQ=255)
                                                                                                                                           tGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGACCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatta         <  1:17509/150‑1 (MQ=255)
                                                                                                                                             gccGCATCGACGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACg       >  2:1643550/3‑150 (MQ=255)
                                                                                                                                               cGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc     <  2:615490/150‑1 (MQ=255)
                                                                                                                                               cGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc     <  1:933157/150‑1 (MQ=255)
                                                                                                                                                gCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAATTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATCAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCa    <  2:336160/150‑1 (MQ=255)
                                                                                                                                                  aTC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt  >  1:368761/1‑150 (MQ=255)
                                                                                                                                                  aTC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGAGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt  <  1:115425/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAT  >  NC_000913/1231808‑1232103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑