Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	526,792	0	A	C	21.0%	89.8 / 24.0	65	S1178R (AGC→CGC)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base C (5/8);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGC  >  NC_000913_3_bme_pgi/526644‑526939                                                                                                                                                     |                                                                                                                                                    gTGGCTTGCGCGTTGCGGGCTGCCGGTGGAGCAACTGGCCAGCCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGa                                                                                                                                                     <  1:302321/149‑1 (MQ=255) gTGGCTTGCGCAGTGCGGGCAGACGGTGGAGCAACTGGCAAGACAGGTGGAGCCGGAATACACACCGGCGCGAACAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGa                                                                                                                                                     <  1:144143/149‑1 (MQ=255)  tGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                                    >  2:569724/1‑149 (MQ=255)  tGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAg                                                                                                                                                    >  1:473109/1‑149 (MQ=255)    gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcg                                                                                                                                                  >  1:124651/1‑149 (MQ=255)      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                                                                                                                                 <  1:81751/149‑1 (MQ=255)                     tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                 <  2:553286/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:656995/149‑1 (MQ=255)                                         aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                             >  1:398238/1‑149 (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)                                                                                                                                                     |                                                                                                                                                    GTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGC  >  NC_000913_3_bme_pgi/526644‑526939

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