Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,412,516 C→A L44I (CTA→ATA)  yjfM → DUF1190 domain‑containing protein YjfM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,412,5160CA99.4% 588.2 / ‑4.1 169L44I (CTA→ATA) yjfMDUF1190 domain‑containing protein YjfM
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (80/88);  total (81/88)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 171 total aligned reads displayed.

GACGATCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGCTAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTGCGCGC  >  NC_000913/4412380‑4412651
                                                                                                                                        |                                                                                                                                       
     tcctGGCCAGTAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGATCGATTAGTCGCATTGTTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  2:2304269/4‑140 (MQ=255)
     tCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  2:329684/1‑140 (MQ=255)
     tCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  2:320097/1‑140 (MQ=255)
     tCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  2:569549/1‑140 (MQ=255)
     tCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  2:2284486/1‑140 (MQ=255)
     tCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGc                                                                                                                                 >  1:1508626/1‑140 (MQ=255)
        tGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGt                                                                                                                              >  1:1089921/1‑140 (MQ=255)
         ggCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTg                                                                                                                             <  1:880045/140‑1 (MQ=255)
               aaaaCGCAAATCAAGAAATAACAGTAAAATTGGTCACGGTGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCtt                                                                                                                       >  2:824103/1‑140 (MQ=255)
               aaaaCGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGc                                                                                                                         >  2:2179893/1‑138 (MQ=255)
               aaaaCGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCtt                                                                                                                       >  2:651219/1‑140 (MQ=255)
               aaaaCGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCtt                                                                                                                       >  1:946791/1‑140 (MQ=255)
               aaaaCGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCtt                                                                                                                       >  1:1303425/1‑140 (MQ=255)
                    gCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTc                                                                                                                  <  2:187747/140‑1 (MQ=255)
                    gCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTc                                                                                                                  <  1:1778179/140‑1 (MQ=255)
                    gCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTc                                                                                                                  <  1:208056/140‑1 (MQ=255)
                     cAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCg                                                                                                                 <  1:1950342/140‑1 (MQ=255)
                              aaaTAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAg                                                                                                        <  2:1575392/140‑1 (MQ=255)
                                  aaCAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTg                                                                                                    >  1:615532/1‑140 (MQ=255)
                                     aGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAg                                                                                                 <  2:1508626/140‑1 (MQ=255)
                                       tAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCg                                                                                               >  1:1262281/1‑140 (MQ=255)
                                       tAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCg                                                                                               >  1:1021502/1‑140 (MQ=255)
                                           aTTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGAtag                                                                                           >  1:2244719/1‑140 (MQ=255)
                                            ttGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGAtagt                                                                                          >  2:2016905/1‑140 (MQ=255)
                                               gTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGc                                                                                       >  1:1775343/1‑140 (MQ=255)
                                               gTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGc                                                                                       >  2:726926/1‑140 (MQ=255)
                                               gTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGc                                                                                       >  2:886483/1‑140 (MQ=255)
                                                tCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCg                                                                                      <  1:601646/140‑1 (MQ=255)
                                                tCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCg                                                                                      <  2:1368422/140‑1 (MQ=255)
                                                tCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCg                                                                                      <  1:1140523/140‑1 (MQ=255)
                                                 cACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGa                                                                                     >  1:845884/1‑140 (MQ=255)
                                                 cACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGa                                                                                     >  2:972621/1‑140 (MQ=255)
                                                  aCGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGAt                                                                                    <  1:1144655/140‑1 (MQ=255)
                                                   cGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTTCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATg                                                                                   >  2:511936/1‑140 (MQ=255)
                                                          aTTAGTCGCATTGGTAGACCGAATAATCCTTTTTAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATaa                                                                            >  1:69982/1‑140 (MQ=255)
                                                                cGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAAcgacgg                                                                      >  1:2335787/1‑140 (MQ=255)
                                                                   aTTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAAcgacggcga                                                                   <  1:2284486/140‑1 (MQ=255)
                                                                     tGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAAcgacggcgacg                                                                 <  2:575994/140‑1 (MQ=255)
                                                                            ccGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGtt                                                          >  2:1058285/1‑140 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAacgac                                                                        <  2:9126/123‑1 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTa                                                       >  1:1247525/1‑140 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTa                                                       >  1:276970/1‑140 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTa                                                       >  2:1274380/1‑140 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTa                                                       >  2:1718961/1‑140 (MQ=255)
                                                                               aataatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGATGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTa                                                       >  2:2354635/1‑140 (MQ=255)
                                                                                 taatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTAcg                                                     <  2:2130399/140‑1 (MQ=255)
                                                                                 taatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTAcg                                                     <  1:1477324/140‑1 (MQ=255)
                                                                                  aatCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTAcgc                                                    <  2:36527/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  2:2214346/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:1924833/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  2:276970/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:2324130/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  2:931020/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:743803/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:726926/140‑1 (MQ=255)
                                                                                    tCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  2:615532/140‑1 (MQ=255)
                                                                                      cTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACa                                                <  2:1946652/140‑1 (MQ=255)
                                                                                      cTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACa                                                >  1:2304097/1‑140 (MQ=255)
                                                                                      cTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACa                                                <  1:511936/140‑1 (MQ=255)
                                                                                        ttttAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGt                                              <  2:69982/140‑1 (MQ=255)
                                                                                         ttGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             <  2:845884/140‑1 (MQ=255)
                                                                                         ttGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             <  1:219515/140‑1 (MQ=255)
                                                                                         ttGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             >  2:2337533/1‑140 (MQ=255)
                                                                                         ttGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             <  1:893082/140‑1 (MQ=255)
                                                                                         ttGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             <  1:842116/140‑1 (MQ=255)
                                                                                           gAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCa                                           <  1:989243/140‑1 (MQ=255)
                                                                                           gAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCa                                           <  2:1775343/140‑1 (MQ=255)
                                                                                              ccGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  2:1085012/140‑1 (MQ=255)
                                                                                              ccGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  1:2337533/140‑1 (MQ=255)
                                                                                              ccGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  1:635601/140‑1 (MQ=255)
                                                                                              ccGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  2:898548/140‑1 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGGACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTg                                     >  2:1508714/1‑140 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGCTAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTg                                     >  1:256297/1‑140 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  >  2:2224776/1‑127 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  >  2:2255855/1‑127 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:2224776/127‑1 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGa                                                  <  1:2255855/127‑1 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTg                                             <  1:916371/132‑1 (MQ=255)
                                                                                                 tGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  2:425193/137‑1 (MQ=255)
                                                                                                  gCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGt                                    >  1:1652966/1‑140 (MQ=255)
                                                                                                       aaTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtga                               <  2:422291/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:2034772/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:2217785/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:1058285/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:320097/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:330993/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  2:1557381/140‑1 (MQ=255)
                                                                                                        aTCGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              >  2:209815/1‑140 (MQ=255)
                                                                                                        aTCGCGACGCACAGAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTCAACGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACAGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              <  1:2380349/140‑1 (MQ=255)
                                                                                                          cGCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgatga                            <  2:2374114/140‑1 (MQ=255)
                                                                                                           gCTACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATtgat                              >  1:393806/1‑137 (MQ=255)
                                                                                                             tACGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACgg                         >  2:853270/1‑140 (MQ=255)
                                                                                                              aCGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGa                                        <  2:125905/124‑1 (MQ=255)
                                                                                                               cGCACAAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGTCGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGta                       <  2:1127640/140‑1 (MQ=255)
                                                                                                                   caAAAATACTTAACGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGtaataa                   >  1:1391005/1‑140 (MQ=255)
                                                                                                                                cGTTGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTgc      >  1:242295/1‑140 (MQ=255)
                                                                                                                                   tGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTgcgcg   >  1:911238/1‑140 (MQ=255)
                                                                                                                                   tGGCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTgcgcg   <  2:1303425/140‑1 (MQ=255)
                                                                                                                                    ggCGATAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTgcgcgc  <  1:1562101/140‑1 (MQ=255)
                                                                                                                                    ggCGATAATGGGCGGTGCCGCTTTCTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTgcgcgc  <  2:71509/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GACGATCATGGCCAGAAAACGCAAATCAAGAAATAACAGTAAAATTGGTCACGGAGCGATTAGTCGCATTGGTAGACCGAATAATCCTTTTGAACCGTGCCGCAATCGCTACGCACAAAAATACTTAACGTTGGCGCTAATGGGCGGTGCCGCTTTTTTCGTATTGAAAGGTTGTAGCGATAGTAGCGATGTCGATAACGACGGCGACGGAACGTTTTACGCGACAGTGCAGGATTGTATTGATGACGGTAATAATGCCGACATTTGCGCGC  >  NC_000913/4412380‑4412651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: