Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0230AC32.1% 84.0 / 10.1 54V153G (GTG→GGG) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base C (0/17);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 1.13e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACAT  >  NC_000913/74789‑75251
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
tAGTTCTCTTTCTTCGCTTCGAGAAGGTGATAAGCAGGAGAGGTGGGGTAACTAAGGACCAGTCCGCTTCCACCTTTAAAAAACAGGCCGTAGGCTTCGCCCCATCCTTTGGTGACCGTGCCCGTTTTCTTCGCCCGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGATTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCg                                                                                                                                                                                                                                  <  2:89695/239‑1 (MQ=255)
        tttcttCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGtttt                                                                                                                                                                                                                          >  1:1184168/1‑239 (MQ=255)
          tcttctcttcGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACA‑GCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGa                                                                                                                                                                                                                        >  1:856608/1‑238 (MQ=255)
          tctcctcttcGAAAATGTGATAAGCGGGAGGGGTGGTGTACCTCAGTCCCAGATCTCTTCCCCCTTTAAAAAACA‑GCCGTAGGCTTCGCCCCAGCTTTTGGTGACCGTGACCGTTTTCTTCGCCATTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGCCTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGAt                                                                                                                                                                                                                       <  2:856608/239‑1 (MQ=255)
                        atgtaagaggCCGGAGAGGTGTGGTACCTCGGTCCCGGTTCGCTTACACTTTTGAAAAGCAGGCCGGGGGCTTCGCTCCGGCTTTGGGTGACCGTGACAGTTTCCTTCGCCAGTTTCTGCCAGGCTTGGGGGGGGTCATCGCCATAGACTTTTTGCATCCAAAGCAACAGACCCAGCCCCGGTGTCCTGGTGCGCGGA‑TCCTAAGAAATCACCCGCCAGTTTTGATCGCTCTCAACCAg                                                                                                                                                                                                          <  2:1035081/230‑1 (MQ=255)
                                 gCCGGAGAGGGGGTGTACCCAGGTCCCAGACCGTTTTCACCTTTAAAAAACAGGCCTTAGGCTTCGCTCCAGCTTTGGGGAACCGTGCCCGTTTTCTTCGCCAGTTTCGGCCAGGCTTGTGGGGCGTCATCGCCATAGCCTTTTGGCATCCATAGAAACAGACCCGGCCCCGGTGTCCTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGctc                                                                                                                                                                                                                  <  2:167689/222‑1 (MQ=255)
                                          gtggtgTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGC‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTtgtg                                                                                                                                                                                         >  1:635630/1‑238 (MQ=255)
                                           tggtgTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGtt                                                                                                                                                                                                        >  1:167689/1‑222 (MQ=255)
                                              tgtAGCCCAGCACCAGATCGCTTTCCCTTTTTAAAAACAGGCGGTAGGCTACGCGCCAGCCTTTGGTGCCCGAGCCCCTTTTCTTCGCAAGTTTCTGCCAGGTTTGTGGGGCGTCTTCGCTATAGACTTTTTGCATCCACAACAACAGACCCAGCCCCGGTGTCCTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCaa                                                                                                                                                                                                              <  2:225757/213‑1 (MQ=255)
                                                     cAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGgtttttc                                                                                                                                                                             >  1:1170095/1‑237 (MQ=255)
                                                        tACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGgttttttgtt                                                                                                                                                                          >  1:38696/1‑234 (MQ=255)
                                                        tACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTctcgt                                                                                                                                                                          >  1:949918/1‑235 (MQ=255)
                                                        tACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGATCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGgttttt                                                                                                                                                                              >  2:1025663/1‑234 (MQ=255)
                                                        tACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGGGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCttt                                                                                                                                                                                                    >  1:225757/1‑213 (MQ=255)
                                                         aCCAGGTCGCTTTCACCTTTTAAAAACAGGCCGTGGGCTTCCCTCCAGCCTTTGGTGCCCGTGACCTTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGt                                                                                                                                                                           <  2:949918/237‑1 (MQ=255)
                                                             gATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTGGCCAGTTTCTGCCAGGCTTGTGGGGGGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCTCCCGc                                                                                                                                                                                                                                 <  1:1140315/179‑1 (MQ=255)
                                                             gATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGc                                                                                                                                                                                                                                 >  2:1140315/1‑179 (MQ=255)
                                                             gATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGtt                                                                                                                                                                                                                            >  1:611818/1‑184 (MQ=255)
                                                                  cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCg                                                                                                                                                                                                                                  <  1:158634/173‑1 (MQ=255)
                                                                  cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCg                                                                                                                                                                                                                                  >  2:158634/1‑173 (MQ=255)
                                                                  cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGctc                                                                                                                                                                                                                  <  1:159496/189‑1 (MQ=255)
                                                                  cTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGctc                                                                                                                                                                                                                  >  2:159496/1‑189 (MQ=255)
                                                                                aaccAGCCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGCCTTTTTGCATCCATAGCACCAGCCCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGa                                                                                                                                                  <  1:442499/236‑1 (MQ=255)
                                                                                      gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGc                                                                                                                                                                                                                                 >  2:1206945/1‑154 (MQ=255)
                                                                                      gCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCACCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGc                                                                                                                                                                                                                                 <  1:1206945/154‑1 (MQ=255)
                                                                                             gCTTCCCTCCAGCCTTTGGTGCCCGGGACTGTTTTTTTCGCCAGTTTCGGCCGGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGAt                                                                                                                                                                                                                       <  2:1081273/157‑1 (MQ=255)
                                                                                               tccccTCCGGCTTTGGGAGCCCGTGCCCTTTTTTTTGGCCGTTTTCGGCCGGGCTTAGGGGGCGCAATCGCCATAGTCTTTTGACACCCATAGCACCAGACCCGGCCCCGGTGTCCTGGTGCGGGTA‑TCTTGAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                    <  2:295906/234‑1 (MQ=255)
                                                                                                       aGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAAc                                                                                                                                                                                                             >  1:1081273/1‑157 (MQ=255)
                                                                                                          ttgtGGTAACCTTGCCCGTTTCCTTCGCCAGTTTCTGCCAGGCTTGGGGGGCGCCATCGCCATAGTCTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGTA‑TCCTAATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTAc                                                                                                                            <  2:814334/232‑1 (MQ=255)
                                                                                                                                                     ttGTGGGGCGTCATCGCCATAGACTTTTTGCACCCATAGCACCAGCCCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCg                                                                                                                                                                                                                     <  2:28693/103‑1 (MQ=255)
                                                                                                                                                               tCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGC‑TCCTGAAAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTc                                                                                                                                                                                                   <  2:852273/111‑1 (MQ=255)
                                                                                                                                                               tCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCa                                                                                                                                                                                                           >  1:28693/1‑103 (MQ=255)
                                                                                                                                                               tCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTc                                                                                                                                                                                                   >  1:852273/1‑111 (MQ=255)
                                                                                                                                                                           aCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAg                                                                                                                                                      <  1:1134368/144‑1 (MQ=255)
                                                                                                                                                                                        aaaGCAACGAACCGACCCCCGGTGTACTGGGGC‑CGGACTCCTGTAAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCg                                                                                <  1:369657/199‑1 (MQ=255)
                                                                                                                                                                                        aTAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                    <  1:1084117/149‑1 (MQ=255)
                                                                                                                                                                                            cAACAGACCCAGCCCGGGTGTGCTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCa                                                                                                                                          <  1:1144064/139‑1 (MQ=255)
                                                                                                                                                                                             acaaaccccAGCCCCGTTGTCGTGGTGCGCGGT‑TCCTAATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTCGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTgg                                         <  1:1002545/229‑1 (MQ=255)
                                                                                                                                                                                                  ccccAGCCCCGGTGTCCTGGTGCGCGGG‑TCCTTATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCaaa                                                                                                                                  <  1:673224/140‑1 (MQ=255)
                                                                                                                                                                                                     cggCCCGGGTTTCCTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACa                             <  1:1018453/237‑1 (MQ=255)
                                                                                                                                                                                                      aGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGAt                                                                                                                                                                                                                       >  2:160279/1‑52 (MQ=255)
                                                                                                                                                                                                      aGCCCCGGTGTACTGGTGCGCGGA‑TCCGGATAAATCACCCGCCAGTTTTGAt                                                                                                                                                                                                                       <  1:160279/52‑1 (MQ=255)
                                                                                                                                                                                                        ccccGGTGTCCTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGt                                                                                                                                                                           <  1:615582/94‑1 (MQ=255)
                                                                                                                                                                                                              tttACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTAc                                            <  1:29111/213‑1 (MQ=255)
                                                                                                                                                                                                              tgtCCTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGAGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATc                                                                                  <  1:1094304/177‑1 (MQ=255)
                                                                                                                                                                                                                tACTGGTGCGCGGT‑TCCTAATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                                                    <  2:712725/125‑1 (MQ=255)
                                                                                                                                                                                                                tACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCa                                                                                                                                          <  1:1114031/119‑1 (MQ=255)
                                                                                                                                                                                                                tACTGGTGCGCGGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcg                                       <  1:1410927/218‑1 (MQ=255)
                                                                                                                                                                                                                      tGCGCGTC‑TCCTGAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGccc             <  2:38696/238‑1 (MQ=255)
                                                                                                                                                                                                                         gcgGA‑TCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGc         <  1:610266/239‑1 (MQ=255)
                                                                                                                                                                                                                          cgGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTTTTTTTTTTTTTTTTTAACCAAAGCGGAGGAGCCCTAATCAAAca                                                                                                                                >  2:1114031/1‑118 (MQ=255)
                                                                                                                                                                                                                           ggA‑TCTTTATAACTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGca        <  1:320607/238‑1 (MQ=255)
                                                                                                                                                                                                                            gA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGcac       <  1:1025663/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ttGAAAAACCCCCCGCCAGTTTGGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAGGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCGGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACAt  <  1:814298/238‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
TAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGA‑TCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACAT  >  NC_000913/74789‑75251

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