Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,3690GA70.6% 61.3 / 32.2 56A96T (GCG→ACG) torCcytochrome c menaquinol dehydrogenase TorC
Reads supporting (aligned to +/- strand):  ref base G (10/5);  new base A (16/20);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTG  >  NC_000913/1058123‑1058615
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
cGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGGCACCAGCACAACAGAATTTTGTATCAGTTGCCCCAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTTCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGa                                                                                                                                                                                                                                                      <  1:12093/250‑1 (MQ=255)
     gTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATg                                                                                                                                                                                                                                                >  2:433826/1‑251 (MQ=255)
         gTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGatat                                                                                                                                                                                                                                            >  1:475223/1‑251 (MQ=255)
          tACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATc                                                                                                                                                                                                                                           >  2:420277/1‑251 (MQ=255)
           aCTGGCGCTGGTCGCAATTGGGATTGTGATGGGCATTTAGCTGATTGTAGTGCCACGCGTTGGGATCAAAGTCCCCAGTACAACCGAATTTTGTGTCAGTTGCCACAGTATGCATCCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCGGAATGACATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGTTATCt                                                                                                                                                                                                                                          <  2:353593/251‑1 (MQ=255)
            cTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGCATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCt                                                                                                                                                                                                                                          >  2:571795/1‑250 (MQ=255)
                    ggTCGCAATTGGGATTGTGATTGGCATAGTGCTGATTGTATTGCCACTCGTTGGGATCAAAGTCACCAGCACATCCGAATTTTGTGTCATTTGCCACAGGATGCAACCGGTGTCTCAAGAATATAAACAGACGTTGCATTTCCAGAACGCCTCCTGCGTGCTAGCTTAATGCCATGACTGTCATATCTCTCCGTATATTCCAGGCATGGTGAAGCGCAAACTGTATTCGAGCAATGATATCTATTAGACCt                                                                                                                                                                                                                                 <  2:318565/251‑1 (MQ=255)
                                     tgATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCCCAACAGAATTTTGTGTCAATTGCCACAGTATGCAACCGGTGTATGAAGAATATACTCAGTCGGTGCATTTACAGCACGCCTCCGGCTTGAGAGACGAATGCCATGAATGTCATATACCGCCGGCTATTCCGGGCATTATGAAGCGCACACTGGCAACGAGCAATGATATCTACCACACCTTTATTGATCACTCCAtt                                                                                                                                                                                                                >  1:374580/1‑251 (MQ=255)
                                                gcgcTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTg                                                                                                                                                                                                      <  1:59767/250‑1 (MQ=255)
                                                         gTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCCATGATATCTACCCGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCg                                                                                                                                                                                             >  2:203813/1‑250 (MQ=255)
                                                              gCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGcc                                                                                                                                                                                        <  1:173224/250‑1 (MQ=255)
                                                                     gTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACGGTATGCAACCGGTGTATGGAGAATATAAACAGGCGGTGCATTTCCAGAACGCCTCCGGCTTGCGGGCTGAATGCTATGACTGTCATAGCCGTCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACTAGCAAGGATATCTACCAGACCGTTATTGCTCAGTCCAGTGATACTCCTGAAATATTCGAAGCCAGAcgcgc                                                                                                                                                                                <  2:567410/251‑1 (MQ=255)
                                                                        gggATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGa                                                                                                                                                                             <  1:571795/251‑1 (MQ=255)
                                                                             cAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTg                                                                                                                                                                        <  2:22285/251‑1 (MQ=255)
                                                                                 gTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGa                                                                                                                                                                    <  1:433826/251‑1 (MQ=255)
                                                                                    aCCAGCACAACCGAAGTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACATTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATGCCAGGCATGGTGAAGCGCAAACTGTAAGCGAGCAATGATATCTACCAGACCTGTATTGCTCACTCCATTGATACACCTGAAAAATTCGTAGCCAAACGCGCGTAACTTGCCGAGCg                                                                                                                                                                 <  2:627791/251‑1 (MQ=255)
                                                                                          acaACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTAACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATg                                                                                                                                                          >  2:406317/1‑251 (MQ=255)
                                                                                             aCCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGgc                                                                                                                                                       <  1:13971/251‑1 (MQ=255)
                                                                                                 aaTTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATCTAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACt                                                                                                                                                                                                                     >  1:324759/1‑186 (MQ=255)
                                                                                                 aaTTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCGACCAGACCTTTATTGCTCACt                                                                                                                                                                                                                     <  2:324759/186‑1 (MQ=255)
                                                                                                  aTTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGa                                                                                                                                                   >  2:52993/1‑250 (MQ=255)
                                                                                                    tttGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGCGCGTGAATGGGCGCGAATg                                                                                                                                                >  2:1664/1‑251 (MQ=255)
                                                                                                             gTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATTGGCGCGAATGAAAGAAaac                                                                                                                                       <  2:351336/251‑1 (MQ=255)
                                                                                                                   acaGTATGCAACCGGTGTATGAAGAATATCAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGCTATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGACGCGACTGAAAGAAAACAACTCg                                                                                                                                 >  2:405239/1‑251 (MQ=255)
                                                                                                                          gCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAAcc                                                                                                                           <  1:138123/250‑1 (MQ=255)
                                                                                                                           cAACCGGTGTAAGAAGAATATAAGCAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGTCATATACCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTCTTGCTCCCTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCGGAGCGTGAATGGGCGCGATTGAAAGAAAACAACTCGGCAACCTg                                                                                                                         >  2:90208/1‑251 (MQ=255)
                                                                                                                                  tgtATGAAGAATATAAAAAGGCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGACGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATCCACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGATCGGGAATGGGCGCGAATGAAAGAAAACAACTCGGCACCCTGCCGCTcc                                                                                                                  >  1:4915/1‑251 (MQ=255)
                                                                                                                                   gtATGCAGAATATAAACAGACGGGGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGCCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAACCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCAAAACGCGCGGACCTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGTCAACCTGCCGCCCCt                                                                                                                 >  2:218419/1‑251 (MQ=255)
                                                                                                                                    tATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAATCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTg                                                                                                                >  1:444806/1‑251 (MQ=255)
                                                                                                                                        aaaaaaaaaacacgTCGGTGAATTTCAAGAACCCCGCCGGCTTGAGTGCGGAAGGCCATGACGGCAATTTCCCGCCGGATATTCCGGGCATGGTGAAGCGCAAACGGGAATCGAGCAATGATATCTCCCAGTCCTTTAGTGTTTACGCCATTGATACCCCTGAAAAATTCGAAGCAAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAAAAACCAACTCGCCAACCTGCCGCCCCTGCCAt                                                                                                            <  2:1229/238‑1 (MQ=255)
                                                                                                                                             tataAACAGTCGGTGAATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGCATGGGCGCGAATGAAAGAAAACAACTAGGCAACCTGCCGCTCCGGCAATAACt                                                                                                        >  2:101909/1‑250 (MQ=255)
                                                                                                                                              ataAAAAGTAGGTGCATTTCAAGAACGGCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCACGCCGGATATTCCAGGCATGGTGACGCGCAAACTGGAAACGAGCAATGATATCTACCAGACATTTATTGATCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGAGCGAATGAAAGAAAACAACTCGGCAAACTGCCGCTCCTGCCATAACtt                                                                                                       >  2:612827/1‑249 (MQ=255)
                                                                                                                                              aaaaaCAGGCGGTGCATTTCCAGAACGCCCCCGGCGTGCGAGCTGAATGCCATGACTGTCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAc                                                                                                      <  2:255033/249‑1 (MQ=255)
                                                                                                                                               taAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAAATTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATACCTAcg                                                                                                     >  1:615878/1‑251 (MQ=255)
                                                                                                                                                 aaCAGTCGGTTCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCGGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTCcgat                                                                                                   <  1:401407/251‑1 (MQ=255)
                                                                                                                                                                           tCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATcc                                                                         >  1:250428/1‑251 (MQ=255)
                                                                                                                                                                                      aGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCAc                                                              >  1:63060/1‑251 (MQ=255)
                                                                                                                                                                                        cTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGt                                                            >  2:516700/1‑251 (MQ=255)
                                                                                                                                                                                           aaTGCCATGACTGTCATATCCCTCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAATCGAGCAATGATATCTACCAGACCTTTATGGCTCACTCCATTGATACCCCTGAAAAATTCGAAGCCAAACGCGCGTAACTTGTCGATCTTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAg                                                         <  1:570614/251‑1 (MQ=255)
                                                                                                                                                                                            aTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGa                                                        <  1:388882/251‑1 (MQ=255)
                                                                                                                                                                                            aTGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTGGCCGAGCGTGAAGGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGa                                                        <  1:52993/251‑1 (MQ=255)
                                                                                                                                                                                             agCCATGACTGTCATATCCCGACGGATATTCAATGCCTGGTGCATCGCAAACTGAAAACGAGCACTGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCTAAGCCAAACGCGCGGAATTTGCCGAGCGTGAATGGGCGAGAATGAAAGATAACAACTCGGCAACCTGCCGCCCCAGCCATAAATACGCTGCGATGTGTAATGCGACGCAGCCTCCTGACGCAGCACGTc                                                           >  2:427577/2‑247 (MQ=255)
                                                                                                                                                                                                cATGACTGTCATATCCCGCCGGCTATTCCAGGCATGGTGAACCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATAAACCTGAAAAATTCGAAGCCAAACGAGCGGAACTTGCCGAGCGTGAATGGGAGAGAATGAAAGAAACCAACTCGGCAACCAGCCGCTCCCGCCATAACTACGATGCGATGGATAATGCGAAGCATAATCCTGAAGCAGCACGTCTGATgc                                                     >  1:372435/1‑249 (MQ=255)
                                                                                                                                                                                                cATGACTGTCATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCGGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGaa                                                    <  1:516700/251‑1 (MQ=255)
                                                                                                                                                                                                cATGACTGACATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATTTCTACCAGACCTTTATTGCTCACTCCATGGATACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGGGAATGGGCGCGAATGAAAGAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGaa                                                    <  1:56728/251‑1 (MQ=255)
                                                                                                                                                                                                    aCTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTg                                                >  2:264400/1‑251 (MQ=255)
                                                                                                                                                                                                                cgccgGATATTCCAGGCATGGTGAAGCGCACACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGaa                                        >  1:181345/1‑247 (MQ=255)
                                                                                                                                                                                                                  ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCTTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaa                                  <  1:420277/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggATATTCCAGGCATGGTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACACGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGAATCATGCGAAGCAGCATACTGAAGCAGCACGTCAGATGAAGGAGGCAGCGAAAGATaatc                                >  2:129971/1‑251 (MQ=255)
                                                                                                                                                                                                                          gccAGGCAGGGTGAAGCGCAAACTGGAAACGAGTAATGATATCTACCAGACTTTTATTGCTCACTCCATTTATACACCGGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGGCAGAGGAAGGTGGCAGCGAAAGATAATCAATcc                           <  1:498248/249‑1 (MQ=255)
                                                                                                                                                                                                                          cccAGGCATGGGGAAGCGCAAACTGGAAACGATCAAGGATATCTACCAGACCTTTATTGCTCACTCCATTGATCCACCGGAAAAATTCGAAGCCAACCGCGCGGACCTGGCCGAGCGTGAATGGGCGCAAATTAAAGAAAACACCTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCt                          <  1:445092/250‑1 (MQ=255)
                                                                                                                                                                                                                           ccAGGCAGGGTAAAGGGCAAACTGAAAACGAGCAATAAGGTCTGCCAGACGTTTATGGCTCACTCCTGTCAGACACCAGAAAAATTCAAAGCAAAACGCGCGGACCGTGGGAAGCTGAAATGGGCGCGAATAAAAAAAAACACCCCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGAGCATGCGAAGCAGGATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGAGAATCAATCCt                          <  2:385362/250‑1 (MQ=255)
                                                                                                                                                                                                                                cATGGTGAAGCGCAAACTGGAAGCGAGCACTGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAACTTCGAAGCCACACGCGCGGCACTTGCCGAGCGTGAATGGGCGAGAATGACAGACAACAACTCGGCAACCTGCCGCTCCTCCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCCGAAGCAGCACGTCAGCTGAAGGCGGAAGCGAAAGATAATCAATCCTGAATCg                    >  1:281609/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggTGAAGCGCAAACTGGAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCATCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACt                 <  1:285175/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aaagtGGAAACGAGCATTTATATCTACCAACCCTTTATTAATCACCCCATGGTTACACCTTAAAACTTCGAACCAAAACGCCAGGAACTTGCCGAGCGTGAATGGGGGCAAATGAAAGAAAAAAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAgg       <  2:344373/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggAAACGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTg  >  1:88317/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
CGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTG  >  NC_000913/1058123‑1058615

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