Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,219 A→G L258P (CTA→CCA)  alsR ← DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,2190AG98.3% 214.5 / NA 60L258P (CTA→CCA) alsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/33);  minor base T (0/1);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATG  >  NC_000913/4311975‑4312458
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
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)
                              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)
                                                                                                        cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCta                                                                                                                                                                                                                <  2:111120/174‑1 (MQ=255)
                                                                                                        cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCta                                                                                                                                                                                                                >  1:111120/1‑174 (MQ=255)
                                                                                                        cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGtt                                                                                                                                   <  1:700577/251‑1 (MQ=255)
                                                                                                          aTGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAGGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCCTGTTCAATGTTTAGCTTGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGATGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGACTCGCTATCGGGTAATGGGAGCTATGGGTta                                                                                                                                  >  2:176179/1‑250 (MQ=255)
                                                                                                                    tACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGCCAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAAc                                                                                                                                                                                                           >  2:158665/1‑167 (MQ=255)
                                                                                                                    tACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAAc                                                                                                                                                                                                           <  1:158665/167‑1 (MQ=255)
                                                                                                                     aCCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCt                                                                                                                      >  1:243582/1‑251 (MQ=255)
                                                                                                                        tttttCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTg                                                                                                                   >  1:331775/1‑251 (MQ=255)
                                                                                                                         ttttCTTATTTCAGCAAGCCTGGTGAGAAAAAATCAACAATTGCGCTGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGACCTGGGCGACAGAGACAAAAAAAGCGACCAGCCAGGTTGGCTGTAATATTATTGCCGAGGCATTACGACCTAATAACGGCGCTTCCGGGGCTGGTGAGGGAATAATATGATCGGCCAGTTTCGCTATCGTGGAATGGTTGCACTGGGTTATACAAATAATCTTTGc                                                                                                                  >  2:86789/1‑251 (MQ=255)
                                                                                                                               tATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGcccc                                                                                                               <  2:477165/248‑1 (MQ=255)
                                                                                                                                  ttCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAAAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACATAATAACGGCGTTTCCGGGGCTGGGGAGCAAATAATATGATCGGCCAGTTTCGCTATCGCAGAATGGTCGCTATGGTTTATACAAATACTCTTTGCCCCGTTCtt                                                                                                         >  2:285952/1‑251 (MQ=255)
                                                                                                                                    cAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtttt                                                                                                       <  1:404224/251‑1 (MQ=255)
                                                                                                                                         cgcCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGcc                                                                                                  >  1:250840/1‑251 (MQ=255)
                                                                                                                                           ccTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAg                                                                                                <  2:22545/251‑1 (MQ=255)
                                                                                                                                              ggTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCATTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTc                                                                                             >  2:512667/1‑251 (MQ=255)
                                                                                                                                                 gagaAGAAATCAACAATTGCGCCGGTTTTTTCCATATTTATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTa                                                                                           <  1:365421/250‑1 (MQ=255)
                                                                                                                                                          acaacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCttt                                                                                 <  1:229326/250‑1 (MQ=255)
                                                                                                                                                           caacaaTTGCGCCGGTTTTTTCCATATTAATAGTAGCTTGTTCAATTTTGAGCTGGGCGACAGAGACAAAAAAAGCTTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTTAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTGCTTTTTTGCCAGTTCTACGGCCGCttt                                                                                 <  2:478076/250‑1 (MQ=255)
                                                                                                                                                           aaacaaGTGCGCCGGTTTGTTCAATATTAATATTAGATTGTTCAATGTTGAGCTGGGCGACAGGGACAAAAAAAGCGTCAAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGTCCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCATTTTCGCTATCGTTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCACGTTATTTGTTGCCAGTGCTACGCCCGCttt                                                                                 <  2:478067/249‑1 (MQ=255)
                                                                                                                                                                   gcgcCGGTTTTTTCCATATTAATATTAGCTTGTTCCATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGGGAATGGTAGCTATGGGTTATACAATTAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATTACt                                                                        >  2:128969/1‑251 (MQ=255)
                                                                                                                                                                   gcgcCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACt                                                                        >  2:278481/1‑251 (MQ=255)
                                                                                                                                                                       cGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTGGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGtt                                                                    <  1:469373/251‑1 (MQ=255)
                                                                                                                                                                              taccATATTAATATTAGCGTTTTCAATGTTGAGCTGTGCGACAGAGACAAAAAAAGCGTCCAGCAAGGATGGTTGTAATATTCTTGCCGATGCAGTACGACATATTAACGGCGTTTCCGGGGGTGGTGGGCTAATAAGATGAACGTCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATACTCTTTGCCCCGTTCTTTTTTGCCAGGTCTACGGCCGCTTTTTCATCTCTGGTTCGCCCg                                                              <  2:414343/247‑1 (MQ=255)
                                                                                                                                                                                ccATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCGAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGa                                                            <  2:550641/249‑1 (MQ=255)
                                                                                                                                                                                                             aGCTGGGCGCCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCCAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATCCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCCCAACATCtccttc                              >  2:423668/1‑251 (MQ=255)
                                                                                                                                                                                                                            ccAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAAGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTCCCAGCACAACATCTCCTTCCTGTAACAACGAAgg               >  2:374463/2‑250 (MQ=255)
                                                                                                                                                                                                                            aCAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGc               <  1:683135/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaagGCGTCAAGCAAGGTTTGTTGTAATATTCTGGCCGAGGCATTACGACCTAATAACGGCGTTTACGGGGCTTGTGAGCAAATAATATGATCGGCCAGTTTCGCTAGCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTTCGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCAGTACCAGCACAACATCTCCTTCCTGTAACAACGCAGCGGa            <  2:545200/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAc           <  1:107767/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAc           >  1:662023/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCAGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAAAAACGAAGCGGAcat         >  1:619406/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaGCGTCCAGCAAGGTTGGGTGTAATATTCTTGCCGAGCCATTACGACCTAATAAAGGCGTTTCCGGTGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAAAAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCAATGGTTCGCCCGGAATGGGTCACTACCCGCACAACATCTCCTTCCTGTAACAACGAAGAGGAcat         >  2:466158/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcatca    <  1:607727/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcatcat   >  1:461131/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ccAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGTGCTGGTGAGCAAATAATAGGCTCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAACTAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGATCGCCCGGAATGGGTCACTACCAGCACCACATCTCCTTCCTGTAACAACGAAGCGGAcatctacatg  >  2:13547/1‑245 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATG  >  NC_000913/4311975‑4312458

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