Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC22.9% 613.2 / 175.2 262S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (86/116);  new base C (27/33);  total (113/149)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-04
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

GCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAA  >  NC_000913_3_bme_pgi/526695‑526890
                                                                                                 |                                                                                                  
gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                  >  1:918489/1‑100 (MQ=255)
  cGGAATACACACCGGCGCGTAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                   <  2:1549331/97‑1 (MQ=255)
  cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcgg                                                                                                <  2:1079033/100‑2 (MQ=255)
   ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                               >  2:1330915/1‑98 (MQ=255)
   ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                               >  1:5131088/1‑98 (MQ=255)
   ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                               >  1:1084668/1‑100 (MQ=255)
    gAATACAGACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                               <  2:1433941/99‑1 (MQ=255)
    gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCgggg                                                                                              <  2:3022409/100‑1 (MQ=255)
    gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                              <  2:1314366/100‑1 (MQ=255)
     attACACACCGGCGCGAAACGCTCATCTTTATCATTGCGACCAGCGGGGACTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCATGGAGCGCGGa                                                                                             <  2:4437553/98‑1 (MQ=17)
     aaTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGTGAAGACGGCAATACGGCGTGGcgcggg                                                                                               >  1:4118779/1‑96 (MQ=255)
      agacacgcCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                            <  2:2396538/98‑1 (MQ=255)
       tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGagt                                                                                           <  1:3058010/100‑3 (MQ=14)
       tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                            >  1:1509114/1‑99 (MQ=255)
         cacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAata                                                                                          <  1:2511510/99‑1 (MQ=255)
          acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatg                                                                                        <  1:1755289/100‑1 (MQ=14)
           cacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                       >  1:3004296/1‑100 (MQ=255)
            acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                      <  1:3227278/100‑1 (MQ=14)
            acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                      <  2:4074931/100‑1 (MQ=255)
             ccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                      <  2:4269307/99‑1 (MQ=11)
              cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATAtgatga                                                                                    <  1:520839/100‑1 (MQ=255)
               ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCNGCTGGCGCTTATCAGCGAAGACGNCAATACGGNGTGGCGCGGNGAGTATGATGaa                                                                                   >  1:427734/1‑100 (MQ=11)
               ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTAGCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                      <  2:3601477/97‑1 (MQ=11)
                gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                  >  1:1349137/1‑100 (MQ=14)
                gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAAt                                                                                  >  1:1867417/1‑100 (MQ=255)
                 cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATg                                                                                 <  1:4204085/100‑1 (MQ=14)
                  gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                >  1:656941/1‑100 (MQ=255)
                  gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                >  1:4995978/1‑100 (MQ=255)
                   cgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                               <  1:3327712/100‑1 (MQ=255)
                    gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                              <  1:4048274/100‑1 (MQ=14)
                    gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                              <  2:2936284/100‑1 (MQ=255)
                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGGAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                             <  2:691519/100‑1 (MQ=11)
                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                              >  1:3136376/1‑99 (MQ=14)
                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGc                                                                             <  1:2357669/100‑1 (MQ=14)
                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGGATGGGGc                                                                             <  2:5245837/100‑1 (MQ=255)
                        aGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                          >  1:104030/1‑100 (MQ=14)
                        aGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                          <  1:3359250/100‑1 (MQ=255)
                         gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCTAATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGAGGAGGAATCGGGCaag                                                                          >  2:3150654/1‑98 (MQ=11)
                          cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                          >  1:4592607/1‑98 (MQ=11)
                          cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                          >  1:4787522/1‑98 (MQ=11)
                             aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCt                                                                     <  1:2141497/100‑1 (MQ=255)
                                tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa                                                                  >  1:1484173/1‑100 (MQ=255)
                                    tCACTGCGGCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                           <  2:2020202/71‑1 (MQ=255)
                                     cACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                              >  2:2867480/1‑99 (MQ=255)
                                      aCTGCGACCACCTGGGACTGCCGCTGGCGCTTATCACCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                             <  1:4422043/99‑1 (MQ=255)
                                      aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                             <  2:4778955/99‑1 (MQ=255)
                                      aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                             <  2:2529058/99‑1 (MQ=255)
                                      aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                             <  2:1364327/99‑1 (MQ=255)
                                      aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                              >  1:1494241/1‑98 (MQ=255)
                                       cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                              >  1:205482/1‑97 (MQ=255)
                                          cGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgaggag                                                           <  1:4496963/97‑1 (MQ=255)
                                          cGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                        <  2:1492817/100‑1 (MQ=255)
                                             ccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                     >  1:1571232/1‑100 (MQ=255)
                                               accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                   <  1:3659477/100‑1 (MQ=255)
                                               accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                   <  2:555396/100‑1 (MQ=255)
                                               accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                   >  2:2093549/1‑100 (MQ=255)
                                               accGGGGACTGCCGCTGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGca                                                   <  2:3864934/100‑1 (MQ=255)
                                                ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                  <  1:1772626/100‑1 (MQ=255)
                                                  ggggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                >  1:4026743/1‑100 (MQ=255)
                                                    ggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                              <  2:935842/100‑1 (MQ=255)
                                                      aCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgtg                                            <  2:784031/100‑1 (MQ=255)
                                                       cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgtg                                            >  1:2688011/1‑99 (MQ=255)
                                                       cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgtg                                            >  2:2729748/1‑99 (MQ=255)
                                                        tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgaggag                                                           >  2:1764346/1‑83 (MQ=255)
                                                         gCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                         >  2:3754151/1‑100 (MQ=255)
                                                         gCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                         >  2:2938723/1‑100 (MQ=255)
                                                          ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                          <  2:4230023/98‑1 (MQ=255)
                                                          ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTa                                          <  1:3569147/98‑1 (MQ=255)
                                                          ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                         >  2:1122315/1‑99 (MQ=255)
                                                            gCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAg                                      <  1:3814936/100‑1 (MQ=255)
                                                                gcgcTTATCAGCGAAGACGGCAGTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCg                                   <  1:5036083/99‑1 (MQ=255)
                                                                 cgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                             >  2:5191767/1‑88 (MQ=255)
                                                                     tATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                             >  2:1272864/1‑100 (MQ=255)
                                                                     tATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                             >  2:4668803/1‑100 (MQ=255)
                                                                     tATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                             >  2:358768/1‑100 (MQ=255)
                                                                      aTCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCg                              <  1:4732605/98‑1 (MQ=255)
                                                                         aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                          <  1:1397632/99‑1 (MQ=255)
                                                                         aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGc                         >  1:2245191/1‑100 (MQ=255)
                                                                             aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                     >  2:2774528/1‑100 (MQ=255)
                                                                             aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                     <  1:2425851/100‑1 (MQ=255)
                                                                             aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                     >  1:58508/1‑100 (MQ=255)
                                                                              aGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                    <  1:1922185/100‑1 (MQ=255)
                                                                               gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                   <  1:310485/100‑1 (MQ=255)
                                                                                 cGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                 >  2:489451/1‑100 (MQ=255)
                                                                                 cGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                 >  1:2756276/1‑100 (MQ=255)
                                                                                 cGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                 <  1:772263/100‑1 (MQ=255)
                                                                                 cGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                 >  1:2563639/1‑100 (MQ=255)
                                                                                   gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca               <  1:3168016/100‑1 (MQ=255)
                                                                                   gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagc                <  2:3279629/99‑1 (MQ=255)
                                                                                    cAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcag              <  2:1931554/100‑1 (MQ=255)
                                                                                      aTACGGCGTGGAGCGCGGAATATGATGAATGGGGGAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagcn            <  1:467773/100‑2 (MQ=255)
                                                                                         cGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga         >  2:2188159/1‑100 (MQ=255)
                                                                                         cGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg          <  1:168433/99‑1 (MQ=255)
                                                                                          ggCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat        <  1:39161/100‑1 (MQ=255)
                                                                                            cGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCat           >  2:453425/1‑95 (MQ=255)
                                                                                             gTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg     <  1:2899698/100‑1 (MQ=255)
                                                                                              tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga      <  2:3364427/98‑1 (MQ=255)
                                                                                              tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg    <  2:892067/100‑1 (MQ=255)
                                                                                              tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg     >  2:2138867/1‑99 (MQ=255)
                                                                                                gAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa  <  1:1243509/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAA  >  NC_000913_3_bme_pgi/526695‑526890

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