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

GCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526647‑526934
                                                                                                                                                 |                                                                                                                                              
gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATCCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                             >  1:298090/1‑149 (MQ=255)
  ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                           >  1:154123/1‑147 (MQ=255)
  ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                           <  2:467641/149‑1 (MQ=255)
  ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTCGCGCTTATCAGCGAAGACGCCAATACGCCGGGCAGcgcg                                                                                                                                           >  2:414197/1‑149 (MQ=14)
     cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                        <  1:769646/149‑1 (MQ=255)
        aGTGCGGGCGAACGGTGTGGCAACTGGCCATACGGGGGAAGGTGGAATACACGCCGTTGCGAATAGCGCAGTTTTAACTCTTCTTACCCCGGGGACTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                     <  2:22459/149‑1 (MQ=9)
           gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                  <  2:154123/149‑1 (MQ=255)
            cGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgatg                                                                                                                                 <  1:525862/149‑1 (MQ=255)
                cTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATg                                                                                                                             <  2:264222/149‑1 (MQ=255)
                 tGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                            <  1:149270/149‑1 (MQ=255)
                 tGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                            <  1:599485/149‑1 (MQ=255)
                  gACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                           <  2:286078/149‑1 (MQ=11)
                  gACGGTGGCGCAACTGGCCACACAGGTGGAGCCGGAATACACACCGGCCCGAAAATCTCATCTTTCTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCTAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                           <  2:180331/149‑1 (MQ=11)
                  gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                           <  1:328207/149‑1 (MQ=255)
                  gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                           <  1:326457/149‑1 (MQ=255)
                          aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCGGCGAAGACGGCAATACGGCGTGGGGCGGGGGGTATGATGAATGGGGTATCCAg                                                                                                                   >  2:546621/1‑149 (MQ=11)
                             aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                <  1:598507/149‑1 (MQ=255)
                             aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGAAGACGCCAATACGGCGTGGCGCGGGGAGGATGAAGAATTGGGCAACCAGCtg                                                                                                                >  2:496535/1‑148 (MQ=11)
                             aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACC‑GGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                       <  2:200406/125‑1 (MQ=255)
                             aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACC‑GGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                       >  1:200406/1‑125 (MQ=255)
                               cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTaa                                                                                                              <  1:772223/149‑1 (MQ=255)
                                   ccAGACAGGTGGAGCTGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                          <  2:615353/149‑1 (MQ=255)
                                   ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                          <  1:154409/149‑1 (MQ=255)
                                   ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                          <  2:570081/149‑1 (MQ=255)
                                   ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGATTATCAACGGAGACGGCAATACGGCGTGGAGCGCGGAATATGGTGAATGGGGCAACCAGCTTAATggg                                                                                                          >  1:180794/1‑147 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATCCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                         <  2:438991/149‑1 (MQ=11)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                     >  2:61318/1‑121 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                     <  1:61318/121‑1 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                         <  2:715268/149‑1 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                         >  2:95959/1‑149 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCCTTGTCACTTCTACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                         <  1:682689/149‑1 (MQ=255)
                                    cAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTTATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                         <  2:57311/149‑1 (MQ=11)
                                     aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGGAACCGGCTTATTgggga                                                                                                        >  2:213720/1‑149 (MQ=255)
                                     aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                        <  1:662387/149‑1 (MQ=255)
                                       aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCCATACAGCGTGGAGCGCGGAATATGATGCATGGGGCAACCAGCTTAGTgaggcgg                                                                                                      >  1:216618/1‑146 (MQ=17)
                                         aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                    <  2:460963/149‑1 (MQ=255)
                                           gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAacc                                                                                                  >  1:707139/1‑149 (MQ=255)
                                            gggAGCCGGAATACTCACCGGCGCGAAAAGGTCATCGTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGCGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                 <  2:758864/148‑1 (MQ=11)
                                               aGCCGGAATACACACCGGCGCGAAAACCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGTCAACACGGCGTGGAGCGCGGAATATGATGAACGGGGCAACCAGCTTAACGAGGAGAACCCGcat                                                                                              <  2:559950/149‑1 (MQ=14)
                                                gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                >  1:716707/1‑130 (MQ=255)
                                                gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                <  2:716707/130‑1 (MQ=255)
                                                 ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatca                                                                                            <  1:241189/149‑1 (MQ=255)
                                                 ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGAGGCGCGGGGAGTATGATTAATGGGGCAACCAGCTTAATGAGGAGAACCGGcataa                                                                                            >  2:753783/1‑147 (MQ=11)
                                                  cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATCGGGCAACCAGCTTAATGAGGAGAACCCGcataac                                                                                           >  2:72693/1‑146 (MQ=255)
                                                    gAATACACGTCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGGTTTGGAATATGATGAATGGGGCAACTAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                         <  1:229544/149‑1 (MQ=14)
                                                    gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                         <  2:98245/149‑1 (MQ=255)
                                                    gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                         <  1:758568/149‑1 (MQ=255)
                                                    gAATACACACCGGCGCGAAAAGCTCATCTTTAACAAAGATATGGCCGGGGACTTCCGGTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                         <  2:660098/149‑4 (MQ=11)
                                                         cacacCGGCCCGAAAAGCTCACCTTTACCACCGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATc                                                                                    <  1:646820/149‑1 (MQ=255)
                                                          acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                   <  1:129668/149‑1 (MQ=255)
                                                             ccGGCGCGAAAAGCTCATCTTTAGCACGGGGACCGCCGGGGACTGCCGGTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGGATATGATGAATGGGGGAACCTGCGTAATTAGGAGAACGCGCATCATGTGTATCAGcc                                                                                <  2:154575/149‑1 (MQ=17)
                                                              cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTACCAGAGAAGCCGGCGGTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGCGAACCCGCATCATGTGTATTAGcga                                                                               >  2:703840/1‑147 (MQ=17)
                                                                  gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAACCGCATCATGTGTATCAGCCGTAcc                                                                           <  1:587119/149‑1 (MQ=255)
                                                                  ctgAAAACGTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAcc                                                                           <  2:216618/147‑1 (MQ=255)
                                                                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAGCCAGCTTAACGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTc                                                                        >  1:681759/1‑149 (MQ=255)
                                                                     aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTc                                                                        <  1:52841/149‑1 (MQ=255)
                                                                           tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGCCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATAAGCCGTACCGCCTGCCAg                                                                  >  2:369070/1‑149 (MQ=255)
                                                                               cTTTATCACTGCGACCAGCGGGGACTGGCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAACGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca                                                              <  2:205937/149‑1 (MQ=255)
                                                                                tttATCACTGCGACCACCGGGGACTGCCGCTAGCCACCATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                                                             <  1:213720/149‑1 (MQ=17)
                                                                                    tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATGCGGCGAGCAGCGTGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGACCTCGCATCATGTGTACCAGACGTACCGTCTGCACGGccagcaga                                                         >  2:577342/1‑148 (MQ=14)
                                                                                      aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAATATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTATATCAGCGGTCCCGTCTGCCAGGGCCGAAGCat                                                       >  2:109129/1‑149 (MQ=17)
                                                                                           gacgagcGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                  <  2:514918/143‑1 (MQ=255)
                                                                                               accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATGTGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGCCGAGGaa                                              >  2:768072/1‑149 (MQ=255)
                                                                                                ccGGTGACTGCCGCTGGCGCGTATCAGCGAAGACGGCAATACGGGGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                             <  2:330030/149‑1 (MQ=255)
                                                                                                 cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc                                            >  1:330030/1‑149 (MQ=255)
                                                                                                     gagtGCGGCTGGCGGTTTGCAGGGAAGAGGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTGAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                        <  2:121108/146‑1 (MQ=18)
                                                                                                        tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTg                                     >  2:576679/1‑149 (MQ=255)
                                                                                                          cttcTGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATTAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                   <  1:482902/146‑1 (MQ=255)
                                                                                                          ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                   >  1:208802/1‑149 (MQ=255)
                                                                                                            gtcttcgcTGATAAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAAATGATGAATGGGGCAACCAACTTAATGAGGAGAACCCTCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAACATGCTGAGGAATCAGGGCTGTACt                                 >  1:697408/6‑149 (MQ=17)
                                                                                                                 cgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAc                            >  1:4286/1‑149 (MQ=255)
                                                                                                                  gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                           >  1:209964/1‑149 (MQ=255)
                                                                                                                       tCAGCGAAGACGGCAATACGGCGTGGAGCTCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                      >  2:599505/1‑149 (MQ=255)
                                                                                                                            gAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCGGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                 <  1:414197/149‑1 (MQ=255)
                                                                                                                               gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtacta              <  1:478067/149‑1 (MQ=255)
                                                                                                                               gACGGCAAGACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtacta              <  1:478093/149‑1 (MQ=255)
                                                                                                                                aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtactac             >  2:547152/1‑149 (MQ=255)
                                                                                                                                aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtactac             >  2:169050/1‑149 (MQ=255)
                                                                                                                                aCGGCAATACGCGGTGGAGCGCGTAATATGATAAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtactac             <  2:373384/149‑1 (MQ=255)
                                                                                                                                  ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa           >  1:28856/1‑149 (MQ=255)
                                                                                                                                     aaTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATcc        >  1:235157/1‑149 (MQ=255)
                                                                                                                                        aCGGCGTTGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt     <  2:135278/149‑1 (MQ=255)
                                                                                                                                           gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa  >  1:598864/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCA  >  NC_000913_3_bme_pgi/526647‑526934

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