Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,611,4550GA23.2% 144.2 / 16.7 69C93C (TGC→TGTyhhSputative transporter YhhS
Reads supporting (aligned to +/- strand):  ref base G (15/38);  new base A (16/0);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.57e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGTTT  >  NC_000913/3611218‑3611704
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
cacGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTTATTCCGGCGGCCAATACCCCAGACCGCCCCAAAAACAGCCGCATAAAc                                                                                                                                                                                                                                                 >  2:130345/1‑251 (MQ=255)
   gacgCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGATAATCTTGCCGCCAGATACCCCACACCGCTCCAAAAGAAGCCGCATAAACCg                                                                                                                                                                                                                                               >  1:256086/1‑250 (MQ=255)
       ccTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGaaaac                                                                                                                                                                                                                                          <  2:106572/251‑4 (MQ=255)
          aaCGGCGCACCCATCGCCATCGCCCTGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCCGCCGTCCAATACCCCCGACCCCCCAAAAAACAGCCGCCTAAACCGAAgacgac                                                                                                                                                                                                                                       >  1:402460/1‑251 (MQ=255)
            cGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACaa                                                                                                                                                                                                                                     <  2:63906/251‑1 (MQ=255)
            cGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACaa                                                                                                                                                                                                                                     <  2:410908/251‑1 (MQ=255)
             ggCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAAACCCGCCCGCAGATAACCCACACCACCCCAAAAACAGACGCACAATCCGAAGAcaacaac                                                                                                                                                                                                                                    >  2:357485/1‑250 (MQ=255)
              gcgcACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTTATTCTGCCGTCAGATACCCCAGAACGCTCAAAAAGCAGCCGCATAAACCGAAGAAGACCATc                                                                                                                                                                                                                                   >  1:365084/1‑251 (MQ=255)
                 cgcccATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCACATAAACCGAAGACGACAATCttt                                                                                                                                                                                                                                <  2:484060/249‑1 (MQ=255)
                 cccccATCGCCATCGCCCCTTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCATCGTGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAACCTCTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATGttt                                                                                                                                                                                                                                <  2:215369/249‑1 (MQ=255)
                  aCCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAAACCTGCCGGCAGATACCCCAGACCGCTCAAAAAGCAGACGCATAAAACGAAGACGAAAATCtttt                                                                                                                                                                                                                               >  1:353263/1‑251 (MQ=255)
                   cccATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCttttt                                                                                                                                                                                                                              <  1:45779/251‑1 (MQ=255)
                         gctatcgccCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGTTATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAAAAACCCCAGAACCCCCAAAAAACAACCCCATAAAACGAAAAAGACAATCTTTTTTGGTcc                                                                                                                                                                                                                       >  2:392555/4‑251 (MQ=255)
                                  ccGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATc                                                                                                                                                                                                              <  1:454345/251‑1 (MQ=255)
                                              aTGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTcc                                                                                                                                                                                                  <  1:535475/251‑1 (MQ=255)
                                                      ccACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATACCGCCGTCAGAAAACCCCGAACGCTCAAAAAACAGCCCCATAAACCCGAGAACACAATCTTTTTTGGGCCCCGCCAACCAGCGTAAcggccggcctg                                                                                                                                                                                            >  2:564451/1‑246 (MQ=255)
                                                         cGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGggc                                                                                                                                                                                       <  1:130345/251‑1 (MQ=255)
                                                         cGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGggc                                                                                                                                                                                       <  1:554865/251‑1 (MQ=255)
                                                               cACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTc                                                                                                                                                                                  <  1:118739/250‑1 (MQ=255)
                                                                  ccgcccgATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCa                                                                                                                                                                              <  1:561884/251‑1 (MQ=255)
                                                                     cccgATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAgg                                                                                                                                                                           <  1:587187/251‑1 (MQ=255)
                                                                        gATAGGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAAGCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAAGCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTgg                                                                                                                                                                        <  1:261721/251‑1 (MQ=255)
                                                                         atatGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGAAACCCCAGAACGCTCAAAAAAAAGCCGCATAAAACGAAGACGGCAAACTTTTTTGGGCCCACCGAAACGGCGTCACGTCCCGCAAGAGGGCGGCTCCGCAAgggggg                                                                                                                                                                       >  2:448864/1‑247 (MQ=255)
                                                                            tGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTTCCGTCAGATACCCCCGAACGCCCCAAAAACAACCGCATAAAACGAAGACGACAATCTTTTTTGGGCCCCGCCAATACGCGTAAAGTCCCGCCAGAGGGCGGGCCAGCAAGGGGGCGa                                                                                                                                                                     >  2:204599/1‑250 (MQ=255)
                                                                              cAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTGCCGCCAGATACCCCAGACCGCTCCAAAAACAGCCGCATAAAACGAAGACGACAATCTTTTTTGGGCCCAGCGAATCGGCGGAACGTCCGGCATGAGGGCGGCTCAGCAAGGGGGCGAAat                                                                                                                                                                  >  1:226756/1‑251 (MQ=255)
                                                                                  gAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAAAAACCCCCGACCGCTCAAAAAACAGCCGCATAAACCGAAGACGACAAAATTTTTGGGGCCCAACGAATAGGCGGAAAGGCCGGCATGAGGGCGGACCAGCAAAGTGGGGAAATATTg                                                                                                                                                              >  1:373049/1‑251 (MQ=255)
                                                                                  gAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTAAATCACGACCTCACATACCCAACAACGCTCAAAAAAAAGCCGCATAAACCCAAGACAACAATATTTTTGGGGCCCAGAGAATCGGCGTAACGGCCGGCCTGAGgggggcgggacgccagaaattaaggg                                                                                                                                                                >  2:206037/1‑229 (MQ=255)
                                                                                   aGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGGCAGATACCCCCGACCGCTCAAAAAGAAGACGCATAAAACCAAGACGAAAATCTTTTTTGGGTCCAGCGAAACGGCGGAACGTCCGGCATGAGGGCGGCCCAGCAAGGGGGGGAAAAAttt                                                                                                                                                              >  1:170086/1‑249 (MQ=255)
                                                                                   aGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTATGCCAAATCCCAAGGATGACGCACCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCTTACAGCCCAAAAACCCAAACAGCTCAAAAAACAACCACCTAAAACCGAGACGACAATCTTTTTGGGGTCCATCGAAACGGCGTGACGCCCGGCATGGGGGCGGCGCCGCAATGTGGCggaagaatg                                                                                                                                                              >  1:568368/1‑242 (MQ=255)
                                                                                     ccAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAg                                                                                                                                                          <  1:392555/251‑1 (MQ=255)
                                                                                          cACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCCGACCGCCCAAAAAACAGCCCCATAAACCGAAGACGACAAACTTTTTTGGTCCCCGCGAATCGGCGGAACGTCCGGCCTGCGGGCGGCTCAGCAAGGGGGCGAAATATTTCAGGCTg                                                                                                                                                      >  2:224086/1‑250 (MQ=255)
                                                                                               cAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAAcc                                                                                                                                               <  1:159845/251‑1 (MQ=255)
                                                                                                     cccATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTTCCGTCAGATACCCCCAACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTTGGTCCCAGCGAATCGGCGTa                                                                                                                                                                                                        >  2:496297/1‑190 (MQ=255)
                                                                                                     cccATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTa                                                                                                                                                                                                        <  1:496297/190‑1 (MQ=255)
                                                                                                     cccATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAAGCCCAAGGATGACTCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGGCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCt                                                                                                                                         <  2:412730/251‑1 (MQ=255)
                                                                                                                 aTCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGGGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGgcg                                                                                                                             >  1:414080/1‑251 (MQ=255)
                                                                                                                   cccGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAACAACCGCATAAACCGAAGACGACAATCTTTTTTGGGCCCAGCGAAACGGCGTAAAGGCCGGCATGAGGGCGGCGCAGCAAGGGGGGGGAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCt                                                                                                                           >  1:384380/1‑251 (MQ=255)
                                                                                                                         ccGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGAAACCCCCGACCGCTCAAAAAACAGCCGCATAAAACGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGGCCGGCATGAGGGGGGCTCAGCCAGGGGGGGAAATATTGACAGCTGATAAACCAACCCGCCCCGAAAGCGCTAAAGcc                                                                                                                     >  1:255591/1‑251 (MQ=255)
                                                                                                                           ggCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAg                                                                                                                                                                                                                    <  1:550907/156‑1 (MQ=255)
                                                                                                                           ggCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGGTAATCCTGCCGTCCAATACCCCCGAACGCCCAAAAAGCAGCCGCCTAAACCCAAGACCACAATCTTTTTTGGTCCCAg                                                                                                                                                                                                                    >  2:550907/1‑156 (MQ=255)
                                                                                                                             cAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCGTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTgg                                                                                                                                                                        <  2:280076/198‑1 (MQ=255)
                                                                                                                             aaaaaCTTTGCCCAATCCCAAGTATGACGCGCCCCAGGCAAAGGAATAACAGGCGGATGACAGGCCGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTATTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCATTCCTGCCCAGAAGGCGCTAAAGCCcatc                                                                                                                 <  2:568368/250‑1 (MQ=255)
                                                                                                                              aaaaCTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTCTGGGTCCCAGCg                                                                                                                                                                                                                  <  1:425148/155‑1 (MQ=255)
                                                                                                                               aaaCTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCcatcac                                                                                                               <  1:55354/251‑1 (MQ=255)
                                                                                                                                      gcgcAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCCGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCCGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCGTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGATGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAATGCCCATCACATCATg                                                                                                         <  2:220157/247‑1 (MQ=255)
                                                                                                                                      gCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATTCTGCCGTCAGATACCCCAGACCCCTAAAAAAAAAGCCGCATAAAACGAAAACGACAAACTTTTTGGGTCCCAGCGAATCGGGGTCACGTCCGGCATGAGGGCGGATCAGCAAGGGGGCGAAATAATGCAGGCTGATAAACAAACCTGCCCACAAGGCGCTAAAACCCATCACAACATgg                                                                                                        >  1:384592/1‑251 (MQ=255)
                                                                                                                                      gCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATgg                                                                                                        <  2:225282/251‑1 (MQ=255)
                                                                                                                                                   gATGGCGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTa                                                                                           <  2:516909/251‑1 (MQ=255)
                                                                                                                                                   gATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGGTAATCCTGCCGTCAGAAACCCCCGGACGACCAAAAAGAAACCGCATAAACCGAAGACGACAAACTTTTTTGGGCCCCGCGAATCGGCGTAAAAGCCGGCATGAGGGGGGCTCAGCAAGGGGGCGAAAAATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCCAAAACCCCTCACCTCCTTGACATAGCCCGGt                                                                                            >  1:516909/1‑250 (MQ=255)
                                                                                                                                                                 ggCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGGCAGATACCCCAAACCGCTCAAAAAACAGGCGCATAAAACCAAGACGACAATCTTTTTGGGGCCCAACGAAACGGCGTAAAGTCCCGCAAGGAGGCGGCGCAGCCAGGGGGCGGAATATTGCAGGCTGATAACCAAACCTTCCCAGAAGGCGCTCAAACCCAACACAACAAGGAAATAGCCCCGCAATACAGCGAGCGGc                                                                             >  2:382202/1‑251 (MQ=255)
                                                                                                                                                                   cAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCaa                                                                           <  1:224086/251‑1 (MQ=255)
                                                                                                                                                                     aaGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAAcc                                                                         <  1:517007/251‑1 (MQ=255)
                                                                                                                                                                       gTAAAAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCAGGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGGGCGGCAACCCg                                                                       <  2:210233/251‑1 (MQ=255)
                                                                                                                                                                          ataaCAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATg                                                                    <  1:579672/251‑1 (MQ=255)
                                                                                                                                                                            aaCAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGCAGACGACAAACTTTTTGGGTCCCAGCGAATCGGCGGAACGTCCGGCATGAGGGCGGCCCAGCCAGGAGGCGAAATATTGCAGGCTGGTAACCAAACCTGCCCAGAAGGCGCTAAAGCCCATCACATCAAGGACATAGCCCGGGAATACAGCGAGCGGCAACCCGagggg                                                                  >  1:198767/1‑247 (MQ=255)
                                                                                                                                                                             aCAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTg                                                                 <  1:204599/251‑1 (MQ=255)
                                                                                                                                                                                     atgaCAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCt                                                         <  2:255591/251‑1 (MQ=255)
                                                                                                                                                                                        aCAGGCAGACTGGCGGTTAATCCTGCCGGCAGATACCCCAGACCGCTCAAAAAGCAACCGCATAAACCGAAGACGCCAATCTTTTTTGGTCCCAGCGAAACGGCGGAACGTCCGGCGTGAGGGCGGGTCAGCAAGGGGGCGAAATATTGCAGGCTGGTAACCCATCCTGCCCAGAAGGCGCTAAAGCCCATCACAACATGGACAAAACCCGGGAATAACGCGAGCGGCAACCCGATGGGGGGGTAGGTGGc                                                      >  1:514417/1‑251 (MQ=255)
                                                                                                                                                                                          aGGCAGACTTGCGGTTAATCCTGCCGGCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAGCCCGATGGTGAGGTAGCTGGCGa                                                    <  2:9313/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAAc                                            >  1:93334/1‑250 (MQ=255)
                                                                                                                                                                                                       gTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGAc                                       >  2:361776/1‑251 (MQ=255)
                                                                                                                                                                                                           aaTCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACtat                                    >  1:432529/1‑251 (MQ=255)
                                                                                                                                                                                                           aaTCCGGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACta                                     <  2:365084/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGAc                              <  2:256086/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAgg              <  1:219565/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTCAAAAAGCAACCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGgcgc           >  1:422977/1‑250 (MQ=255)
                                                                                                                                                                                                                                         aaaaaGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATcc      >  1:387401/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aaaaaGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATGCAGGCGCAATcc      <  2:226756/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAAGACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCg     <  2:468631/251‑1 (MQ=255)
                                                                                                                                                                                                                                             agcagcCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGttt  >  2:308505/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                             
CACGACGCCTAACGGCGCACCCATCGCCATCGCCCCGTAAGTGACAATGCCGTTCCACGAAATCACCCGCCCGATATGCAGCGAGCCAACCACGCCAACGCCCCATAGGGTCGATCCCGTTCCGGCAAAACTTTGCCCAATCCCAAGGATGACGCGCCCCAGGCAAAGTAATAACAGGCTGATGACAGGCAGACTGGCGGTTAATCCTGCCGTCAGATACCCCAGACCGCTCAAAAAGCAGCCGCATAAACCGAAGACGACAATCTTTTTGGGTCCCAGCGAATCGGCGTAACGTCCGGCATGAGGGCGGCTCAGCAAGGTGGCGAAATATTGCAGGCTGATAACCAATCCTGCCCAGAAGGCGCTAAAGCCCATCACATCATGGACATAGCCCGGTAATACAGCGAGCGGCAACCCGATGGTGAGGTAGCTGGCGAAGTTAAACATGACTATAGAGACAATGCGCAAATTCAGGCGCAATCCGTTT  >  NC_000913/3611218‑3611704

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