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

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