Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,6980CA24.6% 118.1 / 12.5 69I120I (ATC→ATAgspEType II secretion system protein GspE
Reads supporting (aligned to +/- strand):  ref base C (36/16);  new base A (1/16);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.49e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAGATGTCCGCGTATCCACACTGCCGTCCATCTATG  >  NC_000913/3458454‑3458948
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
atatTCCACGTAAGGTTGGCCAGAGCCAGACCATGAACAGCTTGTCACCGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGGGTCGCAACAGATTGCGCAGGATATCGACCAATCCGTCGTTCTTCTCTCGCTTTCGGCAGAGATGACGGCAAATGAAGAACTCCTGAATGAAGTTTCAGCGGCACCGGTTATCCGATTGATCAATGCGATTTTAAGTGAGGCCATCAAAGa                                                                                                                                                                                                                                                       >  2:506430/1‑250 (MQ=255)
 tatTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGCTTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGCTTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGCTTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGCGTGGGGCCATCAACGaa                                                                                                                                                                                                                                                      >  2:399733/1‑250 (MQ=255)
  atTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAAc                                                                                                                                                                                                                                                    >  2:460913/1‑251 (MQ=255)
     cAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCg                                                                                                                                                                                                                                                  <  1:204941/250‑1 (MQ=255)
          tCAGGTTGGCCAGTGCCAGAACCTGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGCTTGCGCGGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGCTTTTGAGTGAGGCAATCAAAGAAACCGCCTCTg                                                                                                                                                                                                                                            >  2:513839/1‑251 (MQ=255)
            gggtgggCCAGTGCCAGACCATGCCCAGTTTGTCAAAGGTTGATTTTGAAGTGCGGCGGGAAGCGGTATTCCATCAGAATACGGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCTTGAAGGAAGATTCAGCGGCACCGGTTACCCGCTTGATCAATGCGATTTTGAGTGAGGCCAGCAAAGAAACCGCCTCTGat                                                                                                                                                                                                                                          <  2:296237/246‑1 (MQ=255)
                     aGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCGGGAAGCGGTATTCCATCAGAATCCCGGTGATTCGCAACAGATGGCGCAGGATATCGATCAATCCGTCGATCTTCCCTCGTTTTCGGAAGAGATGCCCGAAAATGAAGATCTCCTGAAGGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGTTTTGGAGTGAGGCCATCAAAGAACCCGCCTCTGATATCCACAtt                                                                                                                                                                                                                                 <  1:300315/251‑1 (MQ=255)
                      gTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGaggttgacgtaattggtggatcc                                                                                                                                                                                                                                 >  2:202401/1‑228 (MQ=255)
                      gTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGATCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTg                                                                                                                                                                                                                                >  1:292307/1‑251 (MQ=255)
                        tccaaaCCATGACCAGGTGTTCACGGGCTAATTTTGAAGCCCGTCGGGAGTCGGTATTCCATCAGAATACCGGTGCTTCCCAAAAGATTGCGCAGATTTGCGATCAACCCGTGGTTCTTCTCTCGCTGCCGGAAGAGATCCCCGCAAATAAAGATCTCCTGAATGAGTATTCCGCGTCACGGTTTATCCCTTTAATCATGGCGATTTGGAGTGGGGCCATAAAAGAAACCGCCTCTGATATCCACATTGa                                                                                                                                                                                                                               <  2:95781/245‑1 (MQ=255)
                          cAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATCCCGGTGAGTCGCAACAGATTGCGCAGGATAGCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAAGGAATATCTCCTGAATGAAGATTCAGCGGCCCGGGTTATCCTCTTGATCAATGCGATTTTGAGGGTGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAAc                                                                                                                                                                                                                            <  1:579290/251‑1 (MQ=255)
                           aGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAAcc                                                                                                                                                                                                                           >  2:216442/1‑251 (MQ=255)
                            gACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCctc                                                                                                                                                                                                                                              <  1:156653/231‑1 (MQ=255)
                            gACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCctc                                                                                                                                                                                                                                              >  2:156653/1‑231 (MQ=255)
                                     cAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGaaaaa                                                                                                                                                                                                                  >  2:68567/1‑250 (MQ=255)
                                       gCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCCGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGCAGAGATGCCCGCAAATGAAGATCTCCTGACTGCAGATTCTGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAGGAAACCGCCTCTGTTATCACCATTGAAACCTATGaaacaaa                                                                                                                                                                                                                >  2:157875/1‑246 (MQ=255)
                                               caGGCTAATTTGTAAGCTCGGTTGAAAGCGTTATTCCATCAGAAGCCCGGTGATTCGCAACAGATGGCGCGGGATTTCGATCAAGCCTGCGATCTTCTCTCGCTTTCGGAAGCGGTGCCCACAAATGAAGAGCTCCGAAATGAAGATTCAGCGGCCCCGGTTAGGCGCTTGTGCAATGCGATTTAAAGTAAGCCTACCAAAAAAACCTCCTCTTAGATCCACATTGAAACCTATGAAAAAACAATGTCGAt                                                                                                                                                                                                       <  1:393045/251‑1 (MQ=255)
                                                     taTTTTGAAGCTCGGCGGGAAGCGGTATTCCATCAGAATACCGGGGAGTCGCAACATAGTGCGAAGGATATCGATAAACCCGTCGATCTTCTCTCGCTTTCGGAAGAGTTGCCCGAAAATGAAGATCTCGTAAATGAGGATTCAGGGCCTCCGGTTAGCCGCTTTATCAATGCGATTTTGAGTGAGGCCATCAAAGAATCCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGtt                                                                                                                                                                                                  <  1:584824/249‑1 (MQ=255)
                                                     gATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGCGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGttt                                                                                                                                                                                                 >  1:55104/1‑251 (MQ=255)
                                                     gATTTTGAAGCTCGGCGGGAAGCGGTATTCCATCAGAATACCGGGGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAAGGCGATTTTGAGGGAGGCCATCAAAGAATCCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGt                                                                                                                                                                                                   <  1:366994/249‑1 (MQ=255)
                                                        ttggaagtccGGCTGGAAGCGGTATTCCACAAAAATACCGGTGAGTCGAAACAGATGGCGAAGGATATCGATCAATCCGTCGATCTCCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTc                                                                                                                                                                                               <  2:136388/241‑1 (MQ=255)
                                                            aaGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATc                                                                                                                                                                                          >  1:530294/1‑251 (MQ=255)
                                                                       gAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGtttt                                                                                                                                                                               >  1:18867/1‑251 (MQ=255)
                                                                        aaGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTg                                                                                                                                                                              >  1:191475/1‑251 (MQ=255)
                                                                            gttATTCCATCAGAATACCGGTGAGCCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCg                                                                                                                                                                            <  1:189533/247‑1 (MQ=255)
                                                                            ggTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCg                                                                                                                                                                            >  2:189533/1‑249 (MQ=255)
                                                                             gTATTCCATCAGCATACCGGTGAGTCGCAACAGATTGCGCAGGTTATCGATCAATCCGTCGATCTTCTCTCACTTTCGGAAGAGGTGCACGCACATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGATATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAACCCGCCTCTGAGATCCACATTGAAACCTATGAACAAACAAGGTCGATCCGTTTTCGCATctcccgctttttgcggcc                                                                                                                                                                         >  2:4804/1‑234 (MQ=255)
                                                                               aTTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGcc                                                                                                                                                                                                                                                <  1:304453/178‑1 (MQ=255)
                                                                               aTTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGcc                                                                                                                                                                                                                                                >  2:304453/1‑178 (MQ=255)
                                                                                 tccagcaaaaTACCGGTGAGTCCCAACAGATTGGGCGGGATATCACACACTCCGTGGTTCTTCTCCCGCTTTCGGAAGAGCGGCCCGCAAAGTAAATTCTCCTGAATGAAGATTCCGGGGCCCGGGTTACCCGCTTGATCATTCCGATTTTAAGGGAGGCAATAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAAt                                                                                                                                                                      <  2:550936/242‑1 (MQ=255)
                                                                                  ccatcagacaacggttgagtccaaaCAGATTGCGCAGAATATGGATCAATCCGTCTAACTTCTCTCGCTTTCGGAAGAGATCCCCGCAAATAAAGATCTCTTGAATAAAGATTCAGCGGCACTTGTTATCCTCTTGACCAATGCAATGTTGGTTGAGCCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCTTTTTCGCATCGACGGCGTTTTGCGGACAAttt                                                                                                                                                                    <  2:217663/228‑1 (MQ=255)
                                                                                            ccgtgtgactagCAACAGATGGCGAAGGATATCCATCAACCCTTCGATCTTCCCCCCCTTTGGGAAGAAATGCCCGCAAATAAAAATCTCCTAAAGGAAGATGCGGCGGCACGGGTTACCCGCTTGATCAAGGCGATTTTGAGGGAGCCCACAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCcaaa                                                                                                                                                          <  2:36423/240‑1 (MQ=255)
                                                                                            aCCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCcaa                                                                                                                                                           >  1:413792/1‑250 (MQ=255)
                                                                                                      cGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCg                                                                                                                                                                                                         >  1:279530/1‑194 (MQ=255)
                                                                                                       gCAACAGATTGTGCCGGATCTCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCAATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGAAAATTTTACAGCCAAACACAAAACTgg                                                                                                                                               >  2:75923/1‑251 (MQ=255)
                                                                                                       gCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGCGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATAGAAACCTATGAAACAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGCCAATTTTACAGCCAAACAAAAAACTgg                                                                                                                                               >  1:21907/1‑251 (MQ=255)
                                                                                                         aaCAGATTGCGCAGGATATCGATAAACCCGTCGATCTTCTCTCGTTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAGTGAAGATCCGGCGCACCCGTTTACCCGCTTGATAAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGGCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTggcg                                                                                                                                             <  1:412493/251‑1 (MQ=255)
                                                                                                                 gcgcAGGATTGCGATCAATCCGTCGATCTTCTCTCGCTGTCGGAAGAGATGCCCGAAAATGAAGATCTCCTGAATGAAGATACAGCGGCACCGGTTATCCGCTTGATAAAGGCGATTTTGAGTGAGACCATAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCt                                                                                                                                     <  2:18867/251‑1 (MQ=255)
                                                                                                                        atatCGATCAATCCGTCGATCATCTCTCGCTTTCGGAAGAGATGCCCGCAAATGGAGATCTACTGACTGCAGATTCCGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTcc                                                                                                                              >  2:181047/1‑251 (MQ=255)
                                                                                                                          atCGATCAACCCGTCGATTTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTccc                                                                                                                             <  1:145645/250‑1 (MQ=255)
                                                                                                                               aTCAACCCGTCGATCTTCTCTCGCTTTCGGAAGAGAGGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCCCCGGTTACCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAAtt                                                                                                                        <  1:460913/251‑1 (MQ=255)
                                                                                                                                tcactCCGTCGAGCTTGTCTCGCTTTCGGTTGAGCTGCCCGCAAATGATGCTCTAATGAATGAAGATTCATCGGCACCGGTTATCCGCTTGATCAATGCGATTTGGAGTGAGACCATCAAAGAAACCGAATATGATATCCACATTGAAACCAATGAAAAAACAATGTCGATCCGTTTTCGAATCGACGGCGTTTTGCGGACAATTTTACAGCCGAACAATAAACTGGGGGCACTGCTTACCTCCGGAAtt                                                                                                                        >  2:116463/5‑250 (MQ=255)
                                                                                                                                tCAATCAGACGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGACGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCCAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAAtt                                                                                                                        >  2:217119/1‑250 (MQ=255)
                                                                                                                                   acccgtctatttcctctcGCGTTGGGAAGAGAGCCCCCCAAAGTAAGATCTCCTGAATAAAGATTCAGCGGCCCCGGTTATCGGTTGTATCAATGGATTTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAgg                                                                                                                    <  2:419900/238‑1 (MQ=255)
                                                                                                                                     ccgcggttcttcccccgtttgcGGAAGAGATGCCCGAAAATGAAAATCTCATGAAGGAAGATACAGGGGCACCGGTTTCCCGCTTGATCAATGCGATTTTGAGTGAGGCCACAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTCCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTc                                                                                                                  <  2:327429/230‑1 (MQ=255)
                                                                                                                                       gTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATCAAGATATCCTGAATGAAGATTCAGCGGCACCGGTTTTCCGCTTGATCAATGCGATTTTGAGTGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAACAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATATTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCAt                                                                                                                >  1:522844/1‑251 (MQ=255)
                                                                                                                                         cGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGTCACTGCTTATCTCCCGATTTAAGGTCATgg                                                                                                              >  2:19617/1‑251 (MQ=255)
                                                                                                                                          gATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCCTGGc                                                                                                             >  2:314216/1‑251 (MQ=255)
                                                                                                                                          gATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGCCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAAAGTCGATCCGTTTTCGAATCGACGGCGTTTTGCGGACAATTTTACAGCCAAAAAAAAAACTGCCGGCACTGCTTATCTCCCGAATTAAGGTCATGGc                                                                                                             >  1:375654/1‑251 (MQ=255)
                                                                                                                                              ttctCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGt                                                                                                         >  1:43076/1‑251 (MQ=255)
                                                                                                                                                 cctcGCTTGCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAAGGAAGATTCAGCGGCCCCGGTTATCCGCTTGATAAATGCGATTTTGAGGGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCtt                                                                                                      <  1:68567/250‑1 (MQ=255)
                                                                                                                                                     gCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGat                                                                                                   >  1:423046/1‑250 (MQ=255)
                                                                                                                                                       tttCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGatatt                                                                                                >  2:560408/1‑250 (MQ=255)
                                                                                                                                                        ttCGGAAGAGATGCCCGAAAATGAAGATCTCCTGAATGAAGATTCAGCGGCCCCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGatat                                                                                                 <  1:159545/249‑1 (MQ=255)
                                                                                                                                                         gcGGAAGAGATGCCCGAAAAGGAAGATCCCCTAAATGAAGATTCGGCGGCACCGGTTACCCGCTTGATAAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCg                                                                                               <  1:489700/249‑1 (MQ=255)
                                                                                                                                                                 aaTGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCCCCGTTTATCCGTTTGATCAATGCGATTTGGAGTGAGGCAATAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAAc                                                                                      <  1:560408/250‑1 (MQ=255)
                                                                                                                                                                  aTGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGGCTTGATATCGGCGGAAAAcc                                                                                     >  2:148804/1‑250 (MQ=255)
                                                                                                                                                                  aTGCCAGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAACGAAACCGCCTCTGATATCCACATTGACACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTACGGTCATGGCTCGTTTTgcgatttccgtaaaacc                                                                                     >  2:520061/1‑235 (MQ=255)
                                                                                                                                                                              aaaGATCTCCTGAATGAAGATTCAGCGGCCCCGGTTATCCGTTTGACAAATGCGATTTTGAGTGAGCCCACAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGa                                                                         <  1:216442/250‑1 (MQ=255)
                                                                                                                                                                                 aatcccataaatgaaGTTTCGGGGGCACCGTTTATCCGCTTGATCATTGCGATTTGGAGTGAGGCCACAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCGTGATATCGCCGAAAAACGTATTCCACAGGATgg                                                                      <  2:43076/242‑1 (MQ=255)
                                                                                                                                                                                           atgaaGATTCAGCGGCACCGGTTATCAGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGATCGGCTGGATAACCCCGGAAAACGAATTCCACCGGGTGGGAGAAtaatt                                                            >  2:409695/1‑247 (MQ=255)
                                                                                                                                                                                                       cGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGCAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCggg                                                >  1:21562/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACCTATTCCACAGGATGGAAGAATTAGTTTGCGGATCGGGCGACGTaa                                        >  2:509802/1‑250 (MQ=255)
                                                                                                                                                                                                                   aTCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACt                                                                                                                                                 >  1:370808/1‑143 (MQ=255)
                                                                                                                                                                                                                      ccgctgaaTCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAg                                   <  1:19617/244‑1 (MQ=255)
                                                                                                                                                                                                                           gATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAGATGTcc                             <  1:212136/251‑1 (MQ=255)
                                                                                                                                                                                                                           gATCAATGCGATTTTGAGTGAGGCCATAAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACGGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGTGCGACGTAACATAGATGTcc                             <  1:217119/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaTGCGATTTTGAGTGAGGCCATCAAAAAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAGATGTCCGCGt                         <  1:174732/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCGATTTTGAGTGAGGCCATAAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAGATGTCCGCGTa                        <  1:166499/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        cAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGGCTTGGTCTCGCCGCAAAACGCATTCCACAGGATGGGAGAATTAGTTTGCGTATCGGGCGACGTACAAGAAATGTCCCCGTTTCCACCCCGCCGTCCATCtctg  >  2:437428/1‑248 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
ATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGCGCAGGATATCGATCAATCCGTCGATCTTCTCTCGCTTTCGGAAGAGATGCCCGCAAATGAAGATCTCCTGAATGAAGATTCAGCGGCACCGGTTATCCGCTTGATCAATGCGATTTTGAGTGAGGCCATCAAAGAAACCGCCTCTGATATCCACATTGAAACCTATGAAAAAACAATGTCGATCCGTTTTCGCATCGACGGCGTTTTGCGGACAATTTTACAGCCAAACAAAAAACTGGCGGCACTGCTTATCTCCCGAATTAAGGTCATGGCTCGTCTTGATATCGCCGAAAAACGTATTCCACAGGATGGAAGAATTAGTTTGCGTATCGGGCGACGTAACATAGATGTCCGCGTATCCACACTGCCGTCCATCTATG  >  NC_000913/3458454‑3458948

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