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.95.5% 191.2 / 2.6 45intergenic (‑27/+32)alsB/alsRD‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (18/24);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

ATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGT  >  NC_000913/4311820‑4312311
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
ggacgattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACagt                                                                                                                                                                                                                                                   <  1:170841‑M2/14‑2 (MQ=255)
 gacgattttgccgttatcgacgaccaggttgattttgcaacatgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACagt                                                                                                                                                                                                                                                   >  1:476573‑M2/237‑249 (MQ=255)
     attttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagtgttgacgtgaGGTAAGGGGACAGTAACTc                                                                                                                                                                                                                                             <  2:191545‑M2/19‑1 (MQ=255)
           ccgttatcgccgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGAt                                                                                                                                                                                                                                       >  2:176654‑M2/227‑251 (MQ=255)
           ccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGAt                                                                                                                                                                                                                                       <  1:152231‑M2/25‑1 (MQ=255)
               tatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGc                                                                                                                                                                                                                                    <  1:176654‑M2/28‑1 (MQ=255)
                atcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctcgcgtctacagtatgcatcttgtatgccacattttatgtggggtcattccctgatagtacgggcactatttattcaaaactctgacgaaaaacaggctgtcggcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCgg                                                                                                                                                                                                                                  <  1:485699‑M2/30‑1 (MQ=255)
                       accaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                                                                                                           <  1:196984‑M2/37‑1 (MQ=255)
                       accaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                                                                                                           <  1:270253‑M2/37‑1 (MQ=255)
                           ggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATTGCGGAATTACCtttt                                                                                                                                                                                                                       <  1:442417‑M2/41‑1 (MQ=255)
                              tgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                                                                                                                     <  1:365904‑M2/43‑1 (MQ=255)
                                   ttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                               >  2:319336‑M2/203‑251 (MQ=255)
                                   ttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattcccggatattacgggcactatttattcaaaactctgacgaaaaacaggctggcgtcagtttggacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                                                                                                                                                                <  2:361316‑M2/48‑1 (MQ=255)
                                    tgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGAGAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                               <  1:215302‑M2/49‑1 (MQ=255)
                                    tgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                               <  2:385009‑M2/49‑1 (MQ=255)
                                     gcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                                                                                                                                                              <  1:76233‑M2/50‑1 (MQ=255)
                                          ctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCtt                                                                                                                                                                                                                         <  1:488602‑M2/39‑1 (MQ=255)
                                          ctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCtt                                                                                                                                                                                                                         >  2:488602‑M2/196‑234 (MQ=255)
                                            gatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                                                                                                                                                                      <  1:319336‑M2/58‑1 (MQ=255)
                                                                  tgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                                                                                                                                                                >  2:23971‑M2/172‑251 (MQ=255)
                                                                   gacgtctcttttcatcggtctcgctccagttagtccaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttcttcaaaactctgacgaaaaacaggctgtcgtcagttttggcgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTATTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                                                                                                                                                               >  2:452043‑M2/171‑251 (MQ=255)
                                                                               catcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                                                                                                                                                                    <  1:529729‑M2/92‑1 (MQ=255)
                                                                                    gtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCata                                                                                                                                                              <  1:286838‑M2/98‑1 (MQ=255)
                                                                                     tctcgctccatgtaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTCTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAGCAACAATTGCGCcggttttttcca                                                                                                                                                                <  2:304092‑M2/96‑4 (MQ=255)
                                                                                     tctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGACATCAACAATTGCGCCGGTTTTTTCCatat                                                                                                                                                             >  2:455156‑M2/153‑251 (MQ=255)
                                                                                     tctcgctccagttaatcaaatcacgcatacgcgctctcgaatacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTCTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCACCAATTGCGCcggttttttcca                                                                                                                                                                >  1:304092‑M2/153‑245 (MQ=255)
                                                                                                  aatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTg                                                                                                                                                <  1:23971‑M2/112‑1 (MQ=255)
                                                                                                      aaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTc                                                                                                                                             >  2:414448‑M2/136‑250 (MQ=255)
                                                                                                        atcacgcatacgcgctctcgactacagtatgcatcttttatgccacagtttatgtgggttcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTGGTCCATATTAATATTAGCGTGTTCAAt                                                                                                                                          <  2:438438‑M2/118‑1 (MQ=255)
                                                                                                            cgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattacctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATTGCGGCATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGAGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTg                                                                                                                                      >  2:457116‑M2/130‑251 (MQ=255)
                                                                                                              catacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTg                                                                                                                                      >  1:368730‑M2/128‑249 (MQ=255)
                                                                                                                                    atgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaa                                                                                                               >  1:52136‑M2/106‑250 (MQ=255)
                                                                                                                                                         ttatgtggggtcattccctggttttacgggcactatttaatcaaaactctgacgaaaaacaggctgtcgtcagttttggcgtgaGGTAAGGGGACCGTAACTCCGCGATGGCGGAATTGCCTTTTTCTTATTTCAGCGCGCCTGGTTAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGGTGTGCTGGGCGACAGAGACAAAAAAAGCGTCCagcaaggttggatg                                                                                          <  2:132331‑M2/166‑6 (MQ=255)
                                                                                                                                                                       tccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCg                                                                            >  1:128771‑M2/71‑250 (MQ=255)
                                                                                                                                                                               attacgggctctatttattcaaacctctcacgaaaaacaggctgtcgccagttttgacgtgaGGTAAGGGGACAGAAACTCAGCGATGGCGGAATTACCTTTTTATTATTTCAGCGCGCCTGGAGCGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAATTAc                                                                   >  2:315018‑M2/63‑251 (MQ=255)
                                                                                                                                                                               attacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGtt                                                                                                                                       >  1:164992‑M2/63‑183 (MQ=255)
                                                                                                                                                                               attacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGtt                                                                                                                                       <  2:164992‑M2/121‑1 (MQ=255)
                                                                                                                                                                                  acgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGa                                                                 <  1:455156‑M2/191‑1 (MQ=255)
                                                                                                                                                                                             ttattcaaaactctgacgaaaaacaggatgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGc                                                     >  1:180308‑M2/49‑251 (MQ=255)
                                                                                                                                                                                                           gacgaaaaacaggctgtcgtcagttgtgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCGTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCt                                         <  1:452043‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                               aaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGttt                                                 <  2:440610‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                               aaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGttt                                                 >  1:440610‑M2/31‑237 (MQ=255)
                                                                                                                                                                                                                  aacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCa                                 >  2:361947‑M2/28‑250 (MQ=255)
                                                                                                                                                                                                                                agttttgacgagcGGTAAGGGGCCAGTCACTCCGCGATGGCGGAATTACATTTTTCTTATTTCCGCGAGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc                  >  2:478532‑M2/14‑251 (MQ=255)
                                                                                                                                                                                                                                       acgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc           <  1:189860‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGTACGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTATCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGCGAGCAAATAATATGATCGGCCAGTTTCGCTAt      >  2:16077/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aaGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGAGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTACTAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGt  >  1:513696/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
ATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGT  >  NC_000913/4311820‑4312311

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