Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,538 Δ1 bp 5.8% 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,5380T. 5.8% 209.2 / 3.3 69coding (8/900 nt)yfeXporphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base . (2/2);  total (33/36)
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.51e-01

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

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: ‑