Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,859,354 T→C 48.7% N255D (AAC→GAC)  yidK ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,859,3540TC48.7% 6.9 / inf 284N255D (AAC→GAC) yidKputative transporter
Reads supporting (aligned to +/- strand):  ref base T (75/70);  new base C (62/76);  total (138/146)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 290 total aligned reads displayed.

CCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGG  >  NC_000913/3859218‑3859493
                                                                                                                                        |                                                                                                                                           
   ggCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAg                                                                                                                                       >  1:908233/1‑140 (MQ=255)
   ggCAGCACCAGAACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCAGATTTGTACACCAGTAAAAGGTGTCCACCAg                                                                                                                                       >  1:1615195/1‑140 (MQ=255)
        cacaAGTACCAGCGGGTCGAGCATTTTCAGCACAGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACTTTGCCCTGATTTTTACACCAGTAATAGGTGTTCACCAGAATCa                                                                                                                                  <  2:1208355/140‑1 (MQ=255)
        caCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCa                                                                                                                                  <  2:173244/140‑1 (MQ=255)
        caCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCa                                                                                                                                  <  1:293600/140‑1 (MQ=255)
         aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCaa                                                                                                                                 >  2:849657/1‑140 (MQ=255)
          ccAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCaaa                                                                                                                                <  2:1394339/140‑1 (MQ=255)
          ccAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCaaa                                                                                                                                <  1:412568/140‑1 (MQ=255)
                   gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                       >  2:1443067/1‑140 (MQ=255)
                   gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                       >  2:1181400/1‑140 (MQ=255)
                   gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAAt                                                                                                                       >  2:456554/1‑140 (MQ=255)
                    cGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATg                                                                                                                      >  1:1357098/1‑140 (MQ=255)
                        tCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTACAAGGTGTCCACCAGAATCAAACCGGTAAATGCcgc                                                                                                                  >  1:528840/1‑140 (MQ=255)
                        tCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCcgc                                                                                                                  >  1:1634959/1‑140 (MQ=255)
                        tCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCcgc                                                                                                                  >  1:1172473/1‑140 (MQ=255)
                        tCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCcgc                                                                                                                  >  1:866269/1‑140 (MQ=255)
                        tCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCcgc                                                                                                                  >  2:1036780/1‑140 (MQ=255)
                           aGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGcc                                                                                                               >  1:1845608/1‑140 (MQ=255)
                           aGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGcc                                                                                                               >  1:757205/1‑140 (MQ=255)
                                  tCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGc                                                                                                        <  2:949998/140‑1 (MQ=255)
                                  tCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGc                                                                                                        <  1:577729/140‑1 (MQ=255)
                                   cAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCa                                                                                                       <  2:1493419/140‑1 (MQ=255)
                                      cACCGCCGTTAACAGCGCCTCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAgg                                                                                                    >  2:709607/1‑140 (MQ=255)
                                      cACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAAACGGCAAgg                                                                                                    >  2:1812171/1‑140 (MQ=255)
                                        ccgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGa                                                                                                  <  2:474346/140‑1 (MQ=255)
                                        ccgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGa                                                                                                  <  1:365975/140‑1 (MQ=255)
                                         cgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAt                                                                                                 >  1:1704918/1‑140 (MQ=255)
                                                     cgcCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCAc                                                                                     >  1:1211540/1‑140 (MQ=255)
                                                     cgcCCCCTTTTGCCCTTCCGCCAGGCTTTGTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGGAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCTAGGGATCGGTCGGGCCAc                                                                                     >  2:1127262/1‑140 (MQ=255)
                                                       cccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCg                                                                                   >  1:1379195/1‑140 (MQ=255)
                                                        ccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGa                                                                                  <  1:675762/140‑1 (MQ=255)
                                                                     tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAttt                                                                     >  1:792037/1‑140 (MQ=255)
                                                                     tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAttt                                                                     >  1:1851280/1‑140 (MQ=255)
                                                                     tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAttt                                                                     >  1:230576/1‑140 (MQ=255)
                                                                     tCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAAttt                                                                     >  2:453319/1‑140 (MQ=255)
                                                                      ccGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtc                                                                    <  1:600764/140‑1 (MQ=255)
                                                                      ccGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtc                                                                    <  2:1102090/140‑1 (MQ=255)
                                                                      ccGACAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtc                                                                    >  2:1693920/1‑140 (MQ=255)
                                                                           aGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGc                                                               >  1:431570/1‑140 (MQ=255)
                                                                             gCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                             <  1:1456151/140‑1 (MQ=255)
                                                                                ttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAACTGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                          <  2:260583/140‑1 (MQ=255)
                                                                                ttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                          <  2:641708/140‑1 (MQ=255)
                                                                                ttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                          <  2:357982/140‑1 (MQ=255)
                                                                                ttttGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                          <  2:1763900/140‑1 (MQ=255)
                                                                                   tGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACgg                                                       <  1:1114177/140‑1 (MQ=255)
                                                                                     aCGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAACTGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGTTTGAGTTTAATTTCTCGGCGTGAACggtg                                                     <  1:449777/140‑1 (MQ=255)
                                                                                     aCGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                     <  1:971889/140‑1 (MQ=255)
                                                                                     aCGCCAGCGTGCGTTGCACGATGCCCTGAATTGTACACCATTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                     <  2:1866607/140‑1 (MQ=255)
                                                                                       gCCAGCGTGCTTTGCACGTTGCCCTGATATGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCAGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                   <  1:1848308/140‑1 (MQ=255)
                                                                                       gCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCGAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAAGcgtggt                                                   >  2:375937/1‑140 (MQ=255)
                                                                                       gCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                   >  1:197826/1‑140 (MQ=255)
                                                                                       gCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                   >  1:475275/1‑140 (MQ=255)
                                                                                        ccAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGgtggtg                                                  >  1:1224023/1‑140 (MQ=255)
                                                                                         cAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAAACTGCAAGGGATCGGTCGGGCCACCGATAGAGTTTAATTTCTCGGCGAGAACGGTGGTGa                                                 >  2:971783/1‑140 (MQ=255)
                                                                                          aGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                <  1:1400276/140‑1 (MQ=255)
                                                                                          aGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                <  2:714771/140‑1 (MQ=255)
                                                                                          aGCGTGCGTTGCACGATGCCCTGATTTGAACACCAGTAAAAGGTGTCCACCAGAATCAAACCGTTAAATGCCACGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                <  2:1770057/140‑1 (MQ=255)
                                                                                            cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTACATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                              <  1:1081472/140‑1 (MQ=255)
                                                                                            cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                              <  1:1267357/140‑1 (MQ=255)
                                                                                            cgtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                              <  1:1267249/140‑1 (MQ=255)
                                                                                             gtgcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTg                                             <  2:631110/140‑1 (MQ=255)
                                                                                               gcgtTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCt                                           <  1:15714/140‑1 (MQ=255)
                                                                                                     cACGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGc                                     >  1:811142/1‑140 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:1135258/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:721186/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:1301909/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:1305722/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:1554803/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:188016/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCTCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  2:154405/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  2:1765684/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  2:899814/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:1025134/140‑1 (MQ=255)
                                                                                                      aCGATGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGcc                                    <  1:384372/140‑1 (MQ=255)
                                                                                                        gATGCCCTGATTTGTACACGCGTAAAATGTGATCACCAGAATCAAACCGGTAAATGCCTCGCCAATAGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGAGAGTTGCTCAATGCCCt                                  >  2:1588663/1‑140 (MQ=255)
                                                                                                        gATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCt                                  >  1:882153/1‑140 (MQ=255)
                                                                                                          tGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGc                                >  2:1843273/1‑140 (MQ=255)
                                                                                                          tGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGc                                <  2:506232/140‑1 (MQ=255)
                                                                                                          tGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGc                                >  2:3189/1‑140 (MQ=255)
                                                                                                          tGCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGc                                >  2:1071348/1‑140 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:1634959/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  1:1009044/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:767236/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:101507/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:1602462/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:1235561/140‑1 (MQ=255)
                                                                                                           gCCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                               <  2:345052/140‑1 (MQ=255)
                                                                                                               tGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATaaa                           >  2:599283/1‑140 (MQ=255)
                                                                                                               tGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATaaa                           >  1:225156/1‑140 (MQ=255)
                                                                                                                 aTTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGc                         <  1:152252/140‑1 (MQ=255)
                                                                                                                    tGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGc                      >  2:799608/1‑140 (MQ=255)
                                                                                                                      tACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGccc                    >  1:308784/1‑140 (MQ=255)
                                                                                                                      tACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGccc                    >  1:1779350/1‑140 (MQ=255)
                                                                                                                              tAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCAt            <  2:1306899/140‑1 (MQ=255)
                                                                                                                                   ggTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATAGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                >  1:700942/1‑99 (MQ=255)
                                                                                                                                   ggTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATAGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAg                                                <  2:700942/99‑1 (MQ=255)
                                                                                                                                     tgtTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt     <  2:414327/140‑1 (MQ=255)
                                                                                                                                     tgtTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt     <  2:1267434/140‑1 (MQ=255)
                                                                                                                                     tgtTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt     >  1:1724134/1‑140 (MQ=255)
                                                                                                                                     tgtTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGAtt     >  1:1600269/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                           
CCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGG  >  NC_000913/3859218‑3859493

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