Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,5210TC24.4% 68.0 / 16.7 42A28A (GCT→GCCyidXputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/18);  new base C (10/0);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 2.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACT  >  NC_000913/3875266‑3875795
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
gtAATGTAGTAAAAACTTATATTGTTTTACTACAATTTAGAACACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAAGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                  <  2:358207/274‑1 (MQ=255)
    tGTAGTAAAAACTTATATTGTTGGACTACAATTTATATCAAAAAAAGAACAAAGCATAAAAAAATGACATTCGTCGGTCAGAAATGTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTGTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGTTTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                             <  2:868768/275‑1 (MQ=255)
         taCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGTGACAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTATAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGttt                                                                                                                                                                                                                                                                             >  1:671071/1‑254 (MQ=255)
                   tatTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCTGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTa                                                                                                                                                                                                                                              >  1:128402/1‑275 (MQ=255)
                        ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTTTGCGCTTGTTTCCTTTACCGCAAATGGTAGCTACGGAAAGTATCat                                                                                                                                                                                                                                          >  1:245103/1‑273 (MQ=255)
                        ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCaa                                                                                                                                                                                                                                                                 >  2:902332/1‑251 (MQ=255)
                          gtACGACAATTTAGATCACAAAAAGAACAATTCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAAAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCg                                                                                                                                                                                                                                       <  2:41970/275‑1 (MQ=255)
                            actacAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                             <  1:902333/251‑1 (MQ=255)
                                cAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAAGTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCTCTACAAAAAAGGCAAGGTCAGGAAAGTATCAAACGCAgccg                                                                                                                                                                                                                                 >  1:945618/1‑272 (MQ=255)
                                  aTTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTAGTATGCCGTCGTTTGCTTGACAGAAAAGAGTGGTTCAGGGAAGTATAAATGGCAGTCAGa                                                                                                                                                                                                                               >  1:329772/1‑275 (MQ=255)
                                        aTCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAATTAATACAGGAAAGTATCAATCGCAGTCAGACACCAt                                                                                                                                                                                                                         >  1:176210/1‑275 (MQ=255)
                                            caAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTGCTTTACCGAAAAGGGTAGTTCAGGAAAGTATAAATCGCAGTCAGACAACCTCACt                                                                                                                                                                                                                     >  1:864495/1‑275 (MQ=255)
                                             aaaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCGGGAAAGTATCAATCGCAGTCAGACACCATAACTg                                                                                                                                                                                                                    >  1:30977/1‑275 (MQ=255)
                                             aaaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTCTGCCTTACAACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTg                                                                                                                                                                                                                    >  1:1035521/1‑275 (MQ=255)
                                                           aTAAAAAAACGGCCAGGGTCGGGAAGAGATCTGCAAAGGGTGTTCAGGCGCTAAACAGGAGTGTAGGAGGAGTATGGTTTGAACGGCTGGGCTACTCGGTTAAAGGAGTAATCATGAGGTGTGATGTTGAGGGATTTTGTAAGGCTGCAGGTTTATTCGCACTGGAGCCGAAGTTTTCAGGCTGTATCTGGTATGCTCTGGTTACCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCCGTCAGACACCATCACTGGGCTATAGCAGGCa                                                                                                                                                                                                      <  2:992147/275‑1 (MQ=255)
                                                            tAAAAAAATGACATGGGTGGGGCAGAAATCTTAAAAGGGATATAAGGCGATAAACAGGGGGGAAAGAAGAGTATGCTTTCAACGGCTTAACTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAATGGATTTTGTAAGGCTGCCGGTTTATTCCCACTGGCGCTAATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaa                                                                                                                                                                                                     <  2:295394/275‑1 (MQ=255)
                                                                        atgctgcGGGGAGAAATCTGAAAAGGGATATCAGGCGCTAAAAAGGAGGGAAAGAAGAGTATGCTTTAAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAAGTTTAAGGGATTTGTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATg                                                                                                                                                                                         <  2:30977/269‑1 (MQ=255)
                                                                                                 ggATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgg                                                                                                                                                                >  1:315620/1‑275 (MQ=255)
                                                                                                        aGGCGCTACACAGGAGGGAAAGAAGAGTATGATTTTAACGGCTTAGCGACTAGTTTAAAGGATTAATCATGAAGTTGAATTTTAATGGATTTTTTAAGGCTGCCGGGTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGAATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCATTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                                                                                                                                         <  2:775484/275‑1 (MQ=255)
                                                                                                              tagagtatagtgAAAGAAGAGTATGCTTTGAACGGCTTAGCTACTCGTTTAAACGATTAATCATGAAGTTGAATTTTAAGGGATGTTTTAAAGCTGCCGGTTTATTCCTACTGGCTCTGATGCGTTCAGGCTGTATCTCGGATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                                                                                                                                                   <  2:57040/264‑1 (MQ=255)
                                                                                                                  cAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAcc                                                                                                                                               <  2:315620/275‑1 (MQ=255)
                                                                                                                    ggaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCTGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                                                                                                                                             <  2:128402/275‑1 (MQ=255)
                                                                                                                    ggaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                                                                                                                                             <  2:760473/275‑1 (MQ=255)
                                                                                                                    ggaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTCTGCTTTAACAAAAAAGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAAGCTATGTTTTTGTAGGGGAATCGCTGGATTAACt                                                                                                                                              >  1:210997/1‑274 (MQ=255)
                                                                                                                        ggAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTTAAAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGCTCGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATAGCTGGATTACCTTATCa                                                                                                                                         >  1:1030277/1‑275 (MQ=255)
                                                                                                                                  gTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGc                                                                                                                               >  2:824844/1‑275 (MQ=255)
                                                                                                                                  gTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGc                                                                                                                               >  1:191413/1‑275 (MQ=255)
                                                                                                                                                      cTACACGGTTAAAGGATTAAGCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCATGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATg                                                                                                            <  2:191413/274‑1 (MQ=255)
                                                                                                                                                                      ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAAATCGTTAAGATGCTCAATGATCCAGCACt                                                                                           >  1:279387/1‑275 (MQ=255)
                                                                                                                                                                          taaTGAAGTTGAATTTTAAGGGTGTTTGTGAGGTTGCCGGTTTATTACCACTAGCTCGTATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGGAAAGGGTAGTTCAGGAAAGGATCAATCGCAGTCAGACAACATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATTCAGGACTTAAc                                                                                       <  2:347486/273‑1 (MQ=255)
                                                                                                                                                                              gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTTTACAAAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATAGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGACTCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTACAGTGCTCAATGGTCCAGCAATAAACCGGc                                                                                   >  1:962522/1‑275 (MQ=255)
                                                                                                                                                                               aaGTTTAAGTTGAAGGGATTTTTTAATGCTTGCGGTTTAGTGCCACGCGCGATGATGCTTTCAGGCTGTATCTCGTATCATCTGCATTCCCATAACGCAAAGGGTAGTTCAGGAAAGTATCAAGCGCAATCAGACACCATCACTGGGCTATAGAAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGTGGAATAGCTGGATTACCTTATCACTGATGGTGCCGATTACATCGTTAAGATGCTCAATGATCCAGCTCTTAACCGGca                                                                                  <  2:1035520/275‑1 (MQ=255)
                                                                                                                                                                                                              ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTGCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTTGGCTATCGCAGGCAAAAGATAGTACTGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTACGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCCGGTTGCCGATGAAGCAAGAttt                                                   >  1:85393/1‑275 (MQ=255)
                                                                                                                                                                                                              ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTCCTTTACAACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGAttt                                                   >  1:369401/1‑275 (MQ=255)
                                                                                                                                                                                                                   ttATTCCCACTGTCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGtttt                                              <  1:824845/275‑1 (MQ=255)
                                                                                                                                                                                                                           aCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCAc                                                                                                                                                                                                                      >  2:779997/1‑99 (MQ=255)
                                                                                                                                                                                                                            cTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTGCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATAGTTAAGATGCTCAATGCTCCAGCACTTAACCGGCACAATATTCAGGTTGCCGCTGAAGCCAGATTGGTTTTAAATGCgg                                      >  1:773883/1‑274 (MQ=255)
                                                                                                                                                                                                                               gcgcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTggg                                                                                                                                                                                                                  <  1:779997/99‑1 (MQ=255)
                                                                                                                                                                                                                                                   tATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTTGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCAAGCACTTAACCGGCACACTATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGAGGGGAAAAAGAAATTTCCCggaccaa              >  2:332071/1‑270 (MQ=255)
                                                                                                                                                                                                                                                    aTCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCAATTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGAACaa              >  1:678955/1‑274 (MQ=255)
                                                                                                                                                                                                                                                     gctcGTATGCTCTGGTTTCCAATACCGGAAAGGGTAGTTCAGGAAAGAAACAATCGCAGTCAAAAACAATCACTGGGCTATCGCAGCCAAAAGATAGTAATGGAACAAAAGGCGATGTTATTGTAGGGGAATCGCTGGGTTACCTTATCGCTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATAATCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGAGGGGAAATAGAAATTTTCCTGCACAata            <  2:876299/274‑1 (MQ=255)
                                                                                                                                                                                                                                                             gCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTta    >  2:99321/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               tctGGTTTCCCATACCGCAAAGGGGAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCAAGCACTTAACCGGCACAATATTCAGGTTGCCGATGAAGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATGTACCGGCCCAATCTCGCTTtact  >  2:92489/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
GTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACT  >  NC_000913/3875266‑3875795

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