Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,536 Δ1 bp 5.4% coding (8/900 nt) yfeX ← porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,5360T. 5.4% 224.5 / 3.1 73coding (8/900 nt)yfeXporphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base . (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

ATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTT  >  NC_000913/2550310‑2550748
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
atgttgcAGCTCTTCCGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑g                                                                                                                                                                                                                     <  2:327816/222‑1 (MQ=255)
acctttcctccctcctGCCCCACCGCCGCCGCGCAGAGCGCGCCAGGGGTTTTTCCCAAAAGCAACCCCCGGTCCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCCGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑g                                                                                                                                                                                                                     <  1:575773/214‑1 (MQ=255)
    ttCAGCTCTTCCGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGAc                                                                                                                                                                                                                 <  2:441541/228‑1 (MQ=255)
           cttctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTtc                                                                                                                                                                                                          <  2:87474/228‑1 (MQ=255)
             tctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCCCTCTAAACCT‑GAGACATTGTtcct                                                                                                                                                                                                        >  2:556133/1‑228 (MQ=255)
              cagCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTAACAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTtcctc                                                                                                                                                                                                       <  2:274217/226‑1 (MQ=255)
                   ccAACGCCGCCGCTCAGAGCGGTCCAGGTGTTGGTACCAAAGGCAACCAGCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:446649/228‑1 (MQ=255)
                   ccAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCCCCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:232220/228‑1 (MQ=255)
                   ccAACGCCCCCGCTCAGGGCGCGCCAGGTTTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:575917/228‑1 (MQ=255)
                   ccAACCCCTCCGCGCAGATCTCGCCAGGTTTCGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAGGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:668014/228‑1 (MQ=255)
                      aCGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGCCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAaat                                                                                                                                                                                               >  2:298204/1‑228 (MQ=255)
                        cccaccggccgggcggcgcGAGGTGTTTTTACCAAAGGCTACCCACGACCCAAGATGCGCGTCCGGGAATTTCCCTTCAAAAGTTGCCAGGTCGTCGGCAAATGTTTTACTGGCCGCCCGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAaataa                                                                                                                                                                                             <  1:318076/213‑1 (MQ=255)
                                  cgggcgagcCGGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACg                                                                                                                                                                                   <  2:248930/225‑1 (MQ=255)
                                  agggggcgcCGGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAGCCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACg                                                                                                                                                                                   <  1:429482/223‑1 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCCCGCGGCAATGTTCTGGCAAAATGCCACCCTGAACCC‑GAGCCATTTTTCCTCCTGAAAATAAAAATGcccccca                                                                                                                                                                                 >  1:387333/1‑223 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATACCCGCGCGGCAATTTTCTGGCAAAATGCCACTCTGAACC‑‑CAGACCTTGTTCCCCCTGAAAATATAAATTCCAcccca                                                                                                                                                                                >  1:591206/1‑224 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTAACGTTGGCTTCGATCCAAATCGCCCCGCGGCAATGTTCTGGCAACAGGCCACCCTGAACC‑‑GAGACATTGTTCCCCCTGAAAATAATAATGCCCCGCaa                                                                                                                                                                                >  1:169428/1‑228 (MQ=255)
                                       acgcCAGGTGTTGGTACCAAAGGCAACCAACGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTTGCGGCAAATGTTTTACGGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATc                                                                                                                                                                              <  2:405942/227‑1 (MQ=255)
                                         gcCAGGTTTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCg                                                                                                                                                                            <  2:480227/228‑1 (MQ=255)
                                                gttgttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCa                                                                                                                                                                       <  2:582189/226‑1 (MQ=255)
                                                   gttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgtt                                                                                                                                                                  >  2:115224/1‑228 (MQ=255)
                                                   gttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAAGCCCACCCAATCCGTGGCAttgtt                                                                                                                                                                  >  2:460874/1‑228 (MQ=255)
                                                    ttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgtt                                                                                                                                                                  <  2:670604/227‑1 (MQ=255)
                                                       ccAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGt                                                                                                                                                                           >  1:39594/1‑215 (MQ=255)
                                                        cAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTg                                                                                                                                                                          <  2:39594/215‑1 (MQ=255)
                                                               aaCCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAtttt                                                                                                                                                      >  2:274622/1‑228 (MQ=255)
                                                                accaccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGGGGCATGTTTGAGCGCTCtttt                                                                                                                                                      >  2:136695/1‑227 (MQ=255)
                                                                       cACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAACTGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGCAATGGTTGAGCGCTATTTTACGTGgtt                                                                                                                                              >  2:334352/1‑228 (MQ=255)
                                                                                      ccGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACCCTGAACCT‑GAGACATTGTTCCTCCTGAAAATATTAATGCCACGCAATCCGTGCCATTTTTGAGCGCTATTTTACGTGTTTTTTAGCGCAGGCGGt                                                                                                                               >  2:691006/1‑228 (MQ=255)
                                                                                         ggAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGTCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAACTAATAATGCCACGCAATCCGTGGAATGGTTTCGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtgtg                                                                                                                            >  2:330586/1‑228 (MQ=255)
                                                                                                               gtgtttaaggtacagtGTTTTCCTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCAAGTTGGCTGCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccg                                                                                                       <  2:680952/212‑1 (MQ=255)
                                                                                                               cAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccgc                                                                                                      <  2:387333/228‑1 (MQ=255)
                                                                                                                 gTTTGGCGGCAAATGTTTTACGGGCCGTCCGCAGGGCGTCAACTTCCCCTTTCCCGTTGGCCTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCa                                                                                                    <  2:680423/228‑1 (MQ=255)
                                                                                                                  tttGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGCCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGGTTTTTCAACGGCAAATTaaaccccccaa                                                                                                   >  1:45652/1‑219 (MQ=255)
                                                                                                                   ttGTCCGTAAAGTTTTTTCCGGCCGCACGCGGGGGGTGAACTTCCCCGTTCCCGTTGGCTCCGGACCAAATCGCCGCGCGGCAATGTGCTGGCAAAATGCCCCTCTGATCCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                                                                                                  <  1:361278/228‑1 (MQ=255)
                                                                                                                   ttGGCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                                                                                                  <  2:266132/228‑1 (MQ=255)
                                                                                                                       cGGCAAATTTTTTACTGGCCGCACGCGGGGCGCCAACTTCCCCTTGCACTGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATa                                                                                                <  1:115224/226‑1 (MQ=255)
                                                                                                                        ggCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTCTTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCGGATAAt                                                                                              >  2:427200/1‑227 (MQ=255)
                                                                                                                            aaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GCGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGGGTTTTTTCGGGGCGGCGGTTTGTTCGGAGGCAAATTAACGCCCCCAGATAAttttcc                                                                                         >  2:502971/1‑226 (MQ=255)
                                                                                                                               gTTTTGCTGGGCGCACGCAGGGCGTCAAATGCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGa                                                                                      <  2:563848/228‑1 (MQ=255)
                                                                                                                                        gCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTTGCTTCGATCCAAATCGCCGCGCGGCACTGTTCTTGCCAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCACTCCGTGGCATTGTTGAGCGCTCTTTTCAGTGTTTTTTACCGGCGGGTGTTTGTTTCGGCGCAATATTACCCCGCACGTTCATTTTACTGACTTTcgctc                                                                             >  1:330735/1‑224 (MQ=255)
                                                                                                                                        gCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCCCCAGATAATTTTACTGACTTTCCAGt                                                                             >  2:154722/1‑228 (MQ=255)
                                                                                                                                         ccGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAACTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGCCGGCGGTTTGTTCAGCCTCAAATTAACGCCTCCCGATAAtttctttctcttttatac                                                                             >  1:680952/1‑212 (MQ=255)
                                                                                                                                         ccGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGt                                                                             >  2:228635/1‑227 (MQ=255)
                                                                                                                                          cgcacgcgGGGCGTCAAGTTCCCCGTTCACCTTGGCTTCGAACCAAATCCGCGCGGGGCAATGGTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGGGGCATTGTTGAGCGCTATTTGACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGACGCCAGAGAATTTTACTGACGTTCCAGttt                                                                           <  1:671523/228‑1 (MQ=255)
                                                                                                                                          cgcacgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GCGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCCTTGTTTAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGGCGCAAATTAACGCCGCCAGATAATTTTTTTGAGTTTCCCGt                                                                             >  1:634605/1‑226 (MQ=255)
                                                                                                                                                  gggCGTCCACTTCCCCTTTCACGTTGGCTTAGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GACAAATTGTTACTCCTGCAACTAATAATGCGTCGCAATCCGTGGCTTTGTTGAGGGCTATTTTACGTGTTTTTTCGCGGCGGCGGGTTTTTCATACGCAAATTAAACCCTCCCTACAATATTCCTGGCTTTCCTGtttttttaggg                                                                   >  1:592826/1‑222 (MQ=255)
                                                                                                                                                  gggCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGGTTTTCAGGGt                                                                   >  2:117647/1‑228 (MQ=255)
                                                                                                                                                   gtggTCAATTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGGGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGGCATTGTTCCTGCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTa                                                                  <  1:274622/225‑1 (MQ=255)
                                                                                                                                                   ggCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGAAAATTAACGCCCCCCAATAATTTTACTGGCTTTCCAGGTTTTAAgggaa                                                                  >  1:446195/1‑226 (MQ=255)
                                                                                                                                                         aaCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGCTTGTTCAGACGCAAATTAACGCCGCCAGATATTTTTACTGACTTTCCAGTTTTTCAGGGAATCggc                                                             >  2:602146/1‑225 (MQ=255)
                                                                                                                                                          aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACC‑‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:591206/228‑1 (MQ=255)
                                                                                                                                                          aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACC‑‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:169428/228‑1 (MQ=255)
                                                                                                                                                            ttCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCCGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:676054/227‑1 (MQ=255)
                                                                                                                                                             tCCCATTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTACTACTGAAAATAATAATGCCACGCAAACCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAAACGCACAATAACGCAGCAAGATAGTTTTATTGGCTTTCCAGTTTTTCAGGGCATCGTTCGAAg                                                        >  2:87592/1‑228 (MQ=255)
                                                                                                                                                                   ttCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCAtt                                                   >  1:625079/1‑227 (MQ=255)
                                                                                                                                                                    ttccGTTGGCTGCGATCCAAATCGCCACGCGGCAATTTTCGGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTTGCATTGTTGAGCTCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTGTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:556133/225‑1 (MQ=255)
                                                                                                                                                                    tctcGTTGGCTGCTAACCAAATCGCCGCGCGGCAATGTTCTGGCAAAAGGCCACTCTGACCCT‑GAGACATTGTTCCTCCTGAAAAGAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:502971/225‑1 (MQ=255)
                                                                                                                                                                    tcccGTTGGCGTCGATCCAAATCGCCGCGCGGCAATGTTCTCGAAAAAGGCCACTCTGAACCT‑TGGACATTGTTCCTCCTGAAAATAATTATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:136695/225‑1 (MQ=255)
                                                                                                                                                                         ttGGATTCGATCCAAATCGCCGTGCTGCAAGTTTCTGTCAAAATGCCACCCTGAACCT‑GAGCCATTGTTCCTCCTGAAAATAATAATTCCACGGAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCTGTTTGTCCAGACGCAAATTAACGCCGCCAGATAATTTTACTTACTTTCCAGTTTGTCAGGGTATCGTCCGAAGGCATTAACCCtt                                            <  1:688260/228‑1 (MQ=255)
                                                                                                                                                                                  aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTAGGTGTTTTTTCGCGGCGGCGGGTTGTTTTGACGCAAATTAACGCCGCCAGATCATTTTACTGGCTTTCCAGGTTTTCAGGTTACCTCCCTAAGGGATTAACCCtaccatccca                                    >  1:677228/1‑218 (MQ=255)
                                                                                                                                                                                  aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCCAAGGCATTAACCCTTCAGGCCCa                                    >  2:621938/1‑227 (MQ=255)
                                                                                                                                                                                                    cAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACt                 >  2:124204/1‑228 (MQ=255)
                                                                                                                                                                                                     aaTGTTCGGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGaa                   <  1:124204/225‑1 (MQ=255)
                                                                                                                                                                                                       ttttCTTGCAAAATTCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:253505/226‑1 (MQ=255)
                                                                                                                                                                                                       tGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:625079/228‑1 (MQ=255)
                                                                                                                                                                                                       tGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:446195/228‑1 (MQ=255)
                                                                                                                                                                                                       tGTTATGGAAAAATGCCCATCTGAACTT‑TGGACTTTTGTCCTCCTGAAAATAATTATGCCTCGCAATCCGTGTCATTGGTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGGCGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  1:691006/228‑1 (MQ=255)
                                                                                                                                                                                                           ttGGGAAAATGCCGCCCGGAACCTGGAGA‑TTTGTTCCTCCTAAAAATCATAATGCCGCGCAATCCGTGGCATTATTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGTGGTTTGTTCAGACGCAAATTAACGCCGCCAAATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCt          <  1:125814/227‑1 (MQ=255)
                                                                                                                                                                                                              gCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACTTGTTTTTTCGCGCCGGAGTTTTGTTCAGACGCAAATTAACGCCGCCAGTTAATTTTACTGACTTTCCAGTTTTTCAGTGAATCGTCCGCAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGCACTg                >  2:410426/1‑219 (MQ=255)
                                                                                                                                                                                                              gCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCATGGCATCGTCCGAAGGCTTTAACCCTTCCGGACCACGCCATTCCCCCCTCAACgggatgc           >  2:125814/1‑220 (MQ=255)
                                                                                                                                                                                                                   aTGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGtt  >  2:166268/1‑228 (MQ=255)
                                                                                                                                                                                                                   aTGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGtt  >  2:227780/1‑228 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
ATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCT‑GAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTT  >  NC_000913/2550310‑2550748

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