Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,070 Δ1 bp intergenic (‑28/+31) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0690G.96.9% 296.0 / 2.3 64intergenic (‑27/+32)alsB/alsRD‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (22/40);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

AATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAAT  >  NC_000913/4311819‑4312315
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
cggacgattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAg                                                                                                                                                                                                                                                        <  2:368334‑M2/13‑1 (MQ=255)
 ggacgattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACagt                                                                                                                                                                                                                                                       >  1:57335‑M2/238‑250 (MQ=255)
 ggacgattttgccggtatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttgtttcatcggtctcgctccagttaagcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacagtttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACagt                                                                                                                                                                                                                                                       <  2:497441‑M2/14‑2 (MQ=255)
    cgatgttgccgttagcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGgacagtaac                                                                                                                                                                                                                                                    <  2:930‑M2/17‑5 (MQ=255)
     gattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACt                                                                                                                                                                                                                                                  >  2:276860‑M2/234‑251 (MQ=255)
      attttgccgttatcgacgcccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcattcgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTc                                                                                                                                                                                                                                                 <  1:637989‑M2/19‑1 (MQ=255)
      attttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTc                                                                                                                                                                                                                                                 <  2:638808‑M2/19‑1 (MQ=255)
        tttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTcc                                                                                                                                                                                                                                                <  1:148443‑M2/20‑1 (MQ=255)
        tttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaacgctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTcc                                                                                                                                                                                                                                                <  2:199658‑M2/20‑1 (MQ=255)
              gttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATgg                                                                                                                                                                                                                                         <  1:340633‑M2/27‑1 (MQ=255)
                tatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGc                                                                                                                                                                                                                                        <  2:12595‑M2/28‑1 (MQ=255)
                tatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCg                                                                                                                                                                                                                                       >  1:539379‑M2/223‑251 (MQ=255)
                 atcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCgg                                                                                                                                                                                                                                      >  1:479320‑M2/222‑251 (MQ=255)
                 atcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaacgctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCg                                                                                                                                                                                                                                       <  2:175199‑M2/29‑1 (MQ=255)
                                 attttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                                                                                                      <  1:6973‑M2/46‑1 (MQ=255)
                                     tgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                                   <  2:192717‑M2/49‑1 (MQ=255)
                                        aacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                                                                                                               >  2:423370‑M2/199‑251 (MQ=255)
                                                     ataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                                                                                                                  <  2:568822‑M2/66‑1 (MQ=255)
                                                        ggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgcctcattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                                                                                                                                                               <  1:132464‑M2/69‑1 (MQ=255)
                                                           gcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                                                                                                            <  2:180779‑M2/72‑1 (MQ=255)
                                                                ccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                                                                                                                                                                       <  1:593578‑M2/77‑1 (MQ=255)
                                                                  gtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                                                                                                                                                     <  1:123568‑M2/79‑1 (MQ=255)
                                                                    gacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                                                                                                                                                                   <  2:357995‑M2/81‑1 (MQ=255)
                                                                         cttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacagactgtcgtcatttttgacgtgaGGTAAGGGGGCAGTAACTCCGCGATGGCGGAAt                                                                                                                                                                                                                                   >  1:620460‑M2/166‑198 (MQ=255)
                                                                         ctttttgcattggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagtgttgacgtgaGGTAAGGGGACAGTAACACCGCGATGGCGGAAt                                                                                                                                                                                                                                   <  2:620460‑M2/33‑1 (MQ=255)
                                                                             tttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAGGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTGCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                                                                                                                                                                           <  2:467362‑M2/89‑1 (MQ=255)
                                                                                     gtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCata                                                                                                                                                                  <  1:37463‑M2/98‑1 (MQ=255)
                                                                                         cgctccagttaatcaaaacacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGTATTTCCTTTTTCTTATTTCATCGCGGCTTGTGAGAATAAATCAGCAATTGCGCCGGTTTTTTCCatatt                                                                                                                                                                <  1:204582‑M2/100‑1 (MQ=255)
                                                                                           ctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaat                                                                                                                                                             >  1:426566‑M2/148‑250 (MQ=255)
                                                                                                        aatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCa                                                                                                                                                >  1:568730‑M2/135‑250 (MQ=255)
                                                                                                                atacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGc                                                                                                                                       >  1:267482‑M2/127‑251 (MQ=255)
                                                                                                                     cgctctcgactacagtatgcatcttttatgccacattttttgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGc                                                                                                                                  >  1:547971‑M2/122‑251 (MQ=255)
                                                                                                                      gctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCg                                                                                                                                 >  2:152405‑M2/121‑251 (MQ=255)
                                                                                                                           cgactacagtatgcatcttttatgccccattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGAGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTACATGTTGAGATGGGCGAc                                                                                                                               >  2:546026‑M2/116‑248 (MQ=255)
                                                                                                                           cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACag                                                                                                                             <  2:479320‑M2/135‑1 (MQ=255)
                                                                                                                             actacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                                                                                                          <  1:361004‑M2/138‑1 (MQ=255)
                                                                                                                                       gcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                                                                                                <  2:618113‑M2/148‑1 (MQ=255)
                                                                                                                                             tttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTc                                                                                                             <  1:643219‑M2/151‑1 (MQ=255)
                                                                                                                                             tttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTc                                                                                                             <  2:57335‑M2/151‑1 (MQ=255)
                                                                                                                                             tttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGAGACAGAGACAAAACAAGCGtca                                                                                                            <  2:568730‑M2/152‑2 (MQ=255)
                                                                                                                                                   ccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                                                                                                                                                       >  1:612605‑M2/92‑136 (MQ=255)
                                                                                                                                                   ccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTGTTCtt                                                                                                                                                                                                                       <  2:612605‑M2/45‑1 (MQ=255)
                                                                                                                                                    cacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAggtt                                                                                                   >  1:45123‑M2/91‑251 (MQ=255)
                                                                                                                                                          ttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGt                                                                                             >  2:244214‑M2/85‑251 (MQ=255)
                                                                                                                                                            atgtggggtcattccctgatattacgggcactagttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTaa                                                                                           <  2:45123‑M2/169‑1 (MQ=255)
                                                                                                                                                                     cattccctgatattacgggcactatgtattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgggaGGTAAGGGGACAGTAACTCCGCGATGGCTGAATTACCTTTTTCTTATTTCAGCGCGCCTGGGGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCCGCAATGTTGGTTGTAATATTCTTGc                                                                                  <  2:381977‑M2/178‑1 (MQ=255)
                                                                                                                                                                        tccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCg                                                                                <  1:409437‑M2/180‑1 (MQ=255)
                                                                                                                                                                         ccctgatattacgggcactatttattccaaactctgacgaaaaacaggctgtcgttaggtgtgaagtgaGGGAAGGGGACAGTACCTCCGCGATGGCTGAATTACCTTTTTCATATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTCCCATATGAAGATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCGTGCCGAg                                                                              <  2:330450‑M2/182‑1 (MQ=255)
                                                                                                                                                                             gatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGc                                                                            <  2:547971‑M2/184‑1 (MQ=255)
                                                                                                                                                                             gatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCAt                                                                          >  1:450336‑M2/66‑251 (MQ=255)
                                                                                                                                                                                attacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATCAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTAc                                                                       >  1:590784‑M2/63‑251 (MQ=255)
                                                                                                                                                                                    cgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGAcc                                                                   >  2:330992‑M2/59‑251 (MQ=255)
                                                                                                                                                                                            atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAAc                                                            <  1:288329‑M2/200‑1 (MQ=255)
                                                                                                                                                                                            atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                           <  1:546026‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                    aaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTcc                                                   <  1:529805‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                        ctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAgg                                                                             <  2:459037‑M2/183‑1 (MQ=255)
                                                                                                                                                                                                        ctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAgg                                                                             >  1:459037‑M2/39‑221 (MQ=255)
                                                                                                                                                                                                           tgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTg                                            >  1:280256‑M2/36‑251 (MQ=255)
                                                                                                                                                                                                                           gtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGAAGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCTTTTCCGGGGCTGGTGAGCAAataatat                             <  2:450336‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                             cgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGcaaataatatgct                          >  2:45641‑M2/18‑248 (MQ=255)
                                                                                                                                                                                                                             cgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGcaaataatatgat                          <  1:152405‑M2/234‑4 (MQ=255)
                                                                                                                                                                                                                                cagttttgacgtgaGTTAAGGGGACAGTAACTCCGCGATGGCGTAATTACCTTTTTCTTATTTCAGCGCGACTGGTGAGAAGAAATCACCAATTGCGCCTGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGTGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCgg                       <  1:493950‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                       gacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTCATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGCGGCATTACGACCTAATAACGGCGTTTCCGCGGCTGGTGCGCAAATAATATGATCGGCCAGttt                >  2:130126‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                       gacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGtt                 <  2:267482‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                                       gacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGtt                 <  2:518314‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                                       gacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATAGGATCGGCCAGttt                >  1:487345‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                                   ggggACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTg     >  2:330869/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     ggACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAAt  <  2:280256/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
AATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAAT  >  NC_000913/4311819‑4312315

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