Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,7150GA34.0% 78.5 / 16.6 50intergenic (‑96/+45)uspF/ompNuniversal stress protein F/outer membrane porin N
Reads supporting (aligned to +/- strand):  ref base G (28/5);  new base A (0/17);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 2.67e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACC  >  NC_000913/1435466‑1435963
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         
aaCCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTTCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCTATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAGCAAAggagga                                                                                                                                                                                                                                                          <  2:530049/251‑1 (MQ=255)
 aCCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTGCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGgaa                                                                                                                                                                                                                                                         <  2:15484/251‑1 (MQ=255)
         agGCAGAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGGGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGATACTTTTTGTGGAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGc                                                                                                                                                                                                                                                 <  1:733554/251‑1 (MQ=255)
             aTAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGa                                                                                                                                                                                                                                              >  1:571207/1‑250 (MQ=255)
                  agGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                                                                                                                                        >  1:196438/1‑251 (MQ=255)
                   gggCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGcc                                                                                                                                                                                                                                       >  2:792347/1‑251 (MQ=255)
                       aGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                                                                                                                  <  2:194876/251‑1 (MQ=255)
                        gTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTGATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGcc                                                                                                                                                                                                                                  >  1:107346/1‑250 (MQ=255)
                                 ttACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCgg                                                                                                                                                                                                                        >  2:498251/1‑251 (MQ=255)
                                  tACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCggg                                                                                                                                                                                                                       >  2:771786/1‑251 (MQ=255)
                                           gAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                                                                                                              >  2:72415/1‑251 (MQ=255)
                                            aaaTGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAtt                                                                                                                                                                                                             >  2:178646/1‑251 (MQ=255)
                                                 aaCCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                                                                                                         >  2:659691/1‑250 (MQ=255)
                                                        gCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCTTTGATTTAATTCGGAATCGAAAATATCGATAGGGAAAAGAATCGTTCTGTTCATAAACCCTCCTGTTTTAGTATCCGAATAAAGTGAAACGCCAGATGACATTTTTTGTGTAAGGACGGAGTTCACATTTTAAATTAAGATAAAAGGAGGAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                                                                                                                   <  1:74735/249‑1 (MQ=255)
                                                        ccatcatCAATTTTTGCCTCTCCCTAAACGGGCATAACCGCGGGTGAATTTAATTCTAAAGCGAAAATATGAATAGGGACAAGAAGCGTTCTGTTAAAAAAACCTCCTGTTTTAGTACCCGAATAAAGTGTAACGCCAGAGGACATTTTTTGTGAAAGGAGGGAGTTCACATTTTTAATTTAGATAAAAGGAGGAAAAAAAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                                                                                                                   <  1:684961/248‑1 (MQ=255)
                                                                aaTCTTGGCCTCTCCCTCACGGGGGCAAATCACGCGTTGAGTTAATTCTAAATCGAAAATATGGAAGGAGAAAAGAATCTTTCTGTAAAAAAACCCTCTTTTTTTAGTACCCGAAAAAAGTATAACGCCAAATGACTCTTTTTGTGTAATGACGAAGTCCATTTTTTTAATTAAGATAAAGGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                                                                                                         <  2:490378/251‑1 (MQ=255)
                                                                     ttcccTCTTCCTAACCTTGGCAAATCACGTGTTGATTTAATTTTGAACCGAAAATATCGATAGGGAAAAAAACCTTTCTTTTAAAAAACCCGCTTTTTTTAGTACCCGAAAAAATTGTAACGCAAGAGAACATTTTTTGTGAAATGACGGAGTTCACTTTTTAATTTTAAATAAAGGGAGAAAGAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                                                                                                                                                                     <  1:784893/247‑1 (MQ=255)
                                                                     ttGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGCCAAAAATCGTTTTGTTAAAAAACCCTCTTGTTTAAGAATCCCAATAAAGAGTACCGCCAGAGGACACTTTTTGGGTAATGAGGGAGTTACCATTTTAATTTTAGATCAAAGGAGGAAAAATAAGCAGAAAAAGCCCGCAATAAAAGCGGGCAGGAGGATTTAGAACTGATAAACAAGCCCTaa                                                                                                                                                                                      >  1:212024/1‑249 (MQ=255)
                                                                       cccttttcccaaaCGTGGCAAATCACGCGTTGAGTTAATGCTGAACGGGAAATATCGATAGGAACAAGATTCGTTTTGTTCATAAACCCTCTTTTTTAAGTATCCGAAAAAAGTGTAACGCCAGATGACATTTTTTGTGTAATGACGGAGTTCACATTTTTAATTAAGATCAAAGGAGAAAGAAGAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                                                                                                                   <  1:638540/239‑1 (MQ=255)
                                                                           cttctaaaccGTGGCTAATCCCGCGTTGAGTTAATTCTAAACCGAAAATTTCGATAGGGAAAAAAATCGTTTTGTTCAAAAAACCTCCTTTTTTAGTATCCGCATAAAGTGAAACGCAAGATGACTTTTTTTGTGAAATGACGGAGTTCACATTTTTATTTAAGATAAAGGGAGAAAAAATAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                                                                                                                                <  1:178646/240‑1 (MQ=255)
                                                                                        gCTAATCACGCGTTGATTTAATTCTAAATAGAAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAACCCTCTTGTTTTAGTACCCGCATAAAGTGTAACGCCAGATGACATTTTTTGTGTAATGACGGAGTTCAAATTTTTATTTAAAAAAAAAGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTaga                                                                                                                                                                 <  1:359201/251‑1 (MQ=255)
                                                                                            aatccggcgTTGAGTTAATTCAGAAGCGAAAATATCAATAGGAAAAAGAATCTTTCTTTTAATAAACCCCCCTTTTTTAGTATCCGAAAAAAGTGTACCCCCAATGGACTTTTTTTGTGTAATGACGGATTTCAATTTTTTAATTTAGAAAAAAGGAGAAAAAATAAGAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                                                                                                                                                                <  1:807241/243‑1 (MQ=255)
                                                                                            aaTCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCCTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                                                                                                                  >  1:686712/1‑247 (MQ=255)
                                                                                                      tGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTACTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCAGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACACTATCATCGGTAGAGATGCCATTGGc                                                                                                                                                     >  2:649462/1‑250 (MQ=255)
                                                                                                           tAATTATGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGAAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGTCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGAAGAGATGCCATTGGCAGCGt                                                                                                                                                >  1:471359/1‑250 (MQ=255)
                                                                                                              ttCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                                                                                                             >  2:875104/1‑250 (MQ=255)
                                                                                                                  gAATCGAAAATATCAAGAGGGAAAAGAATCGTTTTGTTCAAAAACCCTCTTGTTTAAGTTCCCGCAAAAAGTGAAACGCCGGAGGACTTTTTTTGTTAAATGCCGGGTTTCACTTTTTTAATTTAAACAAAAGGAGAAAAAAAAAGAAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCt                                                                                                                                         <  1:771786/250‑1 (MQ=255)
                                                                                                                            tatCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                                                                                                              >  2:380867/1‑251 (MQ=255)
                                                                                                                                   tagagACAAGAATCGTTCTGTTCATAAAACCTCATGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGATCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATAGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                                         >  2:60333/5‑250 (MQ=255)
                                                                                                                                     ggACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACGGATAAACCAGACCTAAAGCGACAATATCATCGATAGAGATGCCATTGGCAGCGTAGAAGATGTCATCTTCATCCAACAGGt                                                                                                                     >  2:42443/1‑251 (MQ=255)
                                                                                                                                                    tttttaataaaccttctttttttAGTAACCACAAAAATTGAAACGCCAGTTAACATTTGTTGTGGAAGGACGGGTTTCATTTTTTAAATTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc                                                                                                        <  2:571207/231‑1 (MQ=255)
                                                                                                                                                              acccccctttttaagtacccgaaaaaaGTGTAACGCCAGATGACTTTTTTTGTTAAATGAGGGATTTCACATTTTTAATTTAAAAAAAGGGAGAAAAAAAAAGCAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg                                                                                             <  1:792347/227‑1 (MQ=255)
                                                                                                                                                               ccctcctGTTTTATTTCCCGCAAAAGTTGAACCGCCAAATGACCTTTTTTGTGTAATGACGGAGTCCAATTTTTTATTTAAAATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGa                                                                                            <  1:807952/250‑1 (MQ=255)
                                                                                                                                                                   cctGTTTTAGTATCCGCATAAAGTGTAACGCCAGCTGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATg                                                                                        >  2:224287/1‑251 (MQ=255)
                                                                                                                                                                        tttattacccgtaaaaaGTGTAACGCCAGATTATATTTTTAGTAAAAGGAGGGAGTTCATTTTTTAAATTAAAAGCAAAGGAGAAAGAATAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGtttt                                                                                    <  2:90161/236‑1 (MQ=255)
                                                                                                                                                                            gTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTg                                                                               >  2:424777/1‑251 (MQ=255)
                                                                                                                                                                              acccgaaTAAAGTGAAACGCAAGATGAACTTTTTTGTGAAAGGACGGATTTAACATTTATAATTAAAAAAAAAGAAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGa                                                                              <  2:196438/244‑1 (MQ=255)
                                                                                                                                                                                cagaataaatgtaaaCGCAAGTTGACTCTTTTTGGATAATGTCGGATTTAACATTTTAAATTAAGAAAAAAGGAGGAAAAAAAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGaaa                                                                            <  2:881510/237‑1 (MQ=255)
                                                                                                                                                                                  gCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAAt                                                                           >  1:277863/1‑249 (MQ=255)
                                                                                                                                                                                    aTAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGt                                                                        >  1:670174/1‑250 (MQ=255)
                                                                                                                                                                                                  ccAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAat                                                           >  1:51151/1‑250 (MQ=255)
                                                                                                                                                                                                    aGATGACACTTTTTGTGTAAAGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCTCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCCTCCAACAGGTTGATTTTATAGTCCACGTAGGTGGACATGTTTTTATTGAAATAGCAAGTCGCGCCAATATc                                                        >  2:200341/1‑251 (MQ=255)
                                                                                                                                                                                                     agaaaactttttttGTGAAATGAGGGAGTTCAAATTTTTATTTTAGATAAAAGAAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATc                                                        <  2:670174/242‑1 (MQ=255)
                                                                                                                                                                                                       tGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc                                                     >  1:228397/1‑251 (MQ=255)
                                                                                                                                                                                                        gACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                                         >  1:402033/1‑182 (MQ=255)
                                                                                                                                                                                                         gactttttttGTTTAAGGAGGTAGTTCACTTTTTTAATTTAGAAAAAAGGAGAAAAAATAAGAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc                                                     <  1:498251/245‑1 (MQ=255)
                                                                                                                                                                                                                     tgtAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGGGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGAAGGTGGACATGTTTTTATTGAATTAGTAAGTCGCGCAAATATCGGCGTATTTAAcc                                           >  2:552321/1‑249 (MQ=255)
                                                                                                                                                                                                                      gtAATGACGGATTTCAAATTTTTAATATAGAACAAAGGAGAAAAAATAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGa                                        <  1:846198/251‑1 (MQ=255)
                                                                                                                                                                                                                         aggagggagtcaacttttttAATTTAGAAAAAAGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGTTCt                                     <  1:380867/236‑1 (MQ=255)
                                                                                                                                                                                                                            aCGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCttt                                   >  1:461244/1‑250 (MQ=255)
                                                                                                                                                                                                                                            tAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCg                  >  2:583473/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     aacaaaaGGAGAAAGAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTc           <  1:875104/247‑1 (MQ=255)
                                                                                                                                                                                                                                                              aggaAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCaccacc  >  1:475759/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         
AACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACC  >  NC_000913/1435466‑1435963

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