Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,778,3640CT22.4% 398.9 / 112.9 206intergenic (+17/+26)ydiO/ydiPputative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (98/61);  new base T (29/17);  total (128/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

GAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTC  >  NC_000913_3_bme_pgi/1778267‑1778463
                                                                                                 |                                                                                                   
 aaCGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttc                                                                                                  >  2:562740/1‑100 (MQ=255)
 aaCGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGctt                                                                                                   <  2:4565226/99‑1 (MQ=255)
   cGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttctt                                                                                                <  1:73170/100‑1 (MQ=255)
   cGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttctt                                                                                                <  2:285571/100‑1 (MQ=255)
     tATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                              >  2:1272335/1‑100 (MQ=255)
      aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                              >  1:3597263/1‑99 (MQ=255)
      aTCGGCGGCGGTAAAGACGGAATTATGACTTACGTAGCAGTTCGTCCTATCCTGAAAGACCATCAGAACAAATAATCTGCACGCGTCGCAGATTCCTaa                                                                                              >  2:1110249/1‑99 (MQ=255)
       tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                             <  1:553470/99‑1 (MQ=255)
       tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                             >  2:3087368/1‑99 (MQ=255)
       tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                             >  1:2393519/1‑99 (MQ=255)
       tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACa                                                                                            <  2:1190988/100‑1 (MQ=255)
         ggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACaaa                                                                                          >  2:1743146/1‑100 (MQ=255)
           cggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACt                                                                                        >  2:164592/1‑100 (MQ=255)
               ggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgcgc                                                                                    <  1:1371315/100‑1 (MQ=255)
                gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgcc                                                                                   <  1:3325068/100‑1 (MQ=38)
                gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgcc                                                                                   >  1:499806/1‑100 (MQ=255)
                gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgcc                                                                                   >  2:688751/1‑100 (MQ=255)
                 tACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                  <  1:2880579/100‑1 (MQ=38)
                   cAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcc                                                                                <  1:3178287/100‑1 (MQ=38)
                    aGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcct                                                                               <  1:2379064/100‑2 (MQ=38)
                    aGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcc                                                                                <  2:3345590/99‑1 (MQ=37)
                    aGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccgc                                                                                 >  2:3550727/1‑98 (MQ=255)
                     gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctg                                                                              <  1:1278885/100‑3 (MQ=37)
                     gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccgc                                                                                 >  2:3138253/1‑97 (MQ=255)
                     gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccgc                                                                                 >  2:578460/1‑97 (MQ=255)
                     gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAg                                                                              <  1:1764723/100‑1 (MQ=255)
                        gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                           >  1:1146164/1‑100 (MQ=255)
                        gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                           <  2:2115092/100‑1 (MQ=255)
                        gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                           >  1:3531967/1‑100 (MQ=255)
                        gAAATTATGATTTACGTAGCAGGTCGGCAGAGCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                           <  2:2615416/100‑1 (MQ=255)
                              aTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATcc                                                                      <  1:3226756/99‑1 (MQ=255)
                               tggtttACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                    <  2:1948028/97‑1 (MQ=255)
                               tGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                    >  1:1910496/1‑100 (MQ=255)
                                gATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                    >  2:4183384/1‑99 (MQ=255)
                                gATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                    >  1:2957403/1‑99 (MQ=255)
                                 aTTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaaca                                                                  <  2:2037129/100‑1 (MQ=255)
                                 aTTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                   >  1:4132298/1‑99 (MQ=255)
                                  tttACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                   >  1:4380971/1‑98 (MQ=255)
                                   ttACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                   >  2:2798567/1‑97 (MQ=255)
                                        tAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTa                                                           >  2:5051541/1‑100 (MQ=255)
                                           cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTAcc                                                         <  2:177217/99‑1 (MQ=255)
                                           cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTAcc                                                         <  2:1438393/99‑1 (MQ=255)
                                           cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCg                                                        >  1:3796115/1‑100 (MQ=255)
                                            aGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATCTATCCACCAAGACTTACCgg                                                       >  2:914025/1‑100 (MQ=255)
                                               tCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTg                                                    <  1:4254225/100‑1 (MQ=255)
                                                 ggCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAg                                                  >  1:4910420/1‑100 (MQ=255)
                                                  gCAGAGCCTGAAAGACTATCAGAACAAGTAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                   <  2:3886286/82‑1 (MQ=255)
                                                      aTCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAt                                             <  1:2227063/100‑1 (MQ=255)
                                                         cTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAtt                                            >  1:4759004/1‑98 (MQ=255)
                                                              aGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAg                                     >  2:814270/1‑100 (MQ=255)
                                                               gACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAgg                                    <  1:2824845/100‑1 (MQ=255)
                                                               gACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAgg                                    <  1:994859/100‑1 (MQ=255)
                                                               gACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAg                                     >  2:2362272/1‑99 (MQ=255)
                                                               gACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAGTTATCGAACAAGACTTACCGGATGAg                                                  >  2:4709064/1‑86 (MQ=255)
                                                                  tATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAACTTCAGATCGGCtt                                 >  2:4345474/1‑100 (MQ=255)
                                                                  tATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAtt                                 >  2:4061882/1‑100 (MQ=255)
                                                                   aTCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAttt                                >  1:4649395/1‑100 (MQ=255)
                                                                   aTCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAttt                                >  1:1558390/1‑100 (MQ=255)
                                                                   aTCAAAACAAATAATCTGCACGCGGCGCAGCTTCTTAACACACTGCGCAGCCAGATTTACCCAACAAGAc                                                              >  2:3359249/1‑70 (MQ=255)
                                                                    tCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTc                               >  1:3108992/1‑100 (MQ=255)
                                                                    tCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTc                               >  1:2180244/1‑100 (MQ=255)
                                                                    tCAGAAAAAATAATCTGCAGACGGCGCAGCCTCTTAACACACTGCGCCGCCAGATTTATcc                                                                      <  2:1548873/61‑1 (MQ=255)
                                                                      aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                             <  1:2307809/100‑1 (MQ=255)
                                                                      aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                             <  2:2412529/100‑1 (MQ=255)
                                                                      aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                             >  1:3414446/1‑100 (MQ=255)
                                                                      aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCt                              >  2:3772757/1‑99 (MQ=255)
                                                                      aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCt                              <  1:4494587/99‑1 (MQ=255)
                                                                          cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTCCCGGTTGAGGCAATTCAGCTAGGATTTATTGTTc                         >  2:5038705/1‑100 (MQ=255)
                                                                          cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTc                         <  1:3901694/100‑1 (MQ=255)
                                                                          cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTc                         <  2:2173470/100‑1 (MQ=255)
                                                                            aataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCg                       <  2:3737913/100‑1 (MQ=255)
                                                                            aataatCTACAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAACACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTcc                        <  2:4206970/99‑1 (MQ=255)
                                                                             ntaatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCg                       >  1:54234/2‑99 (MQ=255)
                                                                             ataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCg                       >  1:4987641/1‑99 (MQ=255)
                                                                             ataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGc                      <  1:2661906/100‑1 (MQ=255)
                                                                                atCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaa                   <  1:1996642/100‑1 (MQ=255)
                                                                                atCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaa                    >  2:1127224/1‑99 (MQ=255)
                                                                                     cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaaa                  <  1:3499519/96‑1 (MQ=255)
                                                                                     cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                 <  1:5081585/97‑1 (MQ=255)
                                                                                     cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                 <  1:58773/97‑1 (MQ=255)
                                                                                      tggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt             >  1:5014522/2‑100 (MQ=37)
                                                                                      tggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt             >  1:3773238/2‑100 (MQ=37)
                                                                                      aGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt             <  2:4668425/100‑1 (MQ=255)
                                                                                       ggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTGCCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttt            >  2:884903/1‑100 (MQ=37)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaaa                  >  2:3976753/1‑92 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  2:4796919/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  1:1104370/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  2:1353179/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  1:2105884/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  1:4673772/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  2:1585685/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  2:2430853/1‑99 (MQ=38)
                                                                                         cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  1:3164992/1‑99 (MQ=38)
                                                                                         cggcgCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt          <  2:4867509/100‑1 (MQ=255)
                                                                                         cggcgCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  1:2798127/1‑99 (MQ=255)
                                                                                          ggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt          >  2:3931623/1‑99 (MQ=38)
                                                                                            cgcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaaa       >  2:1806363/1‑100 (MQ=255)
                                                                                              cAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaaa       <  1:2778775/98‑1 (MQ=255)
                                                                                              cAGCTTCTTAACAAAATGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGAAAATTCAGATAGGATTTCTTGTTCCGCAAGACAGttt            >  2:5191509/1‑93 (MQ=255)
                                                                                                 cttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc  >  2:2118989/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTC  >  NC_000913_3_bme_pgi/1778267‑1778463

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