Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,258,9960TG30.8% 131.0 / 20.0 78V69G (GTG→GGG) ychHDUF2583 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (11/43);  new base G (24/0);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 5.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGC  >  NC_000913/1258767‑1259217
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
cccAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGGCCGGGTGTGGGGGGGTTGAAAAGGTGTCGGAcc                                                                                                                                                                                                                        >  2:200523/1‑239 (MQ=255)
cccAGAAAGTTCCGGAGGAGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGCGGTAATGGTGGTCGGCGTGGGGGATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTCCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCGGCTGGAGGGCATTGACAGTTGTGCg                                                                                                                                                                                                                           >  1:937584/1‑235 (MQ=255)
 ccAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGt                                                                                                                                                                                                                      <  1:837523/239‑1 (MQ=255)
    gAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGGGTGTACAGGTGTGGGGCCGttt                                                                                                                                                                                                                    >  2:1143041/1‑237 (MQ=255)
           ccGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTggg                                                                                                                                                                                                            <  1:200523/239‑1 (MQ=255)
                 gTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGCGTTTGCGGGGATGGACAAGGGGGGTGCCGTTACTgggggggttgc                                                                                                                                                                                                       >  2:1285897/1‑233 (MQ=255)
                                         ttCGTTACTCGGTAACGTGCTCAGGGGGTTGGGTCTGGTGGTAATGGTGGTCGGGGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCGGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGGGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGt                                                                                                                                                                               <  1:247802/238‑1 (MQ=255)
                                                     tAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGTCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:1401735/237‑1 (MQ=255)
                                                        cGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACa                                                                                                                                                                                             >  1:331728/1‑209 (MQ=255)
                                                        cGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACa                                                                                                                                                                                             <  2:331728/209‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGGTGAACAGGTGGGCGGCCGTTACTGGGGGGGTCGCCACCATGGCAAACGCTGCCGCCGGCGCGGTAAACGCCGTCATAGCt                                                                                                                                                         >  2:1401735/1‑238 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  2:905964/239‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1091921/239‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1047396/239‑1 (MQ=255)
                                                                   ggTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCTCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1058974/235‑1 (MQ=255)
                                                                    gTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGCTTAACAGGGGTGCGACCGTTACTGGTGGGGTCGCCCCTATGACCAACGGTGCCGGCGTAGCGGTAATTGCCGTCCTAGCTaaaaaa                                                                                                                                                   >  1:752064/1‑239 (MQ=255)
                                                                           tggggggtatgggtggtCGGCGTGGGTTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTTTTTCGCCCATGGTGCAGTGCTAAGTATTTTCGTCGGGGCCATTCTCTGGCTGGCGGGTGCCCGGGTTGGCGGGCATGAACAGGTGGGCGACCGTTACGGGTGGGTTCGCCACTATGAAAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCa                                                                                                                                            <  2:146074/228‑1 (MQ=255)
                                                                                         ggtCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcag                                                                                                                               <  1:1285897/238‑1 (MQ=255)
                                                                                                        ttCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  1:715942/237‑1 (MQ=255)
                                                                                                         tCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTg                                                                                                              <  1:694077/239‑1 (MQ=255)
                                                                                                           aaTCCTCACCCAGTTACCCAAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  2:1407980/239‑1 (MQ=255)
                                                                                                           aaTCCTCAACCAGTTACCACATTTTAATATGCCCCAGTTTTTCGCACTTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  2:752064/239‑1 (MQ=255)
                                                                                                             tCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGGGTGCGACCGTTACTGGGGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAAAAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCa                                                                                                          >  1:567204/1‑239 (MQ=255)
                                                                                                                  aaCCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGATGAAACGGTGTGCGACCGTTACTGGGGGGTTCGCCCCTATGGCAAACGTTGCCGCCGGAGCGATAATCGCCGGCATAGCTAACAAATAATGCAGGTTTCTGAACAGTCAGTTTGCGCCGACCGGGCaa                                                                                                         >  2:715942/1‑235 (MQ=255)
                                                                                                                                 tttAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGa                                                                                                                                                                                                                          <  1:650671/107‑1 (MQ=255)
                                                                                                                                 tttAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAaata                                                                                                                                                 <  1:392734/180‑1 (MQ=255)
                                                                                                                                     atttGCCCCAGTATTTCGCACATGGTGCGGTGCTAAGTTTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCAt                                                                                  <  1:1069184/236‑1 (MQ=255)
                                                                                                                                           cccAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGCATGAAAAGGGGTGCGACCGTTACTGGGGGGGTCGCCACTATGACAAACGGTGCCGCCCTAGCGATAATCGCCGGCATAGCTTAAAAAAAATGCAGGTTg                                                                                                                                       >  2:392734/1‑180 (MQ=255)
                                                                                                                                           cccAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCTGTTGGCGGGCATGAAAAGGGGTGCGACCCGTACTGGTGGGGTCGCCCCCATGACAAAAGTTGCCGCCGGCGCGCTAATCGCCGTCATAGCTAACAAAAAATTCCGGTTTCTTACCAGTCAGGTTGCGCTGACTGGGCAATCCCGTACTTTTAAAATCGcccaccccc                                                                             >  2:1380229/1‑231 (MQ=255)
                                                                                                                                                tATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgc                                                                                                                                                                                     <  1:1129817/129‑1 (MQ=255)
                                                                                                                                                  tttCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCt                                                                                                                                     <  1:639364/175‑1 (MQ=255)
                                                                                                                                                  tttCGCACATGGGGCAGTGCTAAGTATTTTCGTCGGGGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGACCGGGTGTGCGACCGTTATTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  1:1021353/200‑1 (MQ=255)
                                                                                                                                                       cacaTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGGTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGGACTTATAAATCCGCCATCGCCGCGCGACGAGTCgg                                                               >  1:494652/1‑239 (MQ=255)
                                                                                                                                                          aTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGGGTTGGCGGGCATGAAAAGGGGGGCGACCGGTACTGGGGGGGTTGCCACTATGACAAAAGGTTGCGCCGGAGCg                                                                                                                                                                           >  2:1129817/1‑129 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGGGGGCATGAAAAGGGGTGCGGCCGTTACTGGTGGGGTTGCCCCCTTGAAAAACGTTTCCGCCGGAGCGGAAATATCCGGCCCTGCTAACAAAAAATGCAGGTTTCTGACCAGTCCGTTTTTGCCGACTTGGCCATCTCGTa                                                                                                  >  2:1021353/1‑200 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTTGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcagtcagt                                                                                                                          >  1:1019033/1‑176 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcagtcagt                                                                                                                          <  2:1019033/176‑1 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGCATGAACAGGGGTGCGACCGGTTTTGGGGGGGTCGCCACCATGGAAAACGTTGCCGCCGGAGCGATAAATGCCGTCAGAACTAACAAAAAATGAAGTTTTCTGACcagtcag                                                                                                                           >  2:639364/1‑175 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:389755/132‑1 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:922504/132‑1 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTgcg                                                                                                                     <  1:427423/179‑1 (MQ=255)
                                                                                                                                                                 aGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:701408/141‑1 (MQ=255)
                                                                                                                                                                   tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGc                                                                                                                                                                                                                            <  1:904891/71‑1 (MQ=255)
                                                                                                                                                                   tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAAc                                                                                                                                                      <  1:1294736/141‑1 (MQ=255)
                                                                                                                                                                    gCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGCGTTGGCGGGGCTGGACCGGGGGGCGGCCGGTACTGGGGGGGTTGCCCCTCCGGCAAACGGTGCCGCCGGGGCGGTAATCGCCGTCCCAGCTAACAAAAAATGCGGTTTGCTGACCAGTAAGGTTGCGCTGACCGGTCAATTCCGGACTTTTAAATCCGCCCTCGCCCCGCGCCGACTTCgggcgaacccccgc                                                   >  2:409525/1‑228 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  1:756880/181‑1 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAgg                                                  <  1:783569/238‑1 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGGGGGTTCGCCACTATGACAAACCTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTTCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAAGc                                                 >  1:529639/1‑239 (MQ=255)
                                                                                                                                                                       aaGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGCTGAACAGGGGGGCGGCCGTTACTGGTGGGTTCGCCCCCATGACAAACCTTGCCGCCGGAGGGATAAACGCCGGCCTAGCTAACAAAAAATGCAGTTTTCTGGCCCATCAGTTTTCGCTGACTGGTCAATCTCGTAATTTTAAATCCCCCATCCCCGCGTGCGCGTTTGGGGAGAAAGCCAGGCg                                                >  2:262825/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGGCCGTGTTGGCGGGCCTGAACCGGGGGGGGGCCGGTACTGGTGGGGTCGCCACTATGACAAACGTTGCCGCCGGAGCGATAATCGCCGGCATAGCTAAAAAAAAAAGCAGTTTTCTGACCAATCCGGTTTGGCTGACTGGTc                                                                                                           >  2:427423/1‑179 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGGGGGGCTGAACAGGGGTGCGACCGTTACTGGGGGGGTCGCCACCATGGCCAACGTTGCCGCCGTAGCGGTAATCGCCGGCCTAGGTAAAAAATAAAGAAGTTTTCTGACCCGTCAGTTTGCGCTGACTGGGCAATCTCGTACTTATAAATCCCCCATCGCCGCCCGGCGCTTCCGGGAGAAAGCCAGGCgg                                               >  2:70210/1‑237 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTc                                              <  1:409525/239‑1 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGGGTGCGGCCGGTACTGGGGGGGTCGCCACTATGACAAACGTTGCCGCCGGAGCGATAATCGCCGTCATAGCTAACAAAAAATGCAGGTTGCTGACCAGGCAGTTTTCGCTGACTGGTCAAACTCGTACTTATAAAACCCCCATCCCCCCCCGACGATTCGGGGAGAAAGcccagcgg                                               >  2:782583/1‑232 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACGGGGTGCGGCCGTTACTGGGGGGTTCGCCCCTATGACAAACGTTGCCGCCGTAGCGATAAACGCCGTCCTAGCTAACAAAAAATGCAGGTTTCTGACCAGGCAGGTTTCGCTGGCTGGTCAAACTCGTACTTTTAAAACCCCCCACTCCGCCCGGCCCATCCGGAAGAAAGCCacgtggc                                              >  2:159691/1‑232 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGGGTTGGCGGGGATGAACAGGGGTGCGGCCGGTACTGGGGGGGTCGCCCCTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGGCATAGCTAACAAATAATTCAGTTTTCTGGCCAGACCGGTTGCGCTGACTGGGCAATCTCGTACTTTTAAATCCGCCCTCCCCCCGCGGCCGTTTGGGAAGAAAGCCCGGGGt                                               >  2:438236/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGCATGAACAGGGGTGCGACCGTTACTGGGGGGGTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCCTCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGt                                               >  1:1038782/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCGGTTGGGGGGCATGAACAGGGGTGCGGCCGTTACTGGTGGGGTCGCCACCATGACAAACGTTGCCGCCGTAGCGATAATCGCCGGCATAGCTCACAAAAAATGCCGGTTGCTGACCAGTCAGTTTTCGCCGACTGGTCAAACTCGTACCTATAAATCCGCCATCGCCGCGCGACGATTCGGGGAGAACGCCCGGGGt                                               >  2:783569/1‑238 (MQ=255)
                                                                                                                                                                          tATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGAc                                                                                                                                                                                                                         <  2:1393811/67‑1 (MQ=255)
                                                                                                                                                                          tATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCa                                                                                                                                            <  1:607317/144‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTccc                                            <  1:782583/238‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTccc                                            <  1:159691/237‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGGTGGCGGGGATGAAAAGGGGTGGGACCGGTACTGGGGGGGTCGCCACTATGACAAAAGGTGCCGCCGGAGCGATAATCGCCGTCATAGCTAACAAAAAATg                                                                                                                                              >  2:701408/1‑141 (MQ=255)
                                                                                                                                                                             tttCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGATGGAAAGGGGTGGGGCCGGTACt                                                                                                                                                                                                                  >  2:904891/1‑71 (MQ=255)
                                                                                                                                                                             tttCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCCTGTTGGCGGGGCTGAAAAGGGGTGGGGCCGGTTCTGGGGGGGTCGCCCCTATGACCAACGTTTCCCCCGTAGCGGTAATCCCCCGCCTAGCCAACaaaaaaatca                                                                                                                                            >  2:1294736/1‑137 (MQ=255)
                                                                                                                                                                               tcgtcgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGGATGAAAAGGGGTGCGACCGGTTCTGGGGGGGTCGCCCCCTTGACAAACGTTGCCGCCGGAGCGATAATCGCCGTCCTAGCTAAAAAATAAAGCAGTTTTCTGACCAGGCAGTTTGCGCTCGCTGGTCAAACTCGTa                                                                                                  >  2:756880/1‑181 (MQ=255)
                                                                                                                                                                                   cgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCgg                                    <  1:438236/238‑1 (MQ=255)
                                                                                                                                                                                    ggTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGGATGAACAGGGGTGCGACCGGTACtggt                                                                                                                                                                                                               >  1:1393811/1‑67 (MQ=255)
                                                                                                                                                                                    ggTGCCATTCTCTGGCTGGCGGGGGCCCGTGTTGGCGGGCATGGAAAGGGGTGCCGCCCGTTCTGGGGGGGTCGCCACTCTGACAAACGTTTCCCCCGTAGGGATAATCGCCGTCCTAGCTAACAAATAATGCAGGTTTCTGAc                                                                                                                                  >  2:607317/1‑144 (MQ=255)
                                                                                                                                                                                     gtgccgtcttctGGCTGGAGGGTGCCGGTGTTGGCGGGCATGAACAGGTGTGCGACCGTGACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCAGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGGCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGtt                                  <  2:1042701/229‑1 (MQ=255)
                                                                                                                                                                                     gtcccATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGtt                                  <  1:1143041/235‑1 (MQ=255)
                                                                                                                                                                                              tcgggctgggGGGTGCCCGTGTTGGCGGGAAGGAAAGGGTGTGGGACGGTTACTGGTGGGTTCGCCAATATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGTGTCCCGGGATCGGTTGAATGCCAgc                        <  2:567204/236‑1 (MQ=255)
                                                                                                                                                                                                      gcGGGTGCCCCTGTTGGCGGGCATGAAAAGGGGTGGGACCGTTACTGGTGGGTTCGCCCCTATGGCAAAAGTTGGCGCCGTAGCGGTAATCGCCGGCATTGCTAACAAATAAAGCAGTTTTCTGACCCGTCAGTTTGCGCTGACTGGTCAAACTCGTACTTTTAAATCCGCCCTCCCCCCGCGGCGGTTCGGGGAGAAAGCCCGGCGTCCCGGGCTCGGTTGAATGCCacacggagcc                 >  2:411221/1‑229 (MQ=255)
                                                                                                                                                                                                               cGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCa       <  1:174621/239‑1 (MQ=255)
                                                                                                                                                                                                                     tGGCGGGCATGAACGGGTGGGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGc  <  2:1038782/239‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCAc                                                                                                                                                                                                    >  2:470219/1‑45 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCAc                                                                                                                                                                                                    <  1:470219/45‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAgg                                                  >  1:300808/1‑190 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAgg                                                  <  2:300808/190‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGTAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  2:274096/113‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 >  1:274096/1‑113 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 >  1:1268744/1‑113 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  2:1268744/113‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAAtctc                                                                                                     >  1:78870/1‑125 (MQ=255)
                                                                                                                                                                                                                                     gtgtGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAAtctc                                                                                                     <  2:78870/125‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
CCCAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGTTGGCGGGCATGAACAGGTGTGCGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGATTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGC  >  NC_000913/1258767‑1259217

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