Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,5210TC23.6% 80.0 / 13.4 55A28A (GCT→GCCyidXputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (25/17);  new base C (13/0);  total (38/17)
Fisher's exact test for biased strand distribution p-value = 5.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCT  >  NC_000913/3875250‑3875790
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
cgcgTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTg                                                                                                                                                                                                                                                                            >  2:388713/1‑275 (MQ=255)
                              aCTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGa                                                                                                                                                                                                                                              >  2:31278/1‑275 (MQ=255)
                              aCTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGa                                                                                                                                                                                                                                              >  1:390467/1‑275 (MQ=255)
                                   tatTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTTTGGTCTCCCATACCGCAAAGGGTAGTTCAGGaaatga                                                                                                                                                                                                                                         >  1:355243/1‑272 (MQ=255)
                                        ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGCAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATTTCGTATGCCGTTGTTTCCCTTAACGCAAAAGGTACGTCAGGAAAGTATaaac                                                                                                                                                                                                                                    >  1:525738/1‑274 (MQ=255)
                                        ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTCTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATTCTTTCAGGTCGTATCTCGTATGCTCTTGTTTCCCATACCGCAAAGGGAAGTTCATGAAATTATCAAt                                                                                                                                                                                                                                    >  1:137347/1‑275 (MQ=255)
                                          gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATTCGTCGGGCAGAAATCTGACAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCTTTTAAAGGATTAATCCTGCAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCCTACCTCacctggtcgttc                                                                                                                                                                                                                                                  >  1:205360/1‑248 (MQ=255)
                                          gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCTTACCGCAAAGGGGAGTTCAGGAAAGTATCa                                                                                                                                                                                                                                      >  1:599583/1‑271 (MQ=255)
                                          gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATAGCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTTCCATACCGCAAAGGGTCGTTCAGGAAAGTACCAATc                                                                                                                                                                                                                                   >  1:212389/1‑274 (MQ=255)
                                            actacAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGc                                                                                                                                                                                                                                 >  2:190382/1‑274 (MQ=255)
                                              taaaaTTTAGATCACAAAAAGAACAATGCAGAAAAAAATGACATGCGTCTGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGAGAAAGAAGAGTATGCTTTCAAAGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCg                                                                                                                                                                                                                                  <  3:599583/268‑1 (MQ=255)
                                                cAATTTAGAACACAGAAATAACAATGAATAAAAAAAAAAAATGTGTATGGCAGAATTATGAAAAGGAATAACAGGCACTAAACAACAAGAAAAAAAGATGAAGATTTAAACGGATTAGCTACTCGTTTAAAGGATTAATAGTGAAGTGGAGTTTTAAGTGAATTGTTAAGACGGCCGGTTTAATCCCACTGTCGCTGATGATTTCAGGCTGTATCTCGTATGTTTTGGTTTCCCCTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCa                                                                                                                                                                                                                            <  3:313715/275‑1 (MQ=255)
                                                     taaataacaAAAAGAAAGATGCATAAAAAAAAGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAATAGTATGCTTTCAACGGCTTAGCTACTCGTTTCAAGGATTAATCATGAAGTTGAGTTTTAAGGGATTATTTAAGGCTGCTGTTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATCCCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGacac                                                                                                                                                                                                                       <  3:365909/269‑1 (MQ=255)
                                                     tAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCCTTAAAGAAAAGGGTAGTGTAGTAAATTATCAATCGCAGTcgcaca                                                                                                                                                                                                                        >  2:419544/1‑269 (MQ=255)
                                                        aTCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTGCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAt                                                                                                                                                                                                                    >  2:93497/1‑275 (MQ=255)
                                                          cagaaaaaGAACAAAGAATAAAAAAATGACATGGTGCGGGCAGAAATATGGAAAGGGATATCAGGCGATAAACAGGAGTGAAAGAAGTGTATGCTTTCAACGGCTTAGCTACTCGTGTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCa                                                                                                                                                                                                                  <  3:11681/272‑1 (MQ=255)
                                                                       aTGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCCGGAAAGTACCAATCGCAGTCAGAAACGATCAATGGCATTTCGCa                                                                                                                                                                                                     >  3:110716/1‑275 (MQ=255)
                                                                       aTGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCgag                                                                                                                                                                                                     >  2:140686/1‑273 (MQ=255)
                                                                          cATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTCTGCCAGAAAGCAAAGGGTGGTTCAGGTAAGTATCAATTGAAGTCAGAAACCATCACTGGGAAATAGCAggg                                                                                                                                                                                                  >  2:487315/1‑274 (MQ=255)
                                                                             aaaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTACACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGCCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCAGTATCTCGTATGCTCTGGTTTCCCCTCCCGCAATGGATAGTTCGGGAACGTATCAATCGAATTCCGACACCATCATGGGTCTATCGTAGGCaaa                                                                                                                                                                                               >  3:78852/1‑275 (MQ=255)
                                                                             aaaaaaaTGACAGGGGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATAAAGTTGAATGTTAAGAGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaaa                                                                                                                                                                                               <  2:192600/275‑1 (MQ=255)
                                                                                aaaaTGAAATGCGACGGGCATAAAGTTGTAAAGAAAGATCAGGGGCTAAAGAGACGGAAAAAAAGAGTATCCTTTCAAAGGTTTAGCGAGTAGTTTAAACGATTAATGATGACGTTGAAATGAAAGGTTTTTTGCAGGAATACGATGTAATTTACTCTTGCGCCGATGCTGTCTGCATGTAGCTGGTATGCTCTTTGTTCCGAGACCGCAACGGGTAATTCAGTAATGTATCAATCAGAGTCAGACACCATCACTGGGCTATCGCAGGCAAgaga                                                                                                                                                                                            <  3:390467/275‑1 (MQ=255)
                                                                                  aaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCCAAAGATa                                                                                                                                                                                          >  1:303196/1‑275 (MQ=255)
                                                                                                          tGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCTTTAAAAAAAAAGGTAGTTCAGGAAAGTATAAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATCGTAATGGAACAAAAGGCTATGtt                                                                                                                                                                   >  2:76528/1‑274 (MQ=255)
                                                                                                                    tatCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCTTACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACAATCACTGGGCTATCGCAGGCAAACGATAGTAATGGAACAACAGGCTATGTTTTTGTAggg                                                                                                                                                          >  2:258485/1‑273 (MQ=255)
                                                                                                                    tatCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGAAGTCAGACACCCTCACTGTGACCTCGCAGGCAAAAGCTAGTAATGCAACAAAAGGATATGATTTTGtgtgt                                                                                                                                                          >  2:329814/1‑269 (MQ=255)
                                                                                                                      tCAGGCGATAAACAGACGCGAAAGCAGAGTATTCTTTCAACGGCTTAACTACTCGTTTAAAAGATTAATCATGAAGTTGAATTTTAAGGTATTTTTTAAGGCTGCCGGTTTATTCCCCCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAGGGGTAGTTCAGGAAAGTATCAATCGCAGTCAAACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAAt                                                                                                                                                      <  3:83109/275‑1 (MQ=255)
                                                                                                                       ctggCGCTAAACAGGAGGGGAATAAGGGTATGGTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTATAAGGTATTTTTTAAGGCTGCAGGTTTATTCCCGCTGGCGCTGATGCTTTCATGCGGTATCTCGTATGCTCTGGTTGCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATc                                                                                                                                                     <  1:110716/273‑1 (MQ=255)
                                                                                                                        agagtcTAAACAGTAGGAAAATATGAGTATACTTTAAACGGCTTAGCTACTCGTTTTAAGGATTAATCTTGAAGTTGAATTTTAATGGATTTTTTAAGGCTGCCGTTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                                                                                                                                    <  4:31278/270‑1 (MQ=255)
                                                                                                                          gcgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATCCCGCAACGGGTCGTTCAGGCAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCCAGGTTTTTGTAGTGGatgcgc                                                                                                                                                   >  2:275966/1‑269 (MQ=255)
                                                                                                                           cgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTTTAACGCAAAGGGTAGTTCAGGAAAGTATCAATTGCAGTCAGACACCATCACTGGGCTATCGTAGGCAAAAGATAGTAATGGAACAAAAGTCTATGTTTTTGTATGGGTATCGc                                                                                                                                                   >  2:41451/1‑273 (MQ=255)
                                                                                                                           cgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTa                                                                                                                                                                                                          >  1:294140/1‑218 (MQ=255)
                                                                                                                               aaaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGc                                                                                                                                                                                                      <  3:294140/218‑1 (MQ=255)
                                                                                                                                accgggatgGAAAGAATGGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAAGTTTAAGGTATTTTGTAAAGCTGCCGGTTTATTCACACTGGCGCTGATGCTTTCAGGCTGTAACTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGCCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                                                                                                                                              <  3:355243/269‑1 (MQ=255)
                                                                                                                                            agaagaGTATGTGTTCAACGGGTTAGCGAATCGTTTAAAGGATGAATCATAAAGTTGAATTTTAAGGGATTTTTTAAGGCGGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTAGCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa                                                                                                                                <  4:140686/275‑1 (MQ=255)
                                                                                                                                                   taagCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTAGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAACATAGTCATGGAACAAAAGGCTAGGTTTTTGTAGGGGACTCGCTGGATTACCTTCTCACTGATGGTGcc                                                                                                                         >  3:133307/4‑275 (MQ=255)
                                                                                                                                                   tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAAGGGAACAAAAGGCTATGTTTTGTGAGGGGAATCGGTGGATTTCCTTTTTCCTGATTGTGcc                                                                                                                         >  4:121795/1‑275 (MQ=255)
                                                                                                                                                   tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTTGAAAACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAACGGCTATGTTTTTGTAGGGGAATAGCTGGATTACCTTATCACTGATGTTGcc                                                                                                                         >  2:265315/1‑275 (MQ=255)
                                                                                                                                                          acaatggtgtAGCTACACGTTGAAAGGATTAATCATGAGGTTTAATTGTCAGTGATGTTTTAAGGGTGCCGGTTTATTACCAATGGGGGTTGTGTTTTAAGGCTGTATATCGTGTGATGTGGTATCTCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTAAGACACCATCACTGGGCTATTTCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGTGGAAGCGCTGGATTACCTTATCACTGATGGTGCCGATGACa                                                                                                                  <  4:258485/266‑1 (MQ=255)
                                                                                                                                                           cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCACAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCt                                                                                                                                                  >  2:31804/1‑242 (MQ=255)
                                                                                                                                                                  ttggatACTCGTTTAAAGGATTAATCATGAAGTTGAATGTTAAGGTATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAg                                                                                                          <  1:133307/270‑1 (MQ=255)
                                                                                                                                                                  ttAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGTCAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGAAAAAGATCGTAATGGAAAAAAAGGCTATGTTTTTGTAGGGGCATCGCTGGATTACCTTATCACTGATGGTGCCTATGACCTCTTTaa                                                                                                           >  4:331849/1‑274 (MQ=255)
                                                                                                                                                                     gatattattttAAAGAGTTTATCATGAAGTTGAATTTTAAGGGATTTGTTAAGGCTGCAGGGTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCCAACCGCAAAGGGTAGTTCAGGAAAGTATCGATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATg                                                                                                       <  3:303196/267‑1 (MQ=255)
                                                                                                                                                                       ttatCGTTGAAAGTATTCATCATGAAGTTGTAATTTGAGGGAGTTGTTAAGGCTGCGGGTTTATTCCCACCGTGGCTGAAGCTGCCAAGATGTATCACGGATGCTATGCGTTCCCAAAGCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAATCAGACACCATGATTGGGCTAACGCAGGCAAAAGATAGTAATGGAACAAAAGGGTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGGTAAGTTGCt                                                                                                     <  4:21763/272‑1 (MQ=255)
                                                                                                                                                                              ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTCGTTTGCCATACCACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACAGGGCTACCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTATTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGACGACATCGTTAAGATTCTCAATGAt                                                                                              >  1:372457/1‑275 (MQ=255)
                                                                                                                                                                              ttAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTCTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAAAGTAATGGAACAAAAGGCTGTGGTTTTGTAGGGGAATAGCTGGATTACCTTATCACGGATGTTCCCGATGAGAGCGTTCAGCTTTTCAATGAt                                                                                              >  1:352835/1‑275 (MQ=255)
                                                                                                                                                                                       tAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCCGGAAAGTATCAATCGCAGTCAGACACTATAACTGGGCTATCGCAGGACAAAGTAAGTACTGGAACAAAAGGATATGTTTTCTTATGGGAATCGCTGGATTTCCATATCACTGATGGAGCAGATGACATCTTGAAGATGATAAATGacacagagatt                                                                                     >  4:110163/1‑265 (MQ=255)
                                                                                                                                                                                           cATGAATTTGTATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGTGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAAc                                                                                  <  2:121795/274‑1 (MQ=255)
                                                                                                                                                                                                    gAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGCCTCGCTGTATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCAAGCACTTAACCGGCACAata                                                                        >  2:488586/1‑275 (MQ=255)
                                                                                                                                                                                                                                                      aTGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTTTACCAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGATCAATTATCCAGCACTTAACCGGACCCATATTCAAGTTTCCGATGACGCAAGTATTGTGTTAAATACGGGGAAAAAg                         >  2:350342/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                    tATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGTTAGTAATGGAACCAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATTACTGATGGTTCCGATGACATCGTTACGATTCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGAGGACCACAGATTTGCTTTACATGCGGGGAAACAGAAATTTACTGGcacca         >  2:372871/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                    tATCTCGTATGCCCTGGTCTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATAGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGACGACGCAAGACTTGTTTCAAATGCGGGGAAAAAGAAATTTACAGGcac           >  2:252816/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                     aTCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGCTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCATGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATg                                                                                                >  1:428481/1‑187 (MQ=255)
                                                                                                                                                                                                                                                                           tATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATAACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGAACAAGCACTTAACCGGAACAATATTAAGGTTGCCGAAGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGACATATACAGGCACaaacacact  >  4:413548/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                            aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACaa                                                                          >  1:67/1‑202 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                             
CGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCT  >  NC_000913/3875250‑3875790

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: