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

ATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTT  >  NC_000913/2550310‑2550748
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
atgttgcAGCTCTTCCGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTg                                                                                                                                                                                                                     <  2:327816/222‑1 (MQ=255)
acctttcctccctcctGCCCCACCGCCGCCGCGCAGAGCGCGCCAGGGGTTTTTCCCAAAAGCAACCCCCGGTCCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCCGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTg                                                                                                                                                                                                                     <  1:575773/214‑1 (MQ=255)
    ttCAGCTCTTCCGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGAc                                                                                                                                                                                                                 <  2:441541/228‑1 (MQ=255)
           cttctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTtc                                                                                                                                                                                                          <  2:87474/228‑1 (MQ=255)
             tctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCCCTCTAAACCTGAGACATTGTtcct                                                                                                                                                                                                        >  2:556133/1‑228 (MQ=255)
              cagCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTAACAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTtcctc                                                                                                                                                                                                       <  2:274217/226‑1 (MQ=255)
                   ccAACGCCGCCGCTCAGAGCGGTCCAGGTGTTGGTACCAAAGGCAACCAGCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:446649/228‑1 (MQ=255)
                   ccAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCCCCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:232220/228‑1 (MQ=255)
                   ccAACGCCCCCGCTCAGGGCGCGCCAGGTTTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:575917/228‑1 (MQ=255)
                   ccAACCCCTCCGCGCAGATCTCGCCAGGTTTCGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAGGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaa                                                                                                                                                                                                  <  2:668014/228‑1 (MQ=255)
                      aCGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGCCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaat                                                                                                                                                                                               >  2:298204/1‑228 (MQ=255)
                        cccaccggccgggcggcgcGAGGTGTTTTTACCAAAGGCTACCCACGACCCAAGATGCGCGTCCGGGAATTTCCCTTCAAAAGTTGCCAGGTCGTCGGCAAATGTTTTACTGGCCGCCCGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataa                                                                                                                                                                                             <  1:318076/213‑1 (MQ=255)
                                  cgggcgagcCGGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACg                                                                                                                                                                                   <  2:248930/225‑1 (MQ=255)
                                  agggggcgcCGGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAGCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACg                                                                                                                                                                                   <  1:429482/223‑1 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCCCGCGGCAATGTTCTGGCAAAATGCCACCCTGAACCCGAGCCATTTTTCCTCCTGAAAATAAAAATGcccccca                                                                                                                                                                                 >  1:387333/1‑223 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATACCCGCGCGGCAATTTTCTGGCAAAATGCCACTCTGAACCCAGACCTTGTTCCCCCTGAAAATATAAATTCCAcccca                                                                                                                                                                                >  1:591206/1‑224 (MQ=255)
                                    agCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTAACGTTGGCTTCGATCCAAATCGCCCCGCGGCAATGTTCTGGCAACAGGCCACCCTGAACCGAGACATTGTTCCCCCTGAAAATAATAATGCCCCGCaa                                                                                                                                                                                >  1:169428/1‑228 (MQ=255)
                                       acgcCAGGTGTTGGTACCAAAGGCAACCAACGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTTGCGGCAAATGTTTTACGGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATc                                                                                                                                                                              <  2:405942/227‑1 (MQ=255)
                                         gcCAGGTTTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCg                                                                                                                                                                            <  2:480227/228‑1 (MQ=255)
                                                gttgttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCa                                                                                                                                                                       <  2:582189/226‑1 (MQ=255)
                                                   gttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgtt                                                                                                                                                                  >  2:115224/1‑228 (MQ=255)
                                                   gttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAAGCCCACCCAATCCGTGGCAttgtt                                                                                                                                                                  >  2:460874/1‑228 (MQ=255)
                                                    ttACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttgtt                                                                                                                                                                  <  2:670604/227‑1 (MQ=255)
                                                       ccAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGt                                                                                                                                                                           >  1:39594/1‑215 (MQ=255)
                                                        cAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTg                                                                                                                                                                          <  2:39594/215‑1 (MQ=255)
                                                               aaCCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAtttt                                                                                                                                                      >  2:274622/1‑228 (MQ=255)
                                                                accaccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGGGGCATGTTTGAGCGCTCtttt                                                                                                                                                      >  2:136695/1‑227 (MQ=255)
                                                                       cACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAACTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGCAATGGTTGAGCGCTATTTTACGTGgtt                                                                                                                                              >  2:334352/1‑228 (MQ=255)
                                                                                      ccGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACCCTGAACCTGAGACATTGTTCCTCCTGAAAATATTAATGCCACGCAATCCGTGCCATTTTTGAGCGCTATTTTACGTGTTTTTTAGCGCAGGCGGt                                                                                                                               >  2:691006/1‑228 (MQ=255)
                                                                                         ggAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGTCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAACTAATAATGCCACGCAATCCGTGGAATGGTTTCGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGtgtg                                                                                                                            >  2:330586/1‑228 (MQ=255)
                                                                                                               gtgtttaaggtacagtGTTTTCCTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCAAGTTGGCTGCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccg                                                                                                       <  2:680952/212‑1 (MQ=255)
                                                                                                               cAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccgc                                                                                                      <  2:387333/228‑1 (MQ=255)
                                                                                                                 gTTTGGCGGCAAATGTTTTACGGGCCGTCCGCAGGGCGTCAACTTCCCCTTTCCCGTTGGCCTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCa                                                                                                    <  2:680423/228‑1 (MQ=255)
                                                                                                                  tttGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGCCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGGTTTTTCAACGGCAAATTaaaccccccaa                                                                                                   >  1:45652/1‑219 (MQ=255)
                                                                                                                   ttGTCCGTAAAGTTTTTTCCGGCCGCACGCGGGGGGTGAACTTCCCCGTTCCCGTTGGCTCCGGACCAAATCGCCGCGCGGCAATGTGCTGGCAAAATGCCCCTCTGATCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                                                                                                  <  1:361278/228‑1 (MQ=255)
                                                                                                                   ttGGCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                                                                                                  <  2:266132/228‑1 (MQ=255)
                                                                                                                       cGGCAAATTTTTTACTGGCCGCACGCGGGGCGCCAACTTCCCCTTGCACTGTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATa                                                                                                <  1:115224/226‑1 (MQ=255)
                                                                                                                        ggCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTCTTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCGGATAAt                                                                                              >  2:427200/1‑227 (MQ=255)
                                                                                                                            aaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGCGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGGGTTTTTTCGGGGCGGCGGTTTGTTCGGAGGCAAATTAACGCCCCCAGATAAttttcc                                                                                         >  2:502971/1‑226 (MQ=255)
                                                                                                                               gTTTTGCTGGGCGCACGCAGGGCGTCAAATGCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGa                                                                                      <  2:563848/228‑1 (MQ=255)
                                                                                                                                        gCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTTGCTTCGATCCAAATCGCCGCGCGGCACTGTTCTTGCCAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCACTCCGTGGCATTGTTGAGCGCTCTTTTCAGTGTTTTTTACCGGCGGGTGTTTGTTTCGGCGCAATATTACCCCGCACGTTCATTTTACTGACTTTcgctc                                                                             >  1:330735/1‑224 (MQ=255)
                                                                                                                                        gCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCCCCAGATAATTTTACTGACTTTCCAGt                                                                             >  2:154722/1‑228 (MQ=255)
                                                                                                                                         ccGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAACTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGCCGGCGGTTTGTTCAGCCTCAAATTAACGCCTCCCGATAAtttctttctcttttatac                                                                             >  1:680952/1‑212 (MQ=255)
                                                                                                                                         ccGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGt                                                                             >  2:228635/1‑227 (MQ=255)
                                                                                                                                          cgcacgcgGGGCGTCAAGTTCCCCGTTCACCTTGGCTTCGAACCAAATCCGCGCGGGGCAATGGTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGGGGCATTGTTGAGCGCTATTTGACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGACGCCAGAGAATTTTACTGACGTTCCAGttt                                                                           <  1:671523/228‑1 (MQ=255)
                                                                                                                                          cgcacgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGCGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCCTTGTTTAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGGCGCAAATTAACGCCGCCAGATAATTTTTTTGAGTTTCCCGt                                                                             >  1:634605/1‑226 (MQ=255)
                                                                                                                                                  gggCGTCCACTTCCCCTTTCACGTTGGCTTAGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGACAAATTGTTACTCCTGCAACTAATAATGCGTCGCAATCCGTGGCTTTGTTGAGGGCTATTTTACGTGTTTTTTCGCGGCGGCGGGTTTTTCATACGCAAATTAAACCCTCCCTACAATATTCCTGGCTTTCCTGtttttttaggg                                                                   >  1:592826/1‑222 (MQ=255)
                                                                                                                                                  gggCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGGTTTTCAGGGt                                                                   >  2:117647/1‑228 (MQ=255)
                                                                                                                                                   gtggTCAATTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGGGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGGCATTGTTCCTGCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTa                                                                  <  1:274622/225‑1 (MQ=255)
                                                                                                                                                   ggCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGAAAATTAACGCCCCCCAATAATTTTACTGGCTTTCCAGGTTTTAAgggaa                                                                  >  1:446195/1‑226 (MQ=255)
                                                                                                                                                         aaCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGCTTGTTCAGACGCAAATTAACGCCGCCAGATATTTTTACTGACTTTCCAGTTTTTCAGGGAATCggc                                                             >  2:602146/1‑225 (MQ=255)
                                                                                                                                                          aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:591206/228‑1 (MQ=255)
                                                                                                                                                          aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:169428/228‑1 (MQ=255)
                                                                                                                                                            ttCCCCTTGCACGTTGGCTTCGATCCAAATCGCCGCGCCGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGa                                                          <  2:676054/227‑1 (MQ=255)
                                                                                                                                                             tCCCATTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTACTACTGAAAATAATAATGCCACGCAAACCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAAACGCACAATAACGCAGCAAGATAGTTTTATTGGCTTTCCAGTTTTTCAGGGCATCGTTCGAAg                                                        >  2:87592/1‑228 (MQ=255)
                                                                                                                                                                   ttCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCAtt                                                   >  1:625079/1‑227 (MQ=255)
                                                                                                                                                                    ttccGTTGGCTGCGATCCAAATCGCCACGCGGCAATTTTCGGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTTGCATTGTTGAGCTCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTGTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:556133/225‑1 (MQ=255)
                                                                                                                                                                    tctcGTTGGCTGCTAACCAAATCGCCGCGCGGCAATGTTCTGGCAAAAGGCCACTCTGACCCTGAGACATTGTTCCTCCTGAAAAGAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:502971/225‑1 (MQ=255)
                                                                                                                                                                    tcccGTTGGCGTCGATCCAAATCGCCGCGCGGCAATGTTCTCGAAAAAGGCCACTCTGAACCTTGGACATTGTTCCTCCTGAAAATAATTATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTaa                                                 <  1:136695/225‑1 (MQ=255)
                                                                                                                                                                         ttGGATTCGATCCAAATCGCCGTGCTGCAAGTTTCTGTCAAAATGCCACCCTGAACCTGAGCCATTGTTCCTCCTGAAAATAATAATTCCACGGAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCTGTTTGTCCAGACGCAAATTAACGCCGCCAGATAATTTTACTTACTTTCCAGTTTGTCAGGGTATCGTCCGAAGGCATTAACCCtt                                            <  1:688260/228‑1 (MQ=255)
                                                                                                                                                                                  aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTAGGTGTTTTTTCGCGGCGGCGGGTTGTTTTGACGCAAATTAACGCCGCCAGATCATTTTACTGGCTTTCCAGGTTTTCAGGTTACCTCCCTAAGGGATTAACCCtaccatccca                                    >  1:677228/1‑218 (MQ=255)
                                                                                                                                                                                  aTCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCCAAGGCATTAACCCTTCAGGCCCa                                    >  2:621938/1‑227 (MQ=255)
                                                                                                                                                                                                    cAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACt                 >  2:124204/1‑228 (MQ=255)
                                                                                                                                                                                                     aaTGTTCGGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGaa                   <  1:124204/225‑1 (MQ=255)
                                                                                                                                                                                                       ttttCTTGCAAAATTCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:253505/226‑1 (MQ=255)
                                                                                                                                                                                                       tGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:625079/228‑1 (MQ=255)
                                                                                                                                                                                                       tGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  2:446195/228‑1 (MQ=255)
                                                                                                                                                                                                       tGTTATGGAAAAATGCCCATCTGAACTTTGGACTTTTGTCCTCCTGAAAATAATTATGCCTCGCAATCCGTGTCATTGGTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGGCGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGt              <  1:691006/228‑1 (MQ=255)
                                                                                                                                                                                                           ttGGGAAAATGCCGCCCGGAACCTGGAGATTTGTTCCTCCTAAAAATCATAATGCCGCGCAATCCGTGGCATTATTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGTGGTTTGTTCAGACGCAAATTAACGCCGCCAAATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCt          <  1:125814/227‑1 (MQ=255)
                                                                                                                                                                                                              gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACTTGTTTTTTCGCGCCGGAGTTTTGTTCAGACGCAAATTAACGCCGCCAGTTAATTTTACTGACTTTCCAGTTTTTCAGTGAATCGTCCGCAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGCACTg                >  2:410426/1‑219 (MQ=255)
                                                                                                                                                                                                              gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCATGGCATCGTCCGAAGGCTTTAACCCTTCCGGACCACGCCATTCCCCCCTCAACgggatgc           >  2:125814/1‑220 (MQ=255)
                                                                                                                                                                                                                   aTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGtt  >  2:166268/1‑228 (MQ=255)
                                                                                                                                                                                                                   aTGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGtt  >  2:227780/1‑228 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                     
ATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTT  >  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: