Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG26.9% 87.1 / 52.2 71A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (6/12);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_bme_pgi/526651‑526944
                                                                                                                                                 |                                                                                                                                                    
gctagccgcgggcggACGGGTGAGCAACTGCGGACACAGGTGGAGCCGGAATAGACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                                   <  2:408606/142‑1 (MQ=11)
gcgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                                   <  1:600667/149‑1 (MQ=255)
     gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGaggagg                                                                                                                                              >  1:272462/1‑144 (MQ=255)
      tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                             >  1:943933/1‑149 (MQ=255)
      tGCGGGCGGACGGGGGAGGAACTGGCGAGACGGGTGGAGCGGGAATACACACCGGGGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                             <  2:285369/149‑1 (MQ=17)
       gcgggcggACGGTGGAGCAACTGGCCTGACAGGTGGAGCCGTAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGGAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                            <  2:636791/149‑1 (MQ=255)
            gtGAGGGGGGGGCAACTGGCCAGGGAGGTGAGGCCGGAATACAGGCCGGCGGGAAGAGCGTATCTTTATCACGGGGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGGGTGGAGCGTTTAATATGATGAATg                                                                                                                                       <  2:374070/148‑1 (MQ=9)
            cTGAGGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATg                                                                                                                                       <  2:350196/149‑1 (MQ=255)
             tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTCTGATGAATgg                                                                                                                                      >  2:738207/1‑149 (MQ=255)
             tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                                      <  2:6914/149‑1 (MQ=255)
             tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                      <  1:1007147/149‑1 (MQ=255)
             tGACGGGGGAGCAACTGGCCAGACAGGTGGAGGCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCAGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGAGCGAGGAATATGAGGAATgg                                                                                                                                      <  2:497269/149‑1 (MQ=17)
              gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                     <  2:819763/149‑1 (MQ=255)
                   tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACGCCGGCGCGAAAAGCTCATCTTTATCACTGCGATCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATCGGGCAAc                                                                                                                                <  1:930596/149‑1 (MQ=11)
                   tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                <  2:777054/149‑1 (MQ=255)
                         aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCtt                                                                                                                          >  1:756986/1‑149 (MQ=255)
                          aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTa                                                                                                                         >  1:862678/1‑149 (MQ=255)
                           ccggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTAATCAACGAAGCCAGCAATCCGGCCGGGACCGCGGACACCGAAGAATGCGGCCACCAGCtgca                                                                                                                        >  2:204647/3‑146 (MQ=14)
                            tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCCGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACCGCGTGGAGCGCGGAATAAGATGAAAAGGACACCCAGCTTaac                                                                                                                       >  2:279971/1‑148 (MQ=17)
                             ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                      <  1:247643/149‑1 (MQ=255)
                              gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                     <  1:837086/149‑1 (MQ=255)
                              gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                     >  1:448607/1‑149 (MQ=255)
                              gCCAAACAGGAGGAGCCGGAACACCCCCCAGCGCGAAAAGATCATCTTTAACACTGCAACCACCGGGGACTGCCGCTGGCGCGTATCAGCCAAGCCGCCAACACGGCGTGTAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                     <  2:422456/149‑1 (MQ=14)
                               ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                    <  1:434156/149‑1 (MQ=255)
                               ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                    <  1:74071/149‑1 (MQ=255)
                               ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                    <  1:802656/149‑1 (MQ=255)
                                cAGACAGGTGGAGCCGGAATACACTCCGGCGCGAAAAGCTCATCTTTATAAATACGAACACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACTAGCTTAATgagg                                                                                                                   >  1:27025/1‑149 (MQ=17)
                                cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                   <  1:939188/149‑1 (MQ=255)
                                cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                   <  1:39479/149‑1 (MQ=255)
                                cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAAATGATAAATGAGGCAACCAGCTTCAtaaag                                                                                                                   >  1:308203/1‑145 (MQ=17)
                                 aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                  <  2:675542/149‑1 (MQ=255)
                                 aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                  >  1:670890/1‑149 (MQ=255)
                                 aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAACACGGCAATACGGCGTGGCGCGGGGAGTCGGATGAATGGAGCAACAGGCTTAATgcgga                                                                                                                  >  2:663706/1‑149 (MQ=11)
                                 aGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                  <  2:382255/149‑1 (MQ=255)
                                  gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCCGGGACTGCCGCTGGCGCTTACCAGCGAAGACGGCAACACGGCCTGGAGCGAGGAATATGATCAATGAGGCGAGCAGGTGATTaaggag                                                                                                                 >  2:100942/1‑149 (MQ=14)
                                       gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                            >  1:381791/1‑149 (MQ=255)
                                          gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                    >  2:469844/1‑122 (MQ=255)
                                          gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                    <  1:469844/122‑1 (MQ=255)
                                            gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcacg                                                                                                       >  1:504598/1‑147 (MQ=255)
                                            gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGACGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTGAATGAGGAGAACCCGcatc                                                                                                       >  1:956031/1‑149 (MQ=255)
                                               ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                                                                                    <  2:984938/149‑1 (MQ=255)
                                                  aTACACACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGAGCGAGGAATACGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                 <  2:345492/149‑1 (MQ=255)
                                                    acacacCGACGCGAAAAGCCCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCATATCAGCGGATACGGCAATCCGCCGGGGGGGGCGGAATATGAGAAACGGGGAAGACACCCTACGGCGGAGAACCCGCAACACGCGTag                                                                                               >  2:204933/1‑148 (MQ=255)
                                                      acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                             >  1:424768/1‑149 (MQ=255)
                                                            gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                       >  1:958380/1‑149 (MQ=255)
                                                              gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCCAACACGGCAATACGGCGTTGAGCCCGGACTATGATGAATGGGGCACCCAGCTTAGCGAAGCAATCCGGCACCATGAGTACCAGCCGTTcc                                                                                     >  2:668776/1‑149 (MQ=255)
                                                                gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                   >  1:225075/1‑149 (MQ=255)
                                                                      cgcaccTTTATCACTGCCACCACCGGGGACGCCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGGAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                             <  2:691143/144‑1 (MQ=17)
                                                                       tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                            >  1:608845/1‑149 (MQ=255)
                                                                        cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                >  1:278492/1‑96 (MQ=11)
                                                                        cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                                <  2:278492/96‑1 (MQ=11)
                                                                              tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAGTACGGCGTGGCGCGCGGATTATGATGAATGGGGCAACAAGCTTACTGAGGAGCACCCGCGTCAGGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                     >  2:564059/1‑149 (MQ=14)
                                                                                     gCGGCCAACGGGGACTGCCGCTGGCGCTGATCAGCGAAGACGGCAATGCGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGTTTAACGAGGCGACCCCGCCTCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                                              <  2:770655/149‑1 (MQ=17)
                                                                                       gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGACGCAGGAGATGATGAATGGGGCAACCAGCTTGATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCGGCAGCatgatta                                                            >  1:196031/1‑147 (MQ=17)
                                                                                       gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                            <  1:156386/69‑1 (MQ=255)
                                                                                       gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                            >  2:156386/1‑69 (MQ=255)
                                                                                       gACCACCGGGGACTGCCGCTAGGGCTTATCACGGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                            <  2:423000/149‑1 (MQ=255)
                                                                                        cccgccGGGGACTGCCGCGGGGGCTTTTCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                           <  1:450422/145‑1 (MQ=17)
                                                                                        agcaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                           <  1:920044/147‑1 (MQ=255)
                                                                                              ggggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa                                                     <  2:958380/149‑1 (MQ=255)
                                                                                               gcggcTGCCGCTGGTGCGTATCGGCGAAGACGGCAATACGGAGTGGAGCGCGGAACATGATGAATGGGGCAACCCGCTTAAGGAGTAGAGCCCGCATCAGGTGTATCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg                                                    <  2:637863/145‑1 (MQ=17)
                                                                                                 gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                  <  1:677041/149‑1 (MQ=255)
                                                                                                 gACGGCCGCGGGCGCGTAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                                  <  1:222679/149‑1 (MQ=255)
                                                                                                   cTGCCGCTGGCGCTTATCAGCGAATACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGGGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                                <  2:399003/149‑1 (MQ=255)
                                                                                                            gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATaa                                       <  1:150482/149‑1 (MQ=255)
                                                                                                              gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                     >  2:804526/1‑149 (MQ=255)
                                                                                                                               gCAATACGGCGTGGGGCGCGGAATATGATGAAGGGGGCAACCAGCGTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                    <  1:100942/149‑1 (MQ=255)
                                                                                                                                  atccGGCGTGTAGCGCGGAATATGATGAATGGGTCAAGCTGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCg                 <  1:259856/146‑1 (MQ=17)
                                                                                                                                   tACGGCGTGGACCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGt                <  1:638560/149‑1 (MQ=255)
                                                                                                                                            gAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc       <  2:538254/149‑1 (MQ=255)
                                                                                                                                             aGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCg      <  2:448607/149‑1 (MQ=255)
                                                                                                                                              gcgcGGAATATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCgg     <  2:424768/149‑1 (MQ=255)
                                                                                                                                                 cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATATCCTTCACCGGTACTACGATCCGTTGCAGGGGCGGtat  <  1:995298/149‑1 (MQ=255)
                                                                                                                                                 cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat  <  2:670890/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTAT  >  NC_000913_3_bme_pgi/526651‑526944

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