Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7920AC20.3% 114.6 / 30.5 71S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (24/31);  new base C (4/10);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_bme_pgi/526648‑526940
                                                                                                                                                |                                                                                                                                                    
cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGGGTGGGgcgg                                                                                                                                                  >  1:681220/1‑148 (MQ=255)
cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGaaggc                                                                                                                                                  >  2:393661/1‑145 (MQ=255)
 ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                 >  2:59619/1‑149 (MQ=255)
    cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACGGGTGGAGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                              <  2:96960/149‑1 (MQ=255)
    cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTAACACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                              <  1:558300/149‑1 (MQ=255)
     gcAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGagt                                                                                                                                             <  1:486343/149‑3 (MQ=255)
      cAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgta                                                                                                                                            <  1:268759/149‑1 (MQ=255)
         ggCGGGCTGACGGTGGCGAAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                         <  1:192529/148‑1 (MQ=11)
         ggCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                         <  1:437009/148‑1 (MQ=255)
          gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                        <  2:318831/149‑1 (MQ=255)
                tggcgttGGAGCAACTGGCCAGACAGGTGGATCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                  <  1:640992/143‑1 (MQ=255)
                tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                  <  1:30614/149‑1 (MQ=255)
                 gACGGTGGCGCAACTGGCCAGACAGGGGGAGCCGGAATACAAACCGGCGCGAAACGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                 <  1:701591/149‑1 (MQ=255)
                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                 <  1:695248/149‑1 (MQ=255)
                  acgtgggAGCAACGGGCCAGACAGGGGGAGCCGGAAGACACACCGGCGCGAAAAGTTCATCTTTAGCACCGGGGCCACCCGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                <  2:23091/144‑1 (MQ=14)
                  aCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                >  2:419066/1‑149 (MQ=255)
                        gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCa                                                                                                                          >  1:124530/1‑149 (MQ=255)
                        gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCa                                                                                                                          >  1:271364/1‑149 (MQ=255)
                              cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa                                                                                                                    <  1:59619/149‑1 (MQ=255)
                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTCTGATGAATGGGGCAACCAGCTTAAt                                                                                                                   >  1:711352/1‑149 (MQ=255)
                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                   <  1:12471/149‑1 (MQ=255)
                                 ggcAGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                 <  2:489364/147‑1 (MQ=255)
                                 gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                           >  2:81930/1‑123 (MQ=255)
                                 gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                           <  1:81930/123‑1 (MQ=255)
                                  cccgACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCCTCTTTATCTCTGCGCCCTCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  2:667482/146‑1 (MQ=17)
                                   cAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  1:124430/149‑1 (MQ=255)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGCATATGATGAATGGGGCAGCCAGCTTAAggagga                                                                                                              >  1:198606/1‑149 (MQ=255)
                                     gACGGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGAAGACGGCAATACGGCGTGGGGCGGGGAATCTGATGAATGGGGCAACCCGCCTAATgaggag                                                                                                             >  2:529160/1‑149 (MQ=14)
                                      aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGgaga                                                                                                            >  1:497372/1‑149 (MQ=255)
                                       cAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                           <  1:187295/149‑1 (MQ=255)
                                          gTGGCGCCGGAATACACACCGGCGCGAAAAGCTCATCCTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                        <  1:385022/149‑1 (MQ=255)
                                            ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                      <  2:497372/149‑1 (MQ=255)
                                             gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                     >  1:497913/1‑149 (MQ=255)
                                              aGCCGGAATACAGACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAGTGAGGAGAACCCGcat                                                                                                    <  2:438525/149‑1 (MQ=11)
                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                   >  1:546455/1‑149 (MQ=255)
                                                 cGGAATACGCACCGGCGAGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCCTAATACGACGAATGGTGCAACCAGCTTAATGAGGAGAACCCGcatcat                                                                                                 <  1:432961/149‑1 (MQ=17)
                                                  ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATTATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcaccacc                                                                                                >  1:156041/1‑147 (MQ=11)
                                                  ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcataacc                                                                                                >  1:156052/1‑145 (MQ=11)
                                                   gAATACACACGGGCGCGAAAAGCTCATTTTTAGCACTGCCACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGCTGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                               <  2:450375/149‑1 (MQ=14)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                               <  1:25381/149‑4 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                               <  2:26284/149‑1 (MQ=255)
                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                               >  1:237849/1‑149 (MQ=255)
                                                     aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                             >  1:449717/1‑149 (MQ=255)
                                                      tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGAGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgca                                                                                            <  2:414532/149‑7 (MQ=11)
                                                       acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                           >  1:711597/1‑149 (MQ=255)
                                                           acCGGCGCGAAAAGCTCATTCTCACCACCGCGACCACCGGGGACTGCCGCTTGCGCTTATCAGCGAACACCTCAAGACAGAAAATATCTCTTAAGATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                       <  1:539875/149‑1 (MQ=12)
                                                                          tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                        <  2:395265/149‑1 (MQ=255)
                                                                            aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTGATGAGGAGAACCCGCATCATGTGTATAAGCCGTACCGTCTGCCAggg                                                                      >  1:530863/1‑149 (MQ=255)
                                                                                     aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                            <  1:61392/134‑1 (MQ=255)
                                                                                     aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                            >  2:61392/1‑134 (MQ=255)
                                                                                          gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                        <  1:419066/149‑1 (MQ=255)
                                                                                          gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAAAACCCGCATCATTTGTAGCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                        <  1:542732/149‑1 (MQ=255)
                                                                                          caccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGAGGATCATCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                        <  2:491814/148‑1 (MQ=255)
                                                                                           accaccGGGGACTGCAGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGTATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTTTATCAGCCGTACTCTCTGCCAGGGCAGCAGCCtgatgaa                                                       >  2:601290/1‑148 (MQ=17)
                                                                                              accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGaa                                                    >  2:559596/1‑149 (MQ=255)
                                                                                                cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATCAGTAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc                                                  <  1:205372/149‑1 (MQ=255)
                                                                                                   ggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg                                               <  1:379349/149‑1 (MQ=255)
                                                                                                    gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                              <  1:217022/149‑1 (MQ=255)
                                                                                                    gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                              <  1:439331/149‑1 (MQ=255)
                                                                                                      cggCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                            <  2:504662/147‑1 (MQ=255)
                                                                                                         ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAACCGCAGCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                         >  2:446324/1‑149 (MQ=255)
                                                                                                             tGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtat                                     >  1:361495/1‑149 (MQ=255)
                                                                                                                 ggttATCGGCGAAGACGGGAATATGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                 <  2:45109/147‑1 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGGATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGTTGAGGAATCAGGGCTGTACTATAAcc                                 >  2:49675/1‑149 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGCATCAGGGCTGTACTATAAcc                                 >  2:137275/1‑149 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCCGGGCTGTACTATAAcc                                 >  2:137151/1‑149 (MQ=255)
                                                                                                                 gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                 >  2:396539/1‑149 (MQ=255)
                                                                                                                    tATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGACAGCAGCATGGTGAGGAATCAGGGCTGTACTATAACCGTc                              >  2:493561/1‑149 (MQ=255)
                                                                                                                         gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGGAACCAGCTTAATGAGGAGAGACCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGACTGAGGGCTGTACTATAACCGTCCCCgg                         >  1:662314/1‑149 (MQ=17)
                                                                                                                              gACGGCAATCCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtacta                    <  2:361773/149‑1 (MQ=255)
                                                                                                                                        cGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCGGTACTATAACCGTCACCGGTACTACGATCCGTTg          <  2:124530/149‑1 (MQ=255)
                                                                                                                                            gTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAgg      >  1:332914/1‑149 (MQ=255)
                                                                                                                                                aGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCg  <  1:559596/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_bme_pgi/526648‑526940

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