Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,767,6700CA24.2% 445.4 / 147.6 216intergenic (+315/‑40)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base C (80/83);  new base A (25/27);  total (105/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 220 total aligned reads displayed.

GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATG  >  NC_000913_3_bme_pgi/2767571‑2767768
                                                                                                   |                                                                                                  
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                                >  2:3006452/1‑100 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                               >  1:704980/1‑100 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                             <  1:4633483/99‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                            <  2:3701441/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                            <  1:2322788/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                            <  1:4161247/100‑1 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGcc                                                                                           <  1:2298722/100‑1 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                          <  1:3235859/100‑1 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                          <  1:3382640/100‑1 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                          <  1:5169653/100‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                        >  2:3518455/1‑100 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                         <  1:1780940/99‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                         <  2:2422756/99‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                         <  1:1396977/99‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                         <  1:4040349/99‑1 (MQ=255)
             aGAATGATTCTGCAATCTATGTGCTTTTAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGAGGAACTGCAGAACCGCGAGGTTCTGCCGt                                                                                         <  2:1177302/98‑1 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                       >  2:1776571/1‑100 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                      >  2:1259086/1‑100 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCACTCTCCTGAACACTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                      >  1:4441074/1‑100 (MQ=255)
              gAAAGAGTCTGCAATCTATGTGGATTGAAAAGATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                        <  2:4898891/98‑1 (MQ=38)
                    ttCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtttt                                                                                <  1:4406013/100‑1 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                              >  1:274380/1‑100 (MQ=255)
                      cTGCAATCTATGTGCATGGAAAGACTCTTCAGGGAGATATCTTACCCCTGAACAATCCGCTGAACGCCGGAACCACCCCGTTCGGCCGTATTCGATTTAc                                                                              >  2:1584713/1‑100 (MQ=255)
                        gCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                            >  2:1384622/1‑100 (MQ=255)
                         cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACttt                                                                           >  1:2896426/1‑100 (MQ=255)
                         cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                            <  1:673796/99‑1 (MQ=255)
                           aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                         <  1:1015684/100‑1 (MQ=255)
                              tATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                        >  1:2204273/1‑98 (MQ=255)
                                tgtgCATTGAAAAATTCTTCAGATAGATATCATTCTGGTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                     <  2:3926538/99‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGTTATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                       <  1:308974/97‑1 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                     <  1:5250645/99‑1 (MQ=255)
                                tgtgCATCGAAGACATCTTCAGATAGATATCATTCTCCTGAACAGTCCGATGCACGGCAGAGCCGCGGGGTTCTGcn                                                                                           <  2:429579/77‑2 (MQ=255)
                                  tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGCACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTGTTGTTTTACTTTTCagag                                                                     >  1:1980269/1‑97 (MQ=255)
                                  tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                  <  1:757009/100‑1 (MQ=255)
                                   gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                                 >  2:1637750/1‑100 (MQ=255)
                                    cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTATTTGTTTTACTTTTCAGAGGAgtg                                                                >  2:4051135/1‑100 (MQ=255)
                                      ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                              >  1:1056714/1‑100 (MQ=255)
                                      ataaaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGtt                                                                                                <  2:3406896/63‑1 (MQ=38)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                            >  1:903166/1‑100 (MQ=255)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                            >  1:3918125/1‑100 (MQ=255)
                                         aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                           <  1:1868222/100‑1 (MQ=255)
                                          aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                          <  1:1583816/100‑1 (MQ=255)
                                          aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                           >  2:895424/1‑99 (MQ=255)
                                          aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                           >  2:2351173/1‑99 (MQ=255)
                                             aaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                           >  1:4542417/1‑96 (MQ=255)
                                              aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaaa                                                      <  1:4128069/100‑1 (MQ=255)
                                                 ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                   >  1:4267732/1‑100 (MQ=255)
                                                 ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                   >  2:1488070/1‑100 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                                  <  1:1291251/99‑1 (MQ=255)
                                                    agttagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGGGATGAAAAATTTTg                                                <  2:347900/97‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                 <  1:1830249/99‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                 <  1:3996721/99‑1 (MQ=255)
                                                    agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                <  1:1409908/100‑1 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                               <  1:2003579/100‑1 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                               <  2:4468026/100‑1 (MQ=255)
                                                           tatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCt                                         >  1:4937278/1‑100 (MQ=255)
                                                             tCATGCTCCTGAACAATCCGTCGAACGGCAGAACCTGGCGGTTGTGCAGTTTTTGTTGTACTTTTCAGAGGAGGGGTGATGAAAAATTTTGAAGTTCt                                         <  2:1554894/98‑1 (MQ=255)
                                                              cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTACCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAACTTTTGAAGTTCTg                                        >  2:3150558/1‑98 (MQ=255)
                                                               aatctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAg                                     <  2:4687722/98‑1 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                      <  1:1179570/99‑1 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                      <  1:1002012/99‑1 (MQ=255)
                                                                  ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCATAGGAGTGGTGATGCAAAATTTTGAAGTTCTGCAGCCa                                  >  2:2487002/1‑100 (MQ=255)
                                                                  ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTGTGAAGTTCTGCAGcc                                   <  2:2891570/99‑1 (MQ=255)
                                                                    ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAATAATTTTGAAGTTCTGCAGCCAc                                 >  2:957831/1‑99 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                               >  2:3232655/1‑99 (MQ=255)
                                                                        aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                             <  1:1937122/99‑1 (MQ=255)
                                                                          cAATCCGATGAACGGCAGAACCTTGCGGTTCTGCCGTTTTAGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACAACAGaa                          >  1:4241116/1‑100 (MQ=255)
                                                                           aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                         <  1:3429718/100‑1 (MQ=255)
                                                                               cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                     <  2:1973411/100‑1 (MQ=255)
                                                                               cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                     >  1:3623029/1‑100 (MQ=255)
                                                                                gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                     <  2:2566560/99‑1 (MQ=255)
                                                                                caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGAAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCATTACAGAACTCCCt                    >  2:3582926/5‑100 (MQ=37)
                                                                                 aaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt                   >  1:4030924/4‑100 (MQ=38)
                                                                                     aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg               <  1:4745236/100‑1 (MQ=255)
                                                                                     aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg               >  1:1096268/1‑100 (MQ=255)
                                                                                     aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg               >  1:3847215/1‑100 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTCTTTTCAGAGGCGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              >  2:690624/1‑100 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              >  1:2840137/1‑100 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              <  1:4977011/100‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              <  1:5086265/100‑1 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              >  1:3725688/1‑100 (MQ=39)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              >  1:257788/1‑100 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              >  1:1711193/1‑100 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACCACAGAACTCCCTTTCCg               >  1:3370055/1‑99 (MQ=39)
                                                                                       gtcAGAACCGCGCGGTTCTGCCGTGTTTGGTTTACCGTGCAGAGGAGTGGTGATGAAAGATGTTGAAGTTCGGCAGCCACTACAGGTCTCCGTTTCCgg              <  2:2641455/97‑1 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTATTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc                >  2:2021939/1‑97 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTTAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt             >  2:5090275/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              <  2:2740958/99‑1 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg              <  2:2213425/99‑1 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt             <  1:3737065/100‑1 (MQ=255)
                                                                                       ggCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt             >  2:2100776/1‑100 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc            >  2:4761701/1‑100 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc            >  2:2396734/1‑100 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc            >  1:2182011/1‑100 (MQ=255)
                                                                                          aGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTa          >  2:642135/1‑100 (MQ=255)
                                                                                            aaCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAcc        <  2:544989/100‑1 (MQ=255)
                                                                                              ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc      >  1:601267/1‑100 (MQ=255)
                                                                                              ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc      >  1:4700664/1‑100 (MQ=255)
                                                                                               cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc      >  2:3777299/1‑99 (MQ=255)
                                                                                               cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc      >  1:1027787/1‑99 (MQ=255)
                                                                                                   |                                                                                                  
GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATG  >  NC_000913_3_bme_pgi/2767571‑2767768

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