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.98.4% 278.0 / ‑3.8 61intergenic (‑27/+32)alsB/alsRD‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (26/34);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTG  >  NC_000913/4311824‑4312312
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
gattttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACt                                                                                                                                                                                                                                               <  1:469790‑M2/18‑1 (MQ=255)
     tgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTCAGGGGACA‑GTAACTCcgcg                                                                                                                                                                                                                                          >  1:679027‑M2/229‑251 (MQ=255)
      gccgttatcgacgaccagtttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaacacaggctgacgtaatttttgacgtgaGGTAAGGGGACA‑GTAACCCcgcg                                                                                                                                                                                                                                          >  2:594202‑M2/228‑250 (MQ=255)
         gttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATgg                                                                                                                                                                                                                                      <  2:626970‑M2/27‑1 (MQ=255)
                acgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgctattacgggcactatttattcaaaactctgacgaaacacatgctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAAtt                                                                                                                                                                                                                               >  1:58001‑M2/218‑251 (MQ=255)
                          tgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacactttatgtggggtcattccctgatattacggacactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTACGCGATGGCGGAATTACCTTTTTCt                                                                                                                                                                                                                     >  2:147395‑M2/208‑251 (MQ=255)
                           gattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatctttgatgccacatttgatgtggggtcattccctgatattacgggtactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                                                                                                                                                     <  1:618546‑M2/44‑1 (MQ=255)
                                  caacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                                                                                                                                                               >  1:323208‑M2/200‑249 (MQ=255)
                                    acctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                                                                                                                                                           >  1:674115‑M2/198‑251 (MQ=255)
                                      ctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                                                                                                                                                                         >  1:271645‑M2/196‑251 (MQ=255)
                                            cgcaataggggcaaaccgtgacgacttttttcagcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCtggc                                                                                                                                                                                                   <  1:117114‑M2/62‑2 (MQ=255)
                                                ataggggcaaaccgtgccgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaacaacaggctgtcgtcagttttgacgtgcGGTAAGGGGACA‑GTAACTCCGAGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCctgggg                                                                                                                                                                                                  >  2:683135‑M2/186‑246 (MQ=255)
                                                  cggggcaaaccgtaacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtccgttttgacgtgaGGGAAGGGGACA‑GTACCTCCGCGATGGCGGTATTACCTTTTTCTTATTTCAGCGCTCCTGgggagca                                                                                                                                                                                              >  2:220086‑M2/184‑248 (MQ=255)
                                                          accgtgacgacttttttcatcggtctcgctccaggtaatcaaatcacgaatacgcgctctcgactacagtatgcatcttttatgacacattttatgtggggtaattccctgacattacgggcactatttattcaaaactctgccgaaaaacaggctgtcgtcagttttgacgtgaGCTCAGGGGACA‑GTAACTCCGCGATGGCGGTATTACCTTTTTCTNATTTCAGCGCGCCtggtgagaagaactccac                                                                                                                                                                                     >  2:365421‑M2/176‑245 (MQ=255)
                                                             gtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                                                                                                                                                  >  1:22545‑M2/173‑251 (MQ=255)
                                                               gatggcttttgtcatctgtcttgctccagttaatcaaatagcgaatatgcgctctcgactacagtatgcatcttttatgccacagtttatgtggggtcattccctgatattacgggcactatttatgcaaaactatgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTAGTTCAGCGCGCCTGGTGAGAAGAAATCACCAAtt                                                                                                                                                                                 <  2:684654‑M2/80‑1 (MQ=255)
                                                                acgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                                                                                                               <  1:673975‑M2/82‑1 (MQ=255)
                                                                  gacttttttcatcgggctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                                                                                                                                                             <  2:387378‑M2/84‑1 (MQ=255)
                                                                   acttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                                                                                                                                                            <  1:473073‑M2/85‑1 (MQ=255)
                                                                       ttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                                                                                                                                                         <  2:244314‑M2/88‑1 (MQ=255)
                                                                         ttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttt                                                                                                                                                                      >  1:477165‑M2/161‑251 (MQ=255)
                                                                                gtttcgctccagttaatcaaatctcgcatacgcgctctcgactacaggatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttagtcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTACCTCCGCGATGGCGGAATTACCGTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATGGCGCCGGTTTTTTCCata                                                                                                                                                               <  2:3‑M2/98‑1 (MQ=255)
                                                                                    cgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaa                                                                                                                                                           <  1:536262‑M2/102‑1 (MQ=255)
                                                                                           gttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGc                                                                                                                                                    >  1:29506‑M2/143‑251 (MQ=255)
                                                                                             taatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCt                                                                                                                                                   <  1:715458‑M2/110‑1 (MQ=255)
                                                                                                 caaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTc                                                                                                                                              >  1:426203‑M2/137‑251 (MQ=255)
                                                                                                  aaatcacgcatacgcgctctcgactacagtattcatcttttatgccacagtttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCa                                                                                                                                             <  1:545088‑M2/116‑1 (MQ=255)
                                                                                                         gcatacgcgctatcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaacgctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGa                                                                                                                                      <  1:425843‑M2/123‑1 (MQ=255)
                                                                                                             acgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCt                                                                                                                                   <  2:249144‑M2/126‑1 (MQ=255)
                                                                                                                              gtatgaatcttttatgcaccattttctgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaccaggctgtcgtcagttttgacgtgaGGCAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaa                                                                                                                  >  2:607727‑M2/108‑250 (MQ=255)
                                                                                                                                    atcttttatgccacattttatgtggggtcattccctgatattacgggcactaattattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGCTGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAAAGAGACAAAAAACGCGt                                                                                                           >  2:695879‑M2/102‑251 (MQ=255)
                                                                                                                                      cttttatgcaacattttatggggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTcc                                                                                                         >  2:700577‑M2/100‑251 (MQ=255)
                                                                                                                                        gttatgccacatgttatggggggtcattcncggatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACCCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTcc                                                                                                         <  2:360926‑M2/152‑1 (MQ=255)
                                                                                                                                              ccacattttatggggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCagcaaggt                                                                                                 <  1:147395‑M2/160‑4 (MQ=255)
                                                                                                                                                       atgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggatgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTa                                                                                         >  1:618162‑M2/83‑250 (MQ=255)
                                                                                                                                                              gtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCtt                                                                                 >  1:461891‑M2/76‑251 (MQ=255)
                                                                                                                                                                        gatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCa                                                                        <  1:74116‑M2/185‑1 (MQ=255)
                                                                                                                                                                        gatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCAt                                                                       <  2:426203‑M2/186‑1 (MQ=255)
                                                                                                                                                                           attacgggctctatttatgcaatactctgacgaaaaaccggctggcttcagttttgtcttgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTTCCTTTTTCTTATTTCAGCGCGGCTGGTGAGAAGAAATCAACAATTGCGCCGGTGTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGGGACAGAGACAAAAAATGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTAc                                                                    <  2:449836‑M2/189‑1 (MQ=255)
                                                                                                                                                                                 tgcactatttattcaacactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGATGGTTGTAATATTCTTGCCGATGCATTACGACCt                                                               >  1:478067‑M2/57‑250 (MQ=255)
                                                                                                                                                                                 ggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCt                                                               >  1:478076‑M2/57‑250 (MQ=255)
                                                                                                                                                                                   cacttgttattcaaaacgcggacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCtaat                                                            <  2:674115‑M2/197‑1 (MQ=255)
                                                                                                                                                                                    actatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                                                                                                                                                       >  1:297737‑M2/54‑143 (MQ=255)
                                                                                                                                                                                    actatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                                                                                                                                                       <  2:297737‑M2/90‑1 (MQ=255)
                                                                                                                                                                                     ctatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCtaata                                                           >  1:133314‑M2/53‑250 (MQ=255)
                                                                                                                                                                                         ttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGc                                                      >  1:550641‑M2/49‑251 (MQ=255)
                                                                                                                                                                                          tattcaaaactcttacgaaaaacaggctgtcgtcagttttgacttgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTATTCTTATTTCAGCGCGCCTGGTTAGAAGAAAGCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGCCATTACTACCTAATAACGGc                                                      <  2:617609‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                  actctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCggg                                             <  1:594202‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                     ctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCt                                          >  1:658763‑M2/37‑251 (MQ=255)
                                                                                                                                                                                                       gacgaaaaacaggctgtcgtcagttttgacgtgaGGGAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGc                                           <  1:633260‑M2/214‑1 (MQ=255)
                                                                                                                                                                                                       gacgaaaaacaggctgtcgtcagttgtgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTgg                                        <  2:292065‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                        acgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCtggg                                       <  2:461891‑M2/218‑2 (MQ=255)
                                                                                                                                                                                                          gaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTTCCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGa                                     <  2:656856‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                            aaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGc                                   <  1:563572‑M2/222‑1 (MQ=255)
                                                                                                                                                                                                                 aggctgtcgtcagttttgtcgttaGGTAAGGGTACA‑GTAACGCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCGTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaata                              <  2:679027‑M2/227‑1 (MQ=255)
                                                                                                                                                                                                                         gtcagttttgacgtgaGGTAAGGGGACA‑GTANCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATc                      <  2:323208‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                         gtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATc                      <  1:606393‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                           cagttttgacgtggGGTAAGGGGGCA‑GTAACTCCGCGATGGAGGAATTATCTTGTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGATGCATTACGACCTAATAACGGCGTTTCCGGGGGTGGTGAGCAAATAATATGATCg                     <  2:68849‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                            agttttgccgtgaGGTAAGGGGACA‑GTAAATCCGCGATGGCGGACTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGCAATCAACAATTGCGCCGGGTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGCGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATCTTCTTGCCGAGGCCTTCCGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc                   >  2:602281‑M2/14‑251 (MQ=255)
                                                                                                                                                                                                                                     gtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGTCTGGTGAGCAAATAATATGATCGGCCAGTTTCGc          <  2:271645‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                                             aGGGGACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTg  <  2:194773/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
CTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTG  >  NC_000913/4311824‑4312312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑