Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,767,6670TG24.8% 398.4 / 161.8 219intergenic (+312/‑43)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (81/83);  new base G (25/29);  total (107/112)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913_3_bme_pgi/2767569‑2767766
                                                                                                  |                                                                                                   
tGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCt                                                                                                     >  2:1626848/1‑99 (MQ=255)
  gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                  <  1:2251249/100‑1 (MQ=255)
     cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                <  2:4828371/99‑1 (MQ=255)
     cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGt                                                                                               <  2:2562635/100‑1 (MQ=255)
      tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                              >  2:3006452/1‑100 (MQ=255)
       gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                             >  1:704980/1‑100 (MQ=255)
        tgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                            <  2:3203380/100‑1 (MQ=255)
          aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                           <  1:4633483/99‑1 (MQ=255)
          aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                          <  2:3700915/100‑1 (MQ=255)
            ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                        <  1:5169653/100‑1 (MQ=255)
            ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                        <  1:230658/100‑1 (MQ=255)
            ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                        <  1:3235859/100‑1 (MQ=255)
            ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                        <  1:3382640/100‑1 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                      >  2:3518455/1‑100 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                       <  1:1780940/99‑1 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                       <  1:4040349/99‑1 (MQ=255)
              cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                       <  1:1396977/99‑1 (MQ=255)
               agaaagaaTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                     >  2:5226780/1‑100 (MQ=255)
               aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                     >  1:2487714/1‑100 (MQ=255)
               aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                     >  2:1776571/1‑100 (MQ=255)
                      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)
                           cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                          <  1:673316/99‑1 (MQ=255)
                           cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtt                                                                          <  1:674741/99‑1 (MQ=255)
                           cAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTATTGTTTTACttt                                                                         >  1:4492234/1‑100 (MQ=255)
                             aTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                       <  1:1015684/100‑1 (MQ=255)
                                tATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCa                                                                      >  1:2203709/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)
                                      cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTATTTGTTTTACTTTTCAGAGGAgtg                                                              >  2:4051135/1‑100 (MQ=255)
                                        ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                            >  1:1056714/1‑100 (MQ=255)
                                        ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                            <  1:4557632/100‑1 (MQ=255)
                                        ataaaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGtt                                                                                              <  2:3406896/63‑1 (MQ=38)
                                          gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                          >  1:903166/1‑100 (MQ=255)
                                           aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                         <  1:1868222/100‑1 (MQ=255)
                                            aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                        <  1:1583816/100‑1 (MQ=255)
                                            aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                         >  2:2351173/1‑99 (MQ=255)
                                            aaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                         >  2:895424/1‑99 (MQ=255)
                                             aaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatga                                                       <  1:2915992/100‑1 (MQ=255)
                                                aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaaa                                                    <  1:4128069/100‑1 (MQ=255)
                                                   ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                 >  1:4267732/1‑100 (MQ=255)
                                                   ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                 >  2:1488070/1‑100 (MQ=255)
                                                   ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtt                                                 >  2:2571328/1‑100 (MQ=255)
                                                    tCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                                >  2:1638436/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)
                                                         tagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                           <  1:3957957/100‑1 (MQ=255)
                                                          agataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGt                                          >  2:2153743/1‑100 (MQ=255)
                                                               tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGAGAGGAGATGaa                                                      >  2:861124/1‑83 (MQ=255)
                                                               tCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTGTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                      <  2:2690187/99‑1 (MQ=255)
                                                                cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTACCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAACTTTTGAAGTTCTg                                      >  2:3150558/1‑98 (MQ=255)
                                                                 aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                    <  1:1179570/99‑1 (MQ=255)
                                                                 aTTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                    <  1:1002012/99‑1 (MQ=255)
                                                                    ctcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTGTGAAGTTCTGCAGcc                                 <  2:2891570/99‑1 (MQ=255)
                                                                      ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAATAATTTTGAAGTTCTGCAGCCAc                               >  2:957831/1‑99 (MQ=255)
                                                                          aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACa                           <  1:1937122/99‑1 (MQ=255)
                                                                          aaCAATCCGATGAACGGCAGAACCGCGAGGTTCGGCCGt                                                                                       <  2:653394/39‑1 (MQ=17)
                                                                            cAATCCGATGAACGGCAGAACCTTGCGGTTCTGCCGTTTTAGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACAACAGaa                        >  1:4241116/1‑100 (MQ=255)
                                                                              aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGTGTGGTGATGAAACATTTTGAAGTTCTGCCGCCACTa                             >  2:1208685/1‑93 (MQ=255)
                                                                                 cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccn                   <  1:3454983/100‑2 (MQ=255)
                                                                                 cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   <  2:1973411/100‑1 (MQ=255)
                                                                                 cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                   >  1:3623029/1‑100 (MQ=255)
                                                                                  gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAGTTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                  <  2:2254978/100‑1 (MQ=255)
                                                                                  caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGAAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCATTACAGAACTCCCt                  >  2:3582926/5‑100 (MQ=37)
                                                                                   aaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCtt                 >  1:4030924/4‑100 (MQ=38)
                                                                                    aaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCttt                >  2:4690788/3‑100 (MQ=39)
                                                                                       aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg             <  1:4745236/100‑1 (MQ=255)
                                                                                       aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg             >  1:1096268/1‑100 (MQ=255)
                                                                                       aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg             >  1:3847215/1‑100 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            <  1:5059159/100‑1 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            >  1:2840137/1‑100 (MQ=255)
                                                                                        cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            <  1:5086265/100‑1 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTTTTTTTCTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            >  1:3725688/1‑100 (MQ=39)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            >  1:1711193/1‑100 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg            >  1:257788/1‑100 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg             >  2:2481759/1‑99 (MQ=255)
                                                                                        cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACCACAGAACTCCCTTTCCg             >  1:3370055/1‑99 (MQ=39)
                                                                                         gtcAGAACCGCGCGGTTCTGCCGTGTTTGGTTTACCGTGCAGAGGAGTGGTGATGAAAGATGTTGAAGTTCGGCAGCCACTACAGGTCTCCGTTTCCgg            <  2:2641455/97‑1 (MQ=255)
                                                                                         ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt           <  1:3737065/100‑1 (MQ=255)
                                                                                          gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc          >  2:2396734/1‑100 (MQ=255)
                                                                                           cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTCTCCGGTc          >  2:5013083/1‑99 (MQ=255)
                                                                                             gAACCTCGCGGTTCTGCCGGTTTTGTTTGACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAGGTTCTGCa                                    <  2:1289001/71‑1 (MQ=255)
                                                                                              aaCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAcc      <  2:544989/100‑1 (MQ=255)
                                                                                                ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    >  1:4700664/1‑100 (MQ=255)
                                                                                                ccTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGCGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    >  1:4232249/1‑100 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    >  1:1027787/1‑99 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    >  2:3777299/1‑99 (MQ=255)
                                                                                                 cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGc    <  2:2064582/99‑1 (MQ=255)
                                                                                                  tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTa  >  1:4969679/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTA  >  NC_000913_3_bme_pgi/2767569‑2767766

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: ‑