Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,659	0	T	G	25.4%	137.2 / 12.3	68	V26G (GTA→GGA)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (15/35);  new base G (17/0);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGC  >  NC_000913/729523‑729793                                                                                                                                         |                                                                                                                                        taAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACg                                                                                                                                      <  2:340684/140‑1 (MQ=255)      cTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCg                                                                                                                                 <  1:216055/140‑1 (MQ=255)       tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCgg                                                                                                                                <  2:676823/140‑1 (MQ=255)         tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTg                                                                                                                              >  1:989927/1‑140 (MQ=255)         tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTg                                                                                                                              >  1:622126/1‑140 (MQ=255)          aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGc                                                                                                                             >  1:279676/1‑140 (MQ=255)          aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCACGGTTCAACCCGGGTACGCA‑CGGTGcc                                                                                                                            >  1:1100608/1‑140 (MQ=14)              tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCCTCGGTGccccc                                                                                                                         >  2:1396520/1‑140 (MQ=255)               ttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCa                                                                                                                        <  1:934841/140‑1 (MQ=255)                tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCAc                                                                                                                       >  1:880901/1‑140 (MQ=255)                 cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCAcc                                                                                                                      >  1:1355207/1‑140 (MQ=255)                 cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGcccccccc                                                                                                                      >  2:862857/1‑140 (MQ=255)                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCCGCCGGGGGACGCATCGGTGCCCCCACCgg                                                                                                                    >  2:859598/1‑140 (MQ=255)                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCCTCCGGGCCCCCACCgg                                                                                                                    >  2:176565/1‑140 (MQ=14)                    cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGc                                                                                                                   <  2:880901/140‑1 (MQ=255)                    cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGc                                                                                                                   <  1:1396520/140‑1 (MQ=255)                    agacaggcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGc                                                                                                                   >  1:1342171/7‑140 (MQ=255)                    agacaggcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGc                                                                                                                   <  2:1342171/134‑1 (MQ=255)                     tGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCg                                                                                                                  <  1:1225181/140‑1 (MQ=255)                      gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGt                                                                                                                 <  1:1394953/140‑1 (MQ=255)                      gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGt                                                                                                                 <  2:1355207/140‑1 (MQ=255)                      gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCgg                                                                                                                 >  1:965251/1‑139 (MQ=255)                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTg                                                                                                                >  1:544699/1‑140 (MQ=255)                       tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGGATTGTCCAGGGGTCAGCCGGGGGAAGCAACGGTGGCCCCACCGGCGTg                                                                                                                >  2:417992/1‑140 (MQ=17)                        cagcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTgg                                                                                                               >  2:465983/1‑140 (MQ=255)                         agcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGc                                                                                                              <  1:176565/140‑1 (MQ=255)                             gACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGt                                                                                                          <  1:266555/140‑1 (MQ=255)                             gACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGt                                                                                                          <  2:1357828/140‑1 (MQ=255)                              aCAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGtt                                                                                                         <  1:915341/140‑1 (MQ=255)                                 aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCgg                                                                                                      >  2:407450/1‑140 (MQ=255)                                 aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTgg                                                                                                               <  2:260517/131‑1 (MQ=255)                                 aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTgg                                                                                                               >  1:260517/1‑131 (MQ=255)                                   tAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGgtg                                                                                                    <  2:515886/140‑1 (MQ=255)                                   tAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGgtg                                                                                                    >  1:200865/1‑140 (MQ=255)                                     aCCCGATAAAACAAGTATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAAGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCCGCGTGGCCTGTTCCGtgtg                                                                                                  >  1:967375/1‑140 (MQ=14)                                        cGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc                                                                                               <  2:355579/140‑1 (MQ=255)                                        cGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc                                                                                               >  2:1077500/1‑140 (MQ=255)                                        cGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc                                                                                               >  1:355579/1‑140 (MQ=255)                                          aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcg                                                                                             >  1:216423/1‑140 (MQ=255)                                          aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccccc                                                                                             >  2:126424/1‑139 (MQ=255)                                               aCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCggg                                                                                        <  1:862857/140‑1 (MQ=255)                                                cAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCA‑CGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGt                                                                                      <  2:1100608/140‑1 (MQ=255)                                                     aTGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCACTTTCCAGGGCTCAGCCGGGGTGAACATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGGCCCCGCGGGGTGaag                                                                                  >  2:465905/1‑138 (MQ=14)                                                        aGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTcc                                                                               >  1:243253/1‑140 (MQ=255)                                                         gCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCg                                                                              <  1:417992/140‑1 (MQ=255)                                                               gAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGGCCGGCCCTc                                                                        >  2:1198663/1‑140 (MQ=14)                                                                     aaaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCa                                                                  <  2:622126/140‑1 (MQ=255)                                                                     aaaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCa                                                                  <  2:989927/140‑1 (MQ=255)                                                                          cggcggCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCCTCGGTGCCCCCACCGGCGTGGCCTGTTTGGGGTGGCCCGGCGGGGGGGCGTCCGGCCCTCCGGGCAAAccc                                                             >  2:1049457/1‑139 (MQ=14)                                                                               gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                        <  1:1198663/140‑1 (MQ=255)                                                                               gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                        <  1:1049457/140‑1 (MQ=255)                                                                                cgcgTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTc                                                       <  2:794858/140‑1 (MQ=255)                                                                                      aGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCa                                                 <  1:465983/140‑1 (MQ=255)                                                                                       gggCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCaa                                                <  2:683514/140‑1 (MQ=255)                                                                                          cGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGGCGGGGGGGCGGCCCGCCCTCCCGGCAAACCCCCGCTCGGGGCAAAgg                                             >  2:789060/1‑140 (MQ=11)                                                                                                   caggTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCg                                   <  1:1162904/137‑1 (MQ=255)                                                                                                       tATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGGCC‑TTCCCGGTGa                               >  1:428966/1‑140 (MQ=255)                                                                                                       tATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGGCC‑TTCCCGGTGa                               >  1:365453/1‑140 (MQ=255)                                                                                                       tATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCCGTCAATCCGCTGCTCGGTGCAAAGGGCC‑TTCCCGGTGa                               >  1:802670/1‑140 (MQ=255)                                                                                                       tATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCCGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGGCCGGCCATCCCGTCAATCCGCTGCTCGGTGCAAAGGGCCTTTCCCggga                                >  2:167810/1‑138 (MQ=14)                                                                                                         tGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGaaa                             <  1:167810/140‑1 (MQ=255)                                                                                                               tAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACa                       <  2:243253/140‑1 (MQ=255)                                                                                                               tAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACa                       <  1:1343591/140‑1 (MQ=255)                                                                                                               tAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGTCCCGCCCTCCCGGCAAACCCCTGCCCGGTGCAAAGGGCC‑TTCCCCGGGAAACCGACa                       >  2:1027051/1‑140 (MQ=11)                                                                                                                   aTTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGc                   <  1:1077500/140‑1 (MQ=255)                                                                                                                      gTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCt                >  1:1205386/1‑140 (MQ=255)                                                                                                                            ggTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTgcccgg          >  2:784570/1‑140 (MQ=40)                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc     >  2:1237159/1‑140 (MQ=40)                                                                                                                                 aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc     >  2:123699/1‑140 (MQ=255)                                                                                                                                 aGCCGGGGCACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCCCGGTGCAAAAGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc     >  1:635193/1‑140 (MQ=255)                                                                                                                                   ccGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTg   <  1:784570/140‑1 (MQ=40)                                                                                                                                    cGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGc  <  1:367108/140‑1 (MQ=40)                                                                                                                                    cGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGc  <  1:1027051/140‑1 (MQ=40)                                                                                                                                         |                                                                                                                                        TAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGC  >  NC_000913/729523‑729793

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑