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

ATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAT  >  NC_000913/1231808‑1232103
                                                                                                                                                     |                                                                                                                                                   
aTACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCa                                                                                                                                                    >  2:253658/1‑150 (MQ=255)
          ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATcc                                                                                                                                          >  1:1067459/1‑150 (MQ=255)
          ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATcc                                                                                                                                          <  2:32865/150‑1 (MQ=255)
             gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTc                                                                                                                                       <  2:1075994/150‑1 (MQ=255)
                tATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGAt                                                                                                                                    >  1:862321/1‑150 (MQ=255)
                  tCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTTAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGTGATGc                                                                                                                                  <  2:505382/150‑1 (MQ=255)
                   ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCg                                                                                                                                 >  2:926703/1‑150 (MQ=255)
                   ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCg                                                                                                                                 >  2:1416942/1‑150 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGTGATGCgtt                                                                                                                               <  2:134552/150‑3 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGc                                                                                                                               >  2:1273307/1‑150 (MQ=255)
                     gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGAGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGc                                                                                                                               >  2:847511/1‑150 (MQ=255)
                      tacaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCa                                                                                                                              <  2:1775636/147‑1 (MQ=255)
                      tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGTGATGCgttc                                                                                                                              <  2:6397/150‑4 (MQ=255)
                        accATTTTAATGTCGTGCTCGACAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAAGTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGAATATCCGTCGATGCGGCAgg                                                                                                                            <  1:1596615/148‑1 (MQ=255)
                        aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAgg                                                                                                                            >  2:769767/1‑150 (MQ=255)
                         aCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGc                                                                                                                           <  2:474621/150‑1 (MQ=255)
                          cATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGcc                                                                                                                          >  1:371631/1‑150 (MQ=255)
                                aaTGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGt                                                                                                                    <  2:1158341/150‑1 (MQ=255)
                                  tGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAg                                                                                                                  <  2:1514311/150‑1 (MQ=255)
                                     cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                               <  1:563441/150‑1 (MQ=255)
                                     cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                               <  2:222437/150‑1 (MQ=255)
                                         cTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAAGTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                           >  1:1732274/1‑150 (MQ=255)
                                          tCGAAAGAACGGTGCGTGAACTGCGCGGCGACCCCTGTTTGCACCTGGAAGAGTTTGCCCCGACGAAGCCGGAAATTATCTGTTCCCGCGCGTTTGGTGAACGCATCACGGATTATCCGTCGATGAGGCAGGCCATGTGGAGTGACGCTg                                                                                                          >  2:85126/1‑150 (MQ=255)
                                          tCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTg                                                                                                          >  2:1361493/1‑150 (MQ=255)
                                           cGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                         >  2:1797784/1‑150 (MQ=255)
                                            gAAAGAACGGTGCGTGAACTGCGCGGCGACCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGATGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTACTTACGCTccc                                                                                                        >  1:1615201/1‑150 (MQ=255)
                                            gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGcc                                                                                                        >  1:615490/1‑150 (MQ=255)
                                              aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCg                                                                                                      >  2:330840/1‑150 (MQ=255)
                                              aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCg                                                                                                      <  2:371631/150‑1 (MQ=255)
                                                  aCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcg                                                                                                  <  1:1462737/150‑1 (MQ=255)
                                                   cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                 >  1:1298328/1‑150 (MQ=255)
                                                   cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                 <  1:727599/150‑1 (MQ=255)
                                                    ggTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcggc                                                                                                >  2:150271/1‑150 (MQ=255)
                                                     gTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGgcggcg                                                                                               <  1:519066/150‑1 (MQ=255)
                                                      tGCGGGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGcggcgg                                                                                              <  2:1805346/150‑1 (MQ=255)
                                                       gCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCGCCGACGAAGCAGGAAATTATCTGTTCCCGCACGTTTGGTGAACGCATCACGGATTAGCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCGGGGCGGCGGa                                                                                             >  1:319001/1‑150 (MQ=255)
                                                        cgggAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaa                                                                                            <  1:900345/147‑1 (MQ=255)
                                                         gTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGGTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaa                                                                                           >  1:1203695/1‑150 (MQ=255)
                                                          tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaa                                                                                          <  2:966215/150‑1 (MQ=255)
                                                           gAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaaa                                                                                         >  2:712372/1‑150 (MQ=255)
                                                            aaCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAc                                                                                        <  2:597081/150‑1 (MQ=255)
                                                              cTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACtt                                                                                      >  2:557687/1‑150 (MQ=255)
                                                               tGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTc                                                                                     <  1:462322/150‑1 (MQ=255)
                                                                gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGATCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCg                                                                                    <  1:511766/150‑1 (MQ=255)
                                                                 cgcgGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGc                                                                                   >  1:1435357/1‑150 (MQ=255)
                                                                        aaCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGc                                                                            <  2:1564567/150‑1 (MQ=255)
                                                                         aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                           <  1:1014405/150‑1 (MQ=255)
                                                                         aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGACCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                           >  2:17509/1‑150 (MQ=255)
                                                                           ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATc                                                                         >  2:115425/1‑150 (MQ=255)
                                                                            cgcggcgCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCa                                                                        <  2:387431/144‑1 (MQ=255)
                                                                                 tGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGTTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                   >  1:464218/1‑150 (MQ=255)
                                                                                   cAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                 <  1:1796709/150‑1 (MQ=255)
                                                                                     aCTGGTATAGTTTGCACCGACGAAGCAGGAATTTATCTGTTCCCGCTCGGTTTGTGAACGCATCACTGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGATGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCg                                                               <  1:309196/150‑1 (MQ=255)
                                                                                         gtagagTTTGCACCGACGAAGCAGGAAATTATGTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTTCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                           <  1:336781/148‑1 (MQ=255)
                                                                                         gatgagTTTGCAGCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCTCGGATTATCCGGCGATGCGGCAGGCCATTTGTAGTTACGGTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCTTCAATATTGTCGGttt                                                           <  2:1484584/147‑1 (MQ=255)
                                                                                           agagTTTGCCCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                         >  2:170262/1‑150 (MQ=255)
                                                                                               tttGCACCGACGAAGCAGGAATTTATGTGTTCCCGCTCGTTTGGTGAACGCATATCGGTTTATCCGTCGATGCGGCAGGCGATTTGTAGTTACGCTGCCAGGGCGGCGGAAAATCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTcc                                                     <  1:916881/150‑1 (MQ=255)
                                                                                                   cACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGt                                                 >  2:646866/1‑150 (MQ=255)
                                                                                                    aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGtt                                                <  1:1559330/150‑1 (MQ=255)
                                                                                                    aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGtt                                                >  2:457948/1‑150 (MQ=255)
                                                                                                        acgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatt                                            <  2:1824045/150‑1 (MQ=255)
                                                                                                         cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatta                                           >  2:623070/1‑150 (MQ=255)
                                                                                                         cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGAGGAAAAACTTCGCAGAGAGCATCAATATTGTCGGTTTATCTCCACGTttatta                                           >  1:1454480/1‑150 (MQ=255)
                                                                                                         cgaAGCAGGAAAGTATCTGATCCCGCTCGTTTGGTTAACGCATCACTGATTATCCGTCGATGCGTCAGGCCTTTTGTAGTTACGATTACCGGTCTGCGGAAAAACTTCGCAGCGAGCAGCAATATTGTCGGATTATCTCCACGTttatta                                           <  1:116399/150‑1 (MQ=255)
                                                                                                             gCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGAc                                       >  2:157814/1‑150 (MQ=255)
                                                                                                             gCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGAc                                       <  2:862321/150‑1 (MQ=255)
                                                                                                              cattaaaTTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCATGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACg                                      <  2:1391509/146‑1 (MQ=255)
                                                                                                              cAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACg                                      >  1:444989/1‑150 (MQ=255)
                                                                                                               aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                     >  1:1756712/1‑150 (MQ=255)
                                                                                                               aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                     >  1:401906/1‑150 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  1:1179830/150‑1 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  1:253658/150‑1 (MQ=255)
                                                                                                                ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                    <  2:1298328/150‑1 (MQ=255)
                                                                                                                 gAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                   >  1:982237/1‑150 (MQ=255)
                                                                                                                 gAAAGTAGCTGTTCCCGCTCGTTTGGTGAACGCATCACGTATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                   <  1:299475/150‑1 (MQ=255)
                                                                                                                     ttATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCCTTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGCGCATCTATATTGTCGGTTTAGCTCCACGTTTATTAAGACGTCACCAt                               >  1:98851/1‑150 (MQ=255)
                                                                                                                      tATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAtt                              >  1:1153490/1‑150 (MQ=255)
                                                                                                                       aTCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAttt                             >  2:348139/1‑150 (MQ=255)
                                                                                                                         cTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgc                           >  1:1089175/1‑150 (MQ=255)
                                                                                                                          tGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgcg                          >  1:669557/1‑150 (MQ=255)
                                                                                                                           gTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgcgc                         <  1:662950/150‑1 (MQ=255)
                                                                                                                             tCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTc                       >  2:642239/1‑150 (MQ=255)
                                                                                                                               ccGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACTCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACAATTTGCGCTCaa                     >  2:1597885/1‑150 (MQ=255)
                                                                                                                                cGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                    >  2:966389/1‑150 (MQ=255)
                                                                                                                                cGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                    >  2:712395/1‑150 (MQ=255)
                                                                                                                                     gTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAcc               <  1:197207/150‑1 (MQ=255)
                                                                                                                                         ggTGAACGCATCACGGATTATCCGTCGATGCGGCATGCCATTTGTAGTTTCGCTGCCCGTGCGGCTGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTAACCATTTGCGCTCAATGAACCttat           <  1:847511/150‑1 (MQ=255)
                                                                                                                                          gTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          >  1:687996/1‑150 (MQ=255)
                                                                                                                                          gTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          <  1:926703/150‑1 (MQ=255)
                                                                                                                                          cctgccGCATCGACGGATAATCCGTGATGCGGCAGGCCATTTGTTGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatt          >  2:782769/6‑150 (MQ=255)
                                                                                                                                           tGAACGCATCACGGATTATCCGTCGATGCGGCAGACCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatta         <  1:17509/150‑1 (MQ=255)
                                                                                                                                             gccGCATCGACGGATAATCCGTGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACg       >  2:1643550/3‑150 (MQ=255)
                                                                                                                                               cGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc     <  2:615490/150‑1 (MQ=255)
                                                                                                                                               cGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc     <  1:933157/150‑1 (MQ=255)
                                                                                                                                                gCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAATTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATCAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCa    <  2:336160/150‑1 (MQ=255)
                                                                                                                                                  aTCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt  >  1:368761/1‑150 (MQ=255)
                                                                                                                                                  aTCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGAGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt  <  1:115425/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAT  >  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: