Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG26.0% 137.6 / 60.5 79A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (23/34);  new base G (6/14);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGT  >  NC_000913_3_bme_pgi/526649‑526942
                                                                                                                                                   |                                                                                                                                                  
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:95959/1‑149 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                 <  2:715268/149‑1 (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)
                                                                                                                                                 cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGt  >  2:402627/1‑149 (MQ=255)
                                                                                                                                                 cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGt  <  1:169050/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGT  >  NC_000913_3_bme_pgi/526649‑526942

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