Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,6680CT26.7% 122.6 / 13.9 60A54T (GCA→ACA) yiaYL‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (36/8);  new base T (0/16);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 4.92e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAA  >  NC_000913/3757467‑3757887
                                                                                                                                                                                                           |                                                                                                                                                                                                                           
cGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTaa                                                                                                                                                                                                  >  1:1436732/1‑229 (MQ=255)
       gagtatagccccgccATTGCGTGCCACCAGCGAAATACCGTTTGCCCAGTCGTGTGGAGACCCACCGCCTAAGGAGATCACGCTATCGAAATTATTCTCTTTAATTAATTACAACCCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCAAAAAAACCGCTAAAATTATTGCGTTCTTCCAGTTTTTTTTGCACATCGCCCGCCATACCTAATTTCGtt                                                                                                                                                                                           <  2:1139908/219‑1 (MQ=255)
                     gatgggcggcCACCGGCGAAATCCCTTTTGCGCAGTCGTGGGGAGAACCACCGCCTAGGGAGATCACGCTATCGCAATTATTCTTTTAAAGTAATTTCAACCCTGCGGCGACTTTTTCCGTGGTGGGTTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCTTTCTTCCATTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGt                                                                                                                                                                             <  2:566508/225‑1 (MQ=255)
                       ttGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGa                                                                                                                                                                           >  1:1873735/1‑229 (MQ=255)
                                          tACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGt                                                                                                                                                           >  1:605719/1‑226 (MQ=255)
                                               tttGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCa                                                                                                                                                    >  1:429438/1‑228 (MQ=255)
                                                  gcgatgTCGTGGGGACACCCCCCGCCTAGGGAGAGCGCGCTATCGCAATTTTTCTTTTTAAGTAATTTCAACCCTGGGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATAAAAAATAACGCAAAAAATTTTTCTTTCTTCCAGTTTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAAt                                                                                                                                                <  2:428672/224‑1 (MQ=255)
                                                      aGTCGTGTGGAGAACCACCGCCTAAGGCGATCAAGCTACCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGCAAATCCATAATCTg                                                                                                                                             >  1:423558/1‑228 (MQ=255)
                                                        tCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCca                                                                                                                                          >  1:1256315/1‑229 (MQ=255)
                                                                  accccccGCCTAAGGAGATCCCGCTCTCGCAATTATTCTCTTTAAGTAATTTCAACCCTGCGTCGCCGTTTTGCGGGGTGGGTTTGGGTTGGTTGCCATCATAAAAACCCCAAAAAATTTTGGTTTCTCCCGGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCg                                                                                                                                                                                             <  2:313980/163‑1 (MQ=255)
                                                                  aaCCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcgcccgcc                                                                                                                                                                                                          >  1:1950590/1‑155 (MQ=255)
                                                                  aaCCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCg                                                                                                                                                                                             >  1:314548/1‑168 (MQ=255)
                                                                  aaCCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTaa                                                                                                                                                                                         >  1:1103754/1‑172 (MQ=255)
                                                                  aaCCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGAAAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCAt                                                                                                                                >  1:1776454/1‑229 (MQ=255)
                                                                       ctgcgaaaGGAGACACCGTTATCGAATTTTTTCTTTTAATTTAATTTCAACCCTGCGGGGTCGTTTTCCGTGGGGGGGTTGGGTGGGTTCCAAAAAAAAAAAACGAAAAAAATTTTGCGTTCTTCCAGTGTTTTTGGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCAt                                                                                                                           <  2:604658/223‑1 (MQ=255)
                                                                          cagagcgcgaTCTCGTTGTCGCATTTTTTCTTTTAAATTATTTCAAAGCCGTGGGATCCTTTTCCGGTGTTGGTGTTAGGTTGGGTTCCTTAAAAAATACCGAAAAAAATTTTTCTTTTTCCCTGTTTTTTTGGCCCTCCACCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCa                                                                                                                         <  2:840123/220‑1 (MQ=255)
                                                                                  gctCACGCTCTCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCCTACCTAATTTCGTTAACATATTGTCAGTGACACTTACGGTACGGGTAAATCCATAATCTGCAATCATATTCATTGCATCAGTCAATgca                                                                                                                >  2:1190268/3‑227 (MQ=255)
                                                                                             cGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCttttt                                                                                                                                                                                                                        >  1:1706896/1‑114 (MQ=255)
                                                                                             cGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCttttt                                                                                                                                                                                                                        <  2:1703984/114‑1 (MQ=255)
                                                                                               cAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGtt                                                                                                                                                                                           >  1:910803/1‑141 (MQ=255)
                                                                                                    tttCTTTTGAAGTAATTACACACGTGCGGAGGCTTTTTCCGTGGGGGGGTTAGGTTGGGTGCCATAATAAATAACACTAAAAATATTGCGTTCTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTc                                                                                              <  2:1955377/228‑1 (MQ=255)
                                                                                                      tctcTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTcac                                                                                            >  1:1201557/1‑229 (MQ=255)
                                                                                                           ttAAGTAATTTAAAACCTGGGGCGACTTTTTCCGTGGTGGGTTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCCCATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAg                                                                                       <  1:1749044/229‑1 (MQ=255)
                                                                                                                   tttCAACCCTGCGGCGATGTTTGCGGTGGTGGGTTGGGGTGGGGTCCAATCAAAAATACCCAAAAAATTTTTTGTTTCTTCCATTGTTTTTGGCCCATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAg                                                                               <  2:1254133/229‑1 (MQ=255)
                                                                                                                    ttCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCg                                                                                                                                                                                             >  1:809802/1‑118 (MQ=255)
                                                                                                                    ttCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATaaa                                                                                >  1:1397505/1‑227 (MQ=255)
                                                                                                                    ttCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATaaa                                                                                >  1:623094/1‑227 (MQ=255)
                                                                                                                                   ttgTTTTTCGTGGTGGGGTTGGGTGGGTTACCAAAATAAAAACCGATAAAAATTTTTCTTCCTACCATTGTTTTTTCAACAGCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCa                                                               <  2:1680748/227‑1 (MQ=255)
                                                                                                                                      agatggggtggtggtgTTTGTGGGGTTGCAAACATAAATACCGCAAAAAATGTTGTTTTTTTCATGTTTTTTTTGCACATGCACCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGc                                                                    <  2:448812/215‑1 (MQ=255)
                                                                                                                                          tcgttggtggGGTTAGGTTGGGTGCCTTCATAAATACCCATAAAAATATTGTTTTTTTCCGTTGTTTTTGCCACATCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacac                                                         <  2:622001/225‑1 (MQ=255)
                                                                                                                                                 tggGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGa                                                                                         >  1:1693483/1‑190 (MQ=255)
                                                                                                                                                       tAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGc                                                                    >  1:449610/1‑205 (MQ=255)
                                                                                                                                                             gggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCt                                                                                                                                                                                                                            >  2:1149951/1‑47 (MQ=255)
                                                                                                                                                             gggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCt                                                                                                                                                                                                                            <  1:1151962/47‑1 (MQ=255)
                                                                                                                                                              ggTGGCATCCTAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATTGCCAGACATACCTAATTTCGTTAACATATTTTCAGTGACAACTACGGAACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCCCACAAGGAATAAAGAACGTTGAAGCAGCCATAACACACTCCTGAAAGTGTaaaaa                                      >  2:1180009/1‑228 (MQ=255)
                                                                                                                                                              ggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGAc                                                                                                                                                                          >  2:358968/1‑96 (MQ=255)
                                                                                                                                                              ggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCcatc                                                                                                                                        <  1:1162554/130‑1 (MQ=255)
                                                                                                                                                              ggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCcatc                                                                                                                                        >  2:1160528/1‑130 (MQ=255)
                                                                                                                                                              ggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGaaaaaa                                      >  1:410147/1‑228 (MQ=255)
                                                                                                                                                              ggTCCCATCAAAAAAAACCCTAAAAATTTTGTTTTCTTCCGGTTCTTTTTACCCATCTCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGAc                                                                                                                                                                          <  1:359625/96‑1 (MQ=255)
                                                                                                                                                               gTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTTAAGCTGCCATAACACACTCCTTAAAAGGTAAAAAAg                                    >  2:924358/1‑229 (MQ=255)
                                                                                                                                                                      aaaaaataaCGAAAAAAATATTGCGTCTTTCCAGTTCTTTTTGCAAACCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAgggcggg                              <  2:1395104/225‑1 (MQ=255)
                                                                                                                                                                             aaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTa                       >  1:1341303/1‑228 (MQ=255)
                                                                                                                                                                               cacaaaaaaTATTGCTTTCTTCAATTTTTTTTTCCCCATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGaa                                                                                        <  2:894737/157‑1 (MQ=255)
                                                                                                                                                                               cacTAAAAATTTTGCGTTCTCCCATTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGt                                                                                                                                                                 <  2:529274/86‑1 (MQ=255)
                                                                                                                                                                               cGCTAAAATTATTCTGTTTTTCTAGTGTTTTTTACAAATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCaa                                                                                                                    <  2:849530/133‑1 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAAt                                                                                                                                                                                                 >  1:1177573/1‑56 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTTTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGAAGCTTACTg                    >  1:1785207/1‑229 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGt                                                                                                                                                                 >  1:530214/1‑88 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCaa                                                                                                                    >  1:851015/1‑133 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGaa                                                                                        >  1:896300/1‑161 (MQ=255)
                                                                                                                                                                               cGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAg                                                                               >  1:872705/1‑170 (MQ=255)
                                                                                                                                                                                         atTGCTTTCTTCAATTTTTTTTTGCCCATCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCa          <  1:1232195/229‑1 (MQ=255)
                                                                                                                                                                                           tGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCTTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAgcgc                                                                                                         >  1:918729/1‑132 (MQ=255)
                                                                                                                                                                                           tGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg                                                                                                                                             <  1:1508614/96‑1 (MQ=255)
                                                                                                                                                                                           tGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTg                                                                                                                                             >  2:1506005/1‑96 (MQ=255)
                                                                                                                                                                                            ggtgtatcaccAGTTTTTTTGACACTTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCAtt        <  2:1758491/220‑1 (MQ=255)
                                                                                                                                                                                            gCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGACGGACAGCTTACTGAGGATTTTCAtt        >  1:1761486/1‑228 (MQ=255)
                                                                                                                                                                                                  tttcCAGTGTTTTTTGCACACCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGAAAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGaaaa  <  1:1965344/228‑1 (MQ=255)
                                                                                                                                                                                                        gTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGc                                                                    >  1:1221322/1‑157 (MQ=255)
                                                                                                                                                                                                           cTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGaaaa  >  2:1961965/1‑220 (MQ=255)
                                                                                                                                                                                                           cTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGaaa   <  2:1821463/219‑1 (MQ=255)
                                                                                                                                                                                                           cTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGaaa   >  1:1824566/1‑219 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                           
CGTAATCGCGAATATCGCCGCCATTGGCTGCCACCAGCGCAATACCTTTTGCGCAGTCGTGTGGAGAACCACCGCCTAAGGAGATCACGCTATCGCAATTATTCTCTTTAAGTAATTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAA  >  NC_000913/3757467‑3757887

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