Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG21.3% 116.8 / 31.2 62A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base G (5/8);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTA  >  NC_000913_3_bme_pgi/526649‑526943
                                                                                                                                                   |                                                                                                                                                   
ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                    >  2:337438/1‑149 (MQ=255)
    gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACATACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGACGACGGCAATACGGCGTGGCGCGGGGagt                                                                                                                                                >  1:521069/1‑147 (MQ=255)
    gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCAGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCACTACGGCGTGGAGCGCGGAat                                                                                                                                                >  1:621554/1‑149 (MQ=255)
      aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                              <  1:569724/149‑1 (MQ=255)
         gCGGGCTGACGGTGCAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                           <  1:301590/149‑1 (MQ=11)
            ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCCAAGACGGCAATACGGCGAGGCGCGGGGCGTATGATGaa                                                                                                                                        >  2:586907/1‑149 (MQ=11)
               tGACGGTGGAGCAACTGGCCAGACAGGTGGCGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGAGGGGCCAGGAATATGATGAAggg                                                                                                                                     >  1:102614/1‑149 (MQ=14)
               tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                     <  1:108221/149‑1 (MQ=255)
               tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                     <  2:553286/149‑1 (MQ=255)
               tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                     <  1:81751/149‑1 (MQ=255)
                        aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAg                                                                                                                            >  1:734655/1‑149 (MQ=255)
                         gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                           <  2:390727/149‑1 (MQ=255)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCAGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGGGGGAGATGGTGGAAGAGGCAACCAGCta                                                                                                                         >  2:111481/1‑148 (MQ=14)
                           aaCTGGCCAAACCGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCAGCGCTTATCAGCGAAGACGGCAATACGGCGGGGACGACGGAAAATGAAAAAAAGGGCAAcccaatt                                                                                                                         >  2:95324/1‑144 (MQ=14)
                             cTGGCCAGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTCCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAACGCGGAATATAATGAAAGGGGCCACCCGCTTaa                                                                                                                       >  2:227824/1‑149 (MQ=17)
                              tgtccAGACGGGGGGAGCCGGAATACACATCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAt                                                                                                                      <  2:669645/146‑1 (MQ=255)
                              tGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCTCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                      <  2:313590/149‑1 (MQ=11)
                              tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                      <  2:314699/149‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                  >  1:323591/1‑149 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                                  >  1:590562/1‑149 (MQ=255)
                                   aGACCGGTGGAGCCGGGATACACGCCGGCGCGAAAAGCTCATCTTTATTACTGCGACCACCGGGGACTGCCGCTGGCGCTGATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  2:27831/149‑1 (MQ=11)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 >  1:398238/1‑149 (MQ=255)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  1:656995/149‑1 (MQ=255)
                                   aGACAGGTGGAGACGGAATACAGACTGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                 <  2:432977/149‑1 (MQ=11)
                                       aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCTGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                   <  1:476992/127‑1 (MQ=255)
                                       aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                   >  2:476992/1‑127 (MQ=255)
                                            gAGCCGGAATACACACCGGCGCGAAAAGCTCGTCTTTATCACTGCGGCCACCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGGTAATACGGCGTGGGGGGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                        <  2:367374/149‑1 (MQ=14)
                                            gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCAAATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCCACCAGCTTCAAGAGGAGAACCCGca                                                                                                        >  2:482603/1‑149 (MQ=17)
                                             aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAGCCCGcat                                                                                                       >  2:243127/1‑149 (MQ=255)
                                              gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGCGGCGCGCGGAATATCAGGAATGGGGCAACCAGCGTAATGAGGAGAACCCGcatg                                                                                                      >  2:579558/1‑148 (MQ=14)
                                               ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACGCGcatca                                                                                                     >  1:411942/1‑149 (MQ=255)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                  <  1:642440/149‑1 (MQ=255)
                                                  gAACACACACCGGCGGGAAAGGCTCAGCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                                  <  2:408171/149‑4 (MQ=11)
                                                   aaTACCCACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgtg                                                                                                 <  1:759047/149‑1 (MQ=255)
                                                    aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgc                                                                                                <  1:119039/149‑6 (MQ=11)
                                                    aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGAACACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGGATCATgggt                                                                                                >  1:277037/1‑149 (MQ=255)
                                                        acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                            >  2:314047/1‑149 (MQ=255)
                                                             ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGt                                                                                       <  2:506990/149‑1 (MQ=255)
                                                              gggcgAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCTGGCGTTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                      <  2:41530/147‑1 (MQ=255)
                                                                   aaaaTCTCGTCTTTATTACAGCCACCACCGGAGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTc                                                                                 <  2:238264/149‑1 (MQ=18)
                                                                    aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                                >  1:70672/1‑149 (MQ=255)
                                                                     aaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                               <  1:525167/149‑1 (MQ=255)
                                                                      aGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAGGCCGGCAATACGGCGGGGTGCGCGGCATATGATGAATGGCGCAACCAGCTTAATGAGGAGCACCCGCAGCACGTGTATCAGCCGTACCGTTTGc                                                                              >  2:774818/1‑149 (MQ=17)
                                                                         tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGGGCGGAATATGATGAATGGGGCAACCAGCCTACTGAGGACAACCCCCATCCTGTGTAGCAGCACTACCGTCTggtat                                                                           >  2:587077/1‑145 (MQ=17)
                                                                         tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTCAATGGGGAGAACCCGCACCCTGTGTATCAGCAGTACCGGCTCCCAg                                                                           >  2:763166/1‑149 (MQ=18)
                                                                           aTCTTTATCACTGCGACCACCGGGGACTGCCGCAGGCGCTTATCAGAGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAGCAGCCTAAGGAGGAGCACCCGGCCACTGTGTATCGGCAGTCCAGCCTGCCAggg                                                                         >  2:26763/1‑149 (MQ=17)
                                                                                tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                    >  1:57982/1‑149 (MQ=255)
                                                                                   cACTGCGACCACCGGGGACTGCCGCTGGCGCTGAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagca                                                                 <  2:371822/149‑1 (MQ=255)
                                                                                    aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAGGACGGCAATACGGCGTGGAGCGCGGAATATGGTGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCat                                                                >  1:582845/1‑149 (MQ=255)
                                                                                        cGACCACCGGGGACTGCCGCGGGCGCTTATCAAGGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                                                            <  2:102614/149‑1 (MQ=255)
                                                                                         gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                           <  1:314047/149‑1 (MQ=255)
                                                                                           cccacGGGGACTGCCGGTGGCGCCTATCAGCGAAGACGGCAATACGGCGTGAAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGCGAACCCGCATCATGTGCTCCTACCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                         <  1:388780/145‑1 (MQ=14)
                                                                                              gcGGGGACTGCCGCTGGCGATTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAACGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                                      <  2:142548/148‑1 (MQ=255)
                                                                                              acGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCCTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                                      <  1:347015/148‑1 (MQ=255)
                                                                                                  ggACTGCCGCTGGCGCGTATCAGCGAAGACGGCAAGACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg                                                  <  2:61487/149‑1 (MQ=255)
                                                                                                        ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                            >  2:433325/1‑149 (MQ=255)
                                                                                                              gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATaa                                      <  1:337438/149‑1 (MQ=255)
                                                                                                                        gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACAATAACCGTCACCgg                            >  2:772746/1‑149 (MQ=255)
                                                                                                                                ggCAATACGGCGTGGAGCGCGGAATATGATGCATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGACGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATCACCGTCACCGGTACTAAGa                    >  1:748462/1‑149 (MQ=255)
                                                                                                                                 gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                   <  2:21594/149‑1 (MQ=255)
                                                                                                                                            tGGAGCGCGGAATATGATGAATGGGCCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCATCAGCATGATGCGGCATCATTGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAggg        <  1:772746/149‑1 (MQ=18)
                                                                                                                                             ggAGCGCGGAATAGGATGAAGGGGGCAACCAGTTTAATGAGGAGAACCCGCATCATGGGTATCAGCCGTGCCGTCTGCCAGGGCAGGAGCATGATGCGGAAGCAGTGCGGTACTATATCCGTCACCGGTACTACGATCCGTTGCAgggg       <  2:400714/149‑1 (MQ=17)
                                                                                                                                              gAGTGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc      <  1:227824/149‑1 (MQ=255)
                                                                                                                                                  gcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCGTGTGTATCATCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGta  >  1:423346/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTA  >  NC_000913_3_bme_pgi/526649‑526943

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