Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,5210TC21.2% 119.5 / 17.3 66A28A (GCT→GCCyidXputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (13/1);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 2.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTAC  >  NC_000913/3875247‑3875794
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                  
gtgtgggTTTGCCGTGCGTGTAATGTAGTACAAACGTATAATATTGTACTACAATTTAGATAACAAAAAGAATAATGCATAAAAAAAAGACATAGGACTGGAAGAAAGCTGAAAAAGGATATCAGGAGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTAGTTTAAAGGATTATTCATGTAGTGGAATTTTAAGGGATTTTTTAAGGGTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGct                                                                                                                                                                                                                                                                                    <  3:616118/269‑1 (MQ=255)
         tGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTccc                                                                                                                                                                                                                                                                          >  2:171144/1‑276 (MQ=255)
                      agttaataaaaaCGTATATTGGTTGAGTAAAAGTTAAATTACAAAAAGAACAAGGAAAAAAAAAATGACATAGGTAGGGAAGAAAGATGAAAAGGGATAGCAGGAGATAAACAGGCGGGAAAGAAGAGTATGCTTTCAACTGCCTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                             <  4:305453/267‑1 (MQ=255)
                            taCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGtt                                                                                                                                                                                                                                                        >  1:224109/1‑275 (MQ=255)
                            taCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGCCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATAACGCAAAGGGTAGtt                                                                                                                                                                                                                                                        >  1:41768/1‑275 (MQ=255)
                                     atatTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTTGTTCAGGAAAGTa                                                                                                                                                                                                                                              >  2:221271/1‑276 (MQ=255)
                                      tatTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTcc                                                                                                                                                                                                                                                                           >  1:501493/1‑246 (MQ=255)
                                                     aTTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAg                                                                                                                                                                                                                                >  2:308439/1‑274 (MQ=255)
                                                           aTCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCTCAGTCAAAAACCaa                                                                                                                                                                                                                         >  3:235658/1‑274 (MQ=255)
                                                                       acaaatcATAAAAAAATGACATGCATCGGGAAGAAATCTGAAAAGGGATATCAGGCGGTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCTGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATc                                                                                                                                                                                                             <  3:186315/269‑1 (MQ=255)
                                                                             aaaaaaaaaaTTGCATGAGGCGGGCAGAAAACAGAAAAGGGATATCAGTCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTAAAAGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGGTTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGGAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCa                                                                                                                                                                                                      <  4:221271/273‑1 (MQ=255)
                                                                              aTAAAAAAAGGACATGAGGCGGGCAGAAATCTGAAAAGGGATATCAGGCGATAAACAAGAGGGAAATAAAAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGAATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGACTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGATGTATCTCGTATGCTCTGGTTTACCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAggta                                                                                                                                                                                                      <  3:224109/275‑3 (MQ=255)
                                                                                            tGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGa                                                                                                                                                                                       >  2:105448/1‑276 (MQ=255)
                                                                                               gtgtgggaaaaaTATGGAAAGGGGTATCAGGAGCTAAACAGGGATGGAAGAAGATTATGGTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATGATGAAGTTTAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCAACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTGCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACa                                                                                                                                                                                    <  4:308439/267‑1 (MQ=255)
                                                                                                           attGAAAAGGGATATCAGGTGCAAAAAAGGAGGGAAAGAAGATTATTCTTTAAAAGGCTTAGATATTCGTTTAAAGAATTAATCATGAAGTTAAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCTCCTGGCGGTGATGCTTTCAGGCTGTATCTCGTATGCTATGGTTTCCCATACCGCAAATGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                         <  3:41768/273‑1 (MQ=255)
                                                                                                                       tatCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTCTGCCTTACAACAAAGAGTAGCTCAGGAAAGTATCAATCGCAGTCACACACCATCACTGGGCTATCGCAGGCAAAAGATAGTTATCAACTAACAGGATATGTTGTTGTCGGGGa                                                                                                                                                             >  1:590076/1‑275 (MQ=255)
                                                                                                                         tCAGGCGATAAACAGGAGGGAAAGAAGATTTTGCTTTCAAGCGCTGACCTACGTGTTTAAAAGATTAATCATGAAATTGAATTTTAAGAGATTATTTAAGGCTACCAGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAAt                                                                                                                                                           <  1:235658/275‑1 (MQ=255)
                                                                                                                          cagtggaaaaaCAGGAGGTAAAGAAGAGAATGTTTGAAACGGTTTAGCTACTCGTTTAAAGGATTAATACATAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTCATTCCCATTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                                                                                                                                         <  4:105448/268‑1 (MQ=255)
                                                                                                                             ttaagaaaCAGGGAGGTAAAGAGTTGTGTGCGTTTAACGGCTTATCTCCTCGTTTAGAGGATTAATTATGAAGTAGAATTTTAAGTGATCTTGTAAGGCTTCCGGTTTATTCCCACTGGCGCTTATGTGTTAAGGCTGTAGGTAGTATGCTCTGGGTTCGCATACAGAAAATGGTAGTTCATGAAAGTATCAATCGCAGTCAGACACCATCAATGGTCTATCGCAGGTAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                                                                                                                                         <  3:608612/269‑1 (MQ=255)
                                                                                                                                   aaCAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAACAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTg                                                                                                                                                      >  2:383820/1‑270 (MQ=255)
                                                                                                                                         gaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGt                                                                                                                                                                                                                                                            >  1:93490/1‑163 (MQ=255)
                                                                                                                                            ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAAAAAAAAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCCCAGGCAAAAGCAAGGAATGGAACAAAAGGCTCAGTTTTTGTAGGGGAATGTCTGGATAACCTTAtag                                                                                                                                         >  1:610878/1‑273 (MQ=255)
                                                                                                                                             gAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGtt                                                                                                                                                                                                                                                        <  3:93490/163‑1 (MQ=255)
                                                                                                                                              aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCaa                                                                                                                                                                                                                                                                 >  1:214244/1‑153 (MQ=255)
                                                                                                                                              aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCATTACCGCACAGGGTTGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCAACAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATGGCTGGAACACCTTATCaata                                                                                                                                      >  4:481789/1‑273 (MQ=255)
                                                                                                                                                agaagaGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa                                                                                                                                     >  1:254491/1‑275 (MQ=255)
                                                                                                                                                  aggagTAGGCTTTTAGCGGGTTAGCTACTCTTTTAAATGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGTTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATg                                                                                                                                   <  3:164895/273‑1 (MQ=255)
                                                                                                                                                  aagaGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                             <  3:214244/153‑1 (MQ=255)
                                                                                                                                                   agagTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTCTGCTAGACCGAAAAGGGTAGTTTAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCAAAAGATAGTCATGGAACAAAAGGCTATGTTTTTGTAGGGGa                                                                                                                                                             >  1:608612/1‑248 (MQ=255)
                                                                                                                                                       tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTTTGCCATAACGCAAAGGGAAGTTCAGGAAAGTATCAATCGCAGTCAGTCACCATCACTGGGCTATCGAAGGCAAAAGATAGTTATGGAATACAAGGCTAAGTTTTTGTAGGTGAATCGCGGAATTACCATATCACTGCTGGTGcc                                                                                                                              >  1:140152/1‑275 (MQ=255)
                                                                                                                                                       tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTTGAAAAAAAAAAGTCGTTCAGGAAAGTATCAATCGCAGTCAGACATCCTCACTGGGCTATCGCAGGAAAAAGATAGTAATGGAACAAAATGCTAGGTTTTTGTAGGGGAATCGATGGATTAACTTATCAATGATGGTGCCg                                                                                                                             >  2:256252/1‑276 (MQ=255)
                                                                                                                                                       tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACCGCAAAAGGTAGTTAAGAAAAGTATCAATCGCAGTCAGACATCATCACTGGGCTATCGCAGGCAAAAGGTCGTAATGGAACACAAGGCTATGGTTTTGTAGGGGAATAGCTGGATTACATTATCACTGTTgtgtccg                                                                                                                             >  2:671657/1‑270 (MQ=255)
                                                                                                                                                               cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACt                                                                                                                                       >  1:164895/1‑258 (MQ=255)
                                                                                                                                                               cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCCTTACCGCAAAGGGCAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACAt                                                                                                                      >  1:197397/1‑275 (MQ=255)
                                                                                                                                                                     cTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTaa                                                                                                                >  2:398108/1‑275 (MQ=255)
                                                                                                                                                                               cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCCATAAAAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGATGACATAGTTAAGATGCTCaa                                                                                                       >  1:574278/1‑274 (MQ=255)
                                                                                                                                                                               cGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATATTAATGGAACAAAAGGCTATGTTTTTTTAAGGGAATCGCTGGCTTCACTTTACACTGATGGTGCCGTTGACATTGTTATGactcttaat                                                                                                      >  2:324592/1‑267 (MQ=255)
                                                                                                                                                                                            aaTCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTCTGCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGAGCTATCGCAGGCAACAGATAGGAATGGAAACAAAGGCTATGTTTTTGAAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATTACATCGTTAAGCTGCTCAATGATCCAGCAATTaa                                                                                        >  2:547189/1‑276 (MQ=255)
                                                                                                                                                                                                   gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGc                                                                                   >  3:99474/1‑275 (MQ=255)
                                                                                                                                                                                                   gAAGTGGAATTTTAAGGGATTTTTTAAGGCGGCCGGTTTATTCCCACTGGCGCGGAGGCTTTCAGGCTGTATCGCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTAACGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGca                                                                                  <  2:481789/276‑1 (MQ=255)
                                                                                                                                                                                                    aaGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATTGTGCCGATGACATAGTTAAGATTCTCAATCAACCAGCACTTAACCGCca                                                                                  >  4:48803/1‑275 (MQ=255)
                                                                                                                                                                                                    aaGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATAATCCAGCACTTAACCGGc                                                                                   >  1:497437/1‑274 (MQ=255)
                                                                                                                                                                                                    aaGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGCTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaaa                                                                                                                                                                                  >  1:143628/1‑179 (MQ=255)
                                                                                                                                                                                                        tGAATTTTAAGGTTTTTTTTAAGGCGGCCGTTTTATTCCCACTGGCTCTGATGCTTTCGTGCTGTATCTCGTATGCTCTGGTTTCCGTTACCGCAAAGGGTAGTTCAGGAAAGTATCCATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                                                                                                                                              <  3:143628/179‑1 (MQ=255)
                                                                                                                                                                                                         gAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCAGGCAAAAGATCATAATGGAAAAAAAGGCCATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTTATGGTTCCGATGACATCGTTAAGTTGCTCCATTATCCAGAACTTACACGGCACCatat                                                                            >  4:369412/1‑276 (MQ=255)
                                                                                                                                                                                                                     aTTGTTTAAGGCTGCCGGTTTATTCCCACTGGCGTTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGa                                                                <  2:369412/276‑1 (MQ=255)
                                                                                                                                                                                                                                  gctggcttATTCCCTATAGCCTTGATTGTTTCAGGCTGTATCTCGTATGGTATGGTTTCCCAACCCGTAAAGGGTAGTTCAGGAAAGTATCAATCGCATTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAt                                                     <  4:563788/268‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGTAATCGCTGGATTACATTAGCACTGATGGTGCCGATGACATCGTTCAGATGCTCAATAATACAGACATTAACCGGAACAATCTTCATGTGCACTATGACGTAAAAttt                                                  >  3:519304/1‑275 (MQ=255)
                                                                                                                                                                                                                                   ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAAGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATAGTGAAGATGACAAATGAGCAAGCAATTAACCGGCACAGTATTCAGGTTGACGCTGACGGATGATTTg                                                 >  4:706812/1‑276 (MQ=255)
                                                                                                                                                                                                                                   ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCGGGTTGCCGATGACGCAAGAtt                                                    >  2:563788/1‑274 (MQ=255)
                                                                                                                                                                                                                                          aTTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGAGGGCTTACCTTATCACTGATTGTGCCGATGACATCGTTAAGATGCTCAATGAGCCAGCACTTAACCGGCACACTATTCAGGTTGTCGATGACGCAAGATTGTGtttaat                                           >  4:667711/1‑273 (MQ=255)
                                                                                                                                                                                                                                             cccACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                                                                                                                                                   >  1:518931/1‑137 (MQ=255)
                                                                                                                                                                                                                                                aCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaaaa                                                                                                                                                                                                   >  2:422019/1‑118 (MQ=255)
                                                                                                                                                                                                                                                 cTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAgg                                                                                                                                                                               <  3:518931/137‑1 (MQ=255)
                                                                                                                                                                                                                                                   tgcgcTGATGCGTTCAGGCTGTATCTCTTATGCTCTGATTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAt                                                                                                                                                                                                <  4:422019/117‑1 (MQ=255)
                                                                                                                                                                                                                                                     cgcTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACCAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCAAGCAATTAACCGGAACAATATTCAGGTTGCCGATGACGCAAGATTTTTTTTAAATTCGGGGtaa                                >  2:348836/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCt                                                                                                                                                       >  3:535114/1‑153 (MQ=255)
                                                                                                                                                                                                                                                         gATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTGCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTAAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGa                            >  1:462839/1‑275 (MQ=255)
                                                                                                                                                                                                                                                         gATGCTTTCAGGCTGTATCTCGTATGCCCTGGTCTCCCATACCGCAAAGGGTAGCTCAGGAAAGTATCCATCGCAGTCCGACACCATCACTGGGCTATCGCAGGCAAAAGACAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTg                                                                                                                                                      >  2:677832/1‑154 (MQ=255)
                                                                                                                                                                                                                                                             cTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                                                                                                                                                   <  1:535114/153‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cggCGGTCTCTCGTATCGTCCGTTTTCCCATACAGCAAAGGGTAGTTCATGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCGATCGAAGGCAAAAGATAGTAGTGGAACAAGAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCgg                  <  2:475675/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACaa                                                                               >  3:63003/1‑216 (MQ=255)
                                                                                                                                                                                                                                                                     gtGTATTTCGTAGGCCTTGGTTTTCCCAACCGAAAAAGGTAGTTCAAGAAAGTATCTATCGCAGTCAGACAACCACACGTGGCTATCGCAAGCAAAAAATAGTAATGGAACAAAATGCTATGTTTTTGTAGGGGAATCGCTGGATTACCGTATCAGTGATGGTGCCGATTACATCGTGAAGATGATCAATGGTCCAGCACTTAACCGGCACAATAGTCAGGTTGCCGATGACGCAAGATTTGGTTTAAATGCGGGGAAAAAGAAATTTACCGGc                 <  3:497437/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ctgcttctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAATTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGca                <  2:706812/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATAtt                                                                           <  1:63003/216‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tATCCCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCAATGTTTTTGTAGGGGAATCGCTGGATCACCTTATCACTGATGGTGTCGCCGACATCGTTAAGTTGCTCAAGTATCCAGCCCTTCACCGCCACAATACTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATTCGGGGAACAAGAAATTGAGAGGcacatt            >  4:475675/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            tcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAata           >  4:208597/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                                   ctctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtac  <  2:48803/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                  
GTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTAC  >  NC_000913/3875247‑3875794

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