Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,5210TC27.3% 93.4 / 21.1 67A28A (GCT→GCCyidXputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (28/20);  new base C (18/0);  total (47/20)
Fisher's exact test for biased strand distribution p-value = 6.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACT  >  NC_000913/3875263‑3875795
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
cgtgTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTAGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCGATACCGc                                                                                                                                                                                                                                                                   >  2:9380/1‑276 (MQ=255)
      atttagtaCAAACTTATAATTTTGTACTAAAATTTAAATCACAAAAAGAACAAATCAAAAAAAAATGAAAGGCGTCGGGAAAAAAACTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAATAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                             <  4:431816/273‑1 (MQ=255)
           gtaCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCAGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTCTCCCTTAACGCAAAGTGTatt                                                                                                                                                                                                                                                         >  1:445108/1‑273 (MQ=255)
               aaaaTTAAATTTGTGTACTACAATTTAAATCAAAAAAAGAACAATGCATAAAAAAAAGACATGCGTCGGGCAGAAATATGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAgg                                                                                                                                                                                                                                                    <  4:28175/276‑1 (MQ=255)
                          gttggAATAAAATTTAGATCACAAAAAGAACAATTCATAAAAAAATTACATGCGTCGGCAATAAATCTGAAAAGGGATATCATGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGACTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAAt                                                                                                                                                                                                                                         <  4:297244/276‑1 (MQ=255)
                           ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGTTTTCCTATACAGCAAAGGGTATTGCAGGAAAGTATAAATc                                                                                                                                                                                                                                        >  2:605827/1‑276 (MQ=255)
                            tgtACTACAAATTAGATAACACAAAGAACAATGCACAAAAAAATGACAGGGGGTGGGTAGACAACTGAAAAGGTATATCAGGAGCTAAGCAGGAGGGAAAGAGGAGTATGAGTTCAACGGGTTAGCTACTAGTTTAATGGATTAATCATGAAGTTGAATTTTAAGGGATTTTATAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGTCTGTATCTCGTATGCTCTGGTTTCCCAGACCGCAAAGGTTAGTTCCGGAAAGTATCAATCg                                                                                                                                                                                                                                       <  4:582382/276‑1 (MQ=255)
                             gtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATc                                                                                                                                                                                                                                        >  1:591046/1‑274 (MQ=255)
                                   cAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAg                                                                                                                                                                                                                                >  2:236316/1‑276 (MQ=255)
                                   cAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTTTCCTTGACAAAAAAGGATAGTACAGGAAAGTAACAATCGCAGTc                                                                                                                                                                                                                                  >  1:726837/1‑274 (MQ=255)
                                   cAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCAACTGGCGCTGAGGCTTTCAGTCTGTATCTCGTATGCCCTGTTTTGCTTTCCCGAAAAGGGTAGTACAGTAAAGTATCAAttggagcaa                                                                                                                                                                                                                                 >  1:666558/1‑267 (MQ=255)
                                      tttAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTAAGac                                                                                                                                                                                                                              >  1:553981/1‑275 (MQ=255)
                                        tAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTGCCTTACAGAAAAGGGTAGTTCAGGAAAGTATCTATCGCATTCATACAcc                                                                                                                                                                                                                           >  2:202646/1‑276 (MQ=255)
                                           aTCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCGCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCTGAGTCAGACACCATc                                                                                                                                                                                                                        >  4:470271/1‑276 (MQ=255)
                                                  aaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTggg                                                                                                                                                                                                                  <  1:491789/275‑1 (MQ=255)
                                                          aTGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTGCCATACCTAAAAGTGAAGTTCAGGAAAGTTTCAATCGCATCCATACACCATCACTGGGcactccga                                                                                                                                                                                                          >  1:498637/1‑268 (MQ=255)
                                                             cATTAAAAAATGAAATGCGTAAGGCAGAAATCTGAAAAGGTATATCAGGCACTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCAGGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCa                                                                                                                                                                                                      <  4:162702/276‑1 (MQ=255)
                                                                aaaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATCCCGCAAAGGGTAGTTCCGGAAAGTATCAAGCGCAGTCAGACACCATCCATGGGGTACCGAAGGCaaa                                                                                                                                                                                                    >  3:472046/1‑275 (MQ=255)
                                                                                tCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTGTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGCTCGTAATGGAACaa                                                                                                                                                                                   >  2:90642/1‑276 (MQ=255)
                                                                                tCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCACAACTGGGCTATCGCAGGCCAAAGATAGTAATGGAACaa                                                                                                                                                                                   >  2:416592/1‑276 (MQ=255)
                                                                                tCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCTTACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATATTAATGGAACa                                                                                                                                                                                    >  1:631207/1‑275 (MQ=255)
                                                                                  gggCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCTCTACCACAAAGGGTAGTTCAGTAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATAGCAGGCAAAAGATAGTAATGGAACaaa                                                                                                                                                                                  >  1:478646/1‑275 (MQ=255)
                                                                                   ggCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATAGCAGTCAGACACCAGCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaaaa                                                                                                                                                                                 >  1:89997/1‑275 (MQ=255)
                                                                                   ggCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTCTCCCTTACAGCAAAGAGTATTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCCGGCAACAGTTAGTACTGGAACAAAAg                                                                                                                                                                                >  2:600904/1‑276 (MQ=255)
                                                                                        aaaTCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCTGGAAAGTATCAATCGCAGTCATACCCCATCACTGGTCTTTCGCAGGCACAAGATAGTAATGGAACAAAAGGCTAt                                                                                                                                                                           >  4:64099/1‑276 (MQ=255)
                                                                                                              cgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCTATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGATATTGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCgt                                                                                                                                                        >  1:160984/1‑272 (MQ=255)
                                                                                                                cTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAGGATAGTAATGGAACAAACGGCTATGTTTTTGTTGGGGAATCGCTGGAt                                                                                                                                                   >  4:152827/1‑276 (MQ=255)
                                                                                                                   aagaggaggGAAAGGAGTGTATGCTTTACACGGCTTAAGTACTGGTTATAAGGATTAATCATGGAGTTGAATTTTAAGGGATTTTATAGGGCTGCCGTTTTATTCCGACTGGCGCTGATGCTTTCAGGCTGTATCGCGTATGCTCTGGGGTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                                                                                                                                                   <  3:45800/270‑1 (MQ=255)
                                                                                                                     cAGGAGGGAAAGAAGGTGATTCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATTAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTAATGCTTTCAGACTGTATCTCGAATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAc                                                                                                                                                <  3:89997/274‑1 (MQ=255)
                                                                                                                             aaaGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTACGTTTTTGTATGGGAATCGCTGGATTACCTTAtgcatg                                                                                                                                      >  4:66951/1‑271 (MQ=255)
                                                                                                                                      tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAcc                                                                                                                                               >  1:421290/1‑258 (MQ=255)
                                                                                                                                      tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCg                                                                                                                             >  2:126712/1‑276 (MQ=255)
                                                                                                                                          cTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTAt                                                                                                                                           <  3:421290/258‑1 (MQ=255)
                                                                                                                                           tttCAAAGGCTTAGATACTAGTTTAAAAGATTAATAATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTTCCGGTTTATTCCCACTGTCGCTGATGCGTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                         <  3:397291/275‑1 (MQ=255)
                                                                                                                                           tgtaaaCGTCTTATCTACTTGTTTAAAGGATGAAGCATGAAGTTTAATGTTAAGGGATTTTTTAAGGCTGCCGGTGTATTCCCACTGTCGCTGATGCTTTCAGGCTGTATCTCGGATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                         <  4:90642/271‑1 (MQ=255)
                                                                                                                                            ttCAACTGCTTAGCTACTTGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTTCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                         <  3:134357/274‑1 (MQ=255)
                                                                                                                                            ttCAACGGGTTAGATTATCGTTTAAAGAATTAATCATGAAGTTGAATTTTAAGGGATTGTTTAAGGCTGCCGATTTATTGCCCCATGGGCTGAGGCTTTCAGGAGATAGCTCGGATGCTCTGGATACCCATAATGTAAAGGGTAGTTCAGGAAAATATCAAACGCAGTCAGACACCAGCACTGGGCTATAGCAGGCAAAAGATAGTAATGGTACAAAAGGCTATGTTTTTGTATGTGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                         <  3:554340/274‑1 (MQ=255)
                                                                                                                                              cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAcc                                                                                                                                                                                                                           >  1:407082/1‑174 (MQ=255)
                                                                                                                                                  ggCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCa                                                                                                                                                                                                                       <  3:407082/174‑1 (MQ=255)
                                                                                                                                                    cTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACAGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTTGATTACCTTATCACCGATGGTGCCGATGACATCGTTaa                                                                                                                >  2:432156/1‑275 (MQ=255)
                                                                                                                                                    cTTAGCTACACGTATAAAGGATTAATCAGGAAGTTGAAGTTTAAGGGATTTTTTATGGCTGCCGGTTTATTACGACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTa                                                                                                                 <  2:64099/274‑1 (MQ=255)
                                                                                                                                                     ctAGCTAATCGTTTAAAGGATTAATCATGAAGTTGAATGTTAAGGGATTTGTTAAGGATTCTGTTTTATTCCCACTGGCGATGATGCTTTCAGGCTGTATCTCGTATGATCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAg                                                                                                               <  3:591046/274‑1 (MQ=255)
                                                                                                                                                                  tAAAGGATTAAGCATGAAGTTGTATTTGAAGGGATTGTTTAAGGCTGCGGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTAAGGAAAGTATCAATCGCAGTCAGACACAAGCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATc                                                                                                  <  2:470271/275‑1 (MQ=255)
                                                                                                                                                                         ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCCGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTACCGATGACATCGTTAAGATGCTAAATTATCCAGAACtt                                                                                          >  4:564500/1‑276 (MQ=255)
                                                                                                                                                                                 gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATAGCTGGCTTACCTTAACACTGATGGTTCCGCAGCCAACGTTAAGATGCTCAATGCTCCCATACCTAACCGGcc                                                                                  >  4:110206/1‑275 (MQ=255)
                                                                                                                                                                                 gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGca                                                                                  >  2:391491/1‑276 (MQ=255)
                                                                                                                                                                                  aaGTTGCAATGTAAGGGATTTTTTCAGGCTCACGGTTTATTACCACTGGGGATGATGCTTTCAGGCTGGATCACGTATGCTCTGGTTTCCAAGACCACAAAGGGTAGGTCAGGAAAGTATCAATCGCAGTCAGACGCCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGcac                                                                                 <  4:416592/276‑1 (MQ=255)
                                                                                                                                                                                       gAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAGAAAAGATTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATAACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAACGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTACGGAGGTCAATGATCCATCACTTACCCGGCACAata                                                                             >  1:272731/1‑275 (MQ=255)
                                                                                                                                                                                       gAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCCCAatat                                                                            >  2:525763/1‑276 (MQ=255)
                                                                                                                                                                                        aaTTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTCTCCCATACAGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGTCGATGACATCGTTAAGATGCTCCATGATCCAGCACTTAACCGGCACAATAtt                                                                           >  2:561898/1‑276 (MQ=255)
                                                                                                                                                                                             tAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGCCATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAAAATTCAgg                                                                       >  1:27695/1‑275 (MQ=255)
                                                                                                                                                                                             tAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGTGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGt                                                                      >  2:254796/1‑276 (MQ=255)
                                                                                                                                                                                                            ggCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCGTCTGCCATACCGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTCTGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGTCGATGACATCCTTAAGATGCTCAATGATCCAGCCCTTAACCGGCACAATATTCAGTTTGCCGCTGACGCAAg                                                       >  2:661478/1‑276 (MQ=255)
                                                                                                                                                                                                                gCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACACTGGTAAGATGCTCAATGATCCAGCACTTAACAGACACAATATTCAGGTTGCCGATGACGCAAGAttg                                                   >  4:56911/1‑275 (MQ=255)
                                                                                                                                                                                                                 ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTTGATTACCTTATCACTTATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGtttt                                                   >  1:430855/1‑275 (MQ=255)
                                                                                                                                                                                                                 ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCTTAAAAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTTCCGATGACGCATGACTTg                                                  >  2:573834/1‑276 (MQ=255)
                                                                                                                                                                                                                              aCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGCTAAGATGCTAAAGGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTTGTTTAAATgggg                                      >  3:291835/1‑275 (MQ=255)
                                                                                                                                                                                                                              aCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTCGGCTATCGCAGGCAAACGATAGTAATGGAACCAAAGGCTATGTTTTTGTAGGGGAATCGATGGATTATCTTATCACTTATGGTGCCGATTACATCGTTAAGCTTCTCAATGCTCCAGCACTTAATCGGCCCAATATTCATGTTGCCTCTTGCCCAAGATTTTTTTTTAATACgtg                                     >  4:298416/1‑274 (MQ=255)
                                                                                                                                                                                                                              aCTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTCTCCCTTACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATGACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGAACAATATTCAGGTTGCCGATGACTCAAGATTTGTTTTAAATGCgg                                      >  1:78641/1‑275 (MQ=255)
                                                                                                                                                                                                                                    gcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGCGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCTAAAGATCGTATTGGAACAAAAGGCTATGTTTTTGTAGGGGAATTGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGAGTCTCAATGAACCAGCACTTAACCGGCACCATATTCAGGATGCCGATGCCGCATGATGTGGTTTACATGCGTGGacaaa                               >  2:154176/1‑276 (MQ=255)
                                                                                                                                                                                                                                        aTGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGa                                                                                                                            >  1:441123/1‑179 (MQ=255)
                                                                                                                                                                                                                                            tttCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAc                                                                                                                        <  3:441123/179‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAAAAATTTACCgg                   >  1:237302/1‑274 (MQ=255)
                                                                                                                                                                                                                                                     gTATCTCGTATGCCGTCGTCTGCCTTACAGAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTTCCGATGACATCGTTAAGATGCTCAATGATCAAGAACTTAACCGGCACACTATACAGGTTGACGCGCACGCAAGAGTTGTTTTAAATGAGGGGAACAAGAAATTTACCGACACaa              >  2:573081/1‑276 (MQ=255)
                                                                                                                                                                                                                                                           cGTATGATATGGTTGCCCATACGGCAAAGTGTAGTTCAGGAAAGTAGCAATCGCAGTCCGACACCATCACTGGGCTATCGCAGTAAAAAGATAGTAATGGAACAAAAGGCTATGTGTTTGTAGGGGAATCGCTGGATGACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGc        <  4:254796/276‑1 (MQ=255)
                                                                                                                                                                                                                                                            gTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATTACATCGGTAAGATGCTAAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGGTTTGTTTTAAATGCGGGGAAAAAAAAATTTAACGGAACAATACCGc        >  3:551174/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                gCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCt                                                                                                                                                       >  3:653099/1‑128 (MQ=255)
                                                                                                                                                                                                                                                                 ctctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtac   >  3:536447/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  tctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAAGGCGGGGAAAAAGAAATTTACCGGCAAAATATCGCTTtact  >  4:16464/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  
CGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACT  >  NC_000913/3875263‑3875795

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