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

ATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTT  >  NC_000913/729522‑729797                                                                                                                                          |                                                                                                                                           ataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTAc                                                                                                                                          <  2:1052646/140‑1 (MQ=255)   aaaTCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGc                                                                                                                                        >  1:718638/1‑140 (MQ=255)        tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGGATTGTCCAGGGAACAGCGAGGGGACGACTCgg                                                                                                                                   >  2:1682670/1‑140 (MQ=14)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTGCCCGGGTTCAGCCCGGGGACACATCggg                                                                                                                                  >  2:1250138/1‑139 (MQ=14)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGt                                                                                                                                  <  1:1293614/140‑1 (MQ=255)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGt                                                                                                                                  >  1:665945/1‑140 (MQ=255)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGt                                                                                                                                  >  1:919422/1‑140 (MQ=255)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATccgg                                                                                                                                  >  2:504495/1‑137 (MQ=255)         gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAATGTCCAGGGCTCAGCCGGGGGACGCATCGGt                                                                                                                                  >  2:1036993/1‑140 (MQ=255)          tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTg                                                                                                                                 >  1:1174085/1‑140 (MQ=255)           aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGc                                                                                                                                >  1:1129396/1‑140 (MQ=255)            gggTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGcc                                                                                                                               <  1:625546/140‑1 (MQ=255)                ttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCa                                                                                                                           <  1:1036993/140‑1 (MQ=255)                ttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCa                                                                                                                           <  1:1698663/140‑1 (MQ=255)                 tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGTGCGGAAAACCGGCGGCGCGTCAGGTCGACAGCACGCAGTATGGCGGTAGCATTGTCAACGGTTGGGGCGGGGGAGGAATAGGTGCCCCCAc                                                                                                                          >  2:696768/1‑140 (MQ=14)                 tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGGATTGTCCAGGGCTTCGCCGGGGAACGCATCGGTGCCCCCAc                                                                                                                          >  1:1045785/1‑140 (MQ=14)                 tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGccccccc                                                                                                                          >  2:495765/1‑138 (MQ=255)                  cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGcccccccc                                                                                                                         >  1:383724/1‑140 (MQ=255)                  cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCAcc                                                                                                                         >  2:1503180/1‑140 (MQ=255)                  cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCAcc                                                                                                                         >  2:1808406/1‑140 (MQ=255)                    ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGTTCAGCCGGGGGACGCATCGGGGCCCCCACCgg                                                                                                                       >  2:1298328/1‑140 (MQ=255)                    ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCCTCGGTGCCCCCACCgg                                                                                                                       >  2:514897/1‑140 (MQ=255)                     cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGc                                                                                                                      <  1:1503180/140‑1 (MQ=255)                     cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACACATCGGTGCCCCCACCGGc                                                                                                                      <  1:1250138/140‑1 (MQ=255)                       gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGt                                                                                                                    <  1:1298328/140‑1 (MQ=255)                        tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTg                                                                                                                   <  1:1511244/140‑1 (MQ=255)                          agcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGc                                                                                                                 <  2:372883/140‑1 (MQ=255)                          agcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGc                                                                                                                 <  2:665945/140‑1 (MQ=255)                              gACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGt                                                                                                             <  1:834808/140‑1 (MQ=255)                              gACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGt                                                                                                             >  1:429621/1‑140 (MQ=255)                               aCAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGtt                                                                                                            >  2:549749/1‑140 (MQ=255)                               aCAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGtt                                                                                                            <  2:1174085/140‑1 (MQ=255)                               aCAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGCAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCCGCGGGGCCCGtt                                                                                                            >  2:1262848/1‑140 (MQ=17)                                  aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCgg                                                                                                         >  1:1529769/1‑140 (MQ=255)                                  aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCCGCGTGGCCCGTTCgg                                                                                                         >  2:346966/1‑140 (MQ=255)                                      aCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGtgtg                                                                                                     <  2:1129396/140‑1 (MQ=255)                                           aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcg                                                                                                >  1:838493/1‑140 (MQ=255)                                           aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcg                                                                                                >  1:469324/1‑140 (MQ=255)                                           aTAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACATTGTCCAGGGGTTAGCCGG‑‑‑‑C‑C‑‑CGGGCTCCCCACCGGCGTGGGCTGTTTGGTGTGCCCCCGCggggg                                                                                         >  2:370775/1‑139 (MQ=9)                                             aaaaCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc                                                                                              <  2:614584/140‑1 (MQ=255)                                                  aaGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGt                                                                                         <  1:1262848/140‑1 (MQ=255)                                                     gATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGGCTGTTCGGTGTGCCCCGGCgggggggc                                                                                      >  2:1449736/1‑136 (MQ=255)                                                         aGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCCGCCGGGGAACGCACCCGTGCCCCCAACCGCGTGGGCTGGTCGGTGTGCCCCGGAGGGGTGACGTcc                                                                                  >  2:1123442/1‑140 (MQ=14)                                                         aGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACATcc                                                                                  >  1:219867/1‑140 (MQ=14)                                                          gCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCg                                                                                 <  1:346523/140‑1 (MQ=255)                                                             gATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGCCGAGGC‑‑‑‑‑T‑‑CCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGt                                                                       <  1:370775/140‑1 (MQ=11)                                                              aTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCa                                                                             <  1:1449736/140‑1 (MQ=255)                                                               tGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCAt                                                                            <  2:412269/140‑1 (MQ=255)                                                                gAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGGCGGCCCGCCacc                                                                           >  2:1458752/1‑138 (MQ=14)                                                                 aGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATcc                                                                          <  1:660119/140‑1 (MQ=255)                                                                 aGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATcc                                                                          <  2:278662/140‑1 (MQ=255)                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGTGGCGTCCCGCCATcccg                                                                        >  1:1396773/1‑138 (MQ=14)                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGGCCCCGCGGGGGGGCGTCCGGCCCTcccg                                                                        >  2:1003331/1‑138 (MQ=14)                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGGCCGGCCATcccg                                                                        >  1:1448714/1‑138 (MQ=255)                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGTCCGGCCCTcccg                                                                        >  2:122643/1‑138 (MQ=14)                                                                   cGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCCTCGGTGCCCCCCCCGGCGTGGGCTGTTCGGTGTGCCCCGGCGGGGGGGCGTCCCGGCATcccg                                                                        >  2:1563916/1‑138 (MQ=14)                                                                     gAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTc                                                                      <  2:1102988/140‑1 (MQ=255)                                                                     gAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGACGTCCGGCCCTCCCGTc                                                                      >  2:132338/1‑140 (MQ=14)                                                                           cggcggCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCg                                                                <  2:429621/140‑1 (MQ=255)                                                                               ggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgc                                                            <  1:1682670/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:1149125/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:469324/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:146493/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:838493/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:1631736/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  1:495765/140‑1 (MQ=255)                                                                                gcgcGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCgctgct                                                           <  2:1719131/140‑1 (MQ=255)                                                                                 cgcgTCAGGGCGGCATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCAGCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGGCGTCCGGCCATCCGGTCAATCCGCTGCTc                                                          <  1:696768/140‑1 (MQ=255)                                                                                       aGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCa                                                    <  2:628617/140‑1 (MQ=255)                                                                                       aGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCa                                                    <  1:175512/140‑1 (MQ=255)                                                                                       aGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCa                                                    <  1:1481539/140‑1 (MQ=255)                                                                                                 gACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTc                                          <  1:549749/140‑1 (MQ=255)                                                                                                   cGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGGGCAAAAGGCCTTccc                                        >  1:1675449/1‑140 (MQ=255)                                                                                                         aTGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGGGACGGCCGGCCCTCCGGGCAATCCCCTGCTCGGGGCAAAGGGCCTTCCCggggga                                  >  2:1165280/1‑136 (MQ=11)                                                                                                           ggCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAc                                <  1:1003331/140‑1 (MQ=255)                                                                                                           ggCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGGCGGGGGGGCGTCCCGCCCTCCCGGCAATCCCCTGCTCGGGGCAAAAGGCCTTCCCCGGGaaaa                                >  2:387532/1‑139 (MQ=11)                                                                                                                tAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACa                           <  2:1103506/140‑1 (MQ=255)                                                                                                                    aTTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGc                       <  2:1055578/140‑1 (MQ=255)                                                                                                                     ttGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGcc                      <  1:387532/140‑1 (MQ=255)                                                                                                                       gTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCt                    >  1:594714/1‑140 (MQ=255)                                                                                                                          cAGGGTTCAGCCGGGGTACGCATCGGTG‑CCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGATGAAACCGACATCGCCCTgccc                >  1:1231090/1‑140 (MQ=255)                                                                                                                             ggTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgcccgg              >  2:972826/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  2:285920/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  2:1004274/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  1:1264123/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  2:721201/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  1:1660228/1‑140 (MQ=40)                                                                                                                                  aGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGc         >  1:516718/1‑140 (MQ=25)                                                                                                                                       gggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCg    >  2:1145841/1‑140 (MQ=255)                                                                                                                                        ggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGt   <  2:1660228/140‑1 (MQ=40)                                                                                                                                         gTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGtt  >  1:1559778/1‑140 (MQ=40)                                                                                                                                          |                                                                                                                                           ATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTT  >  NC_000913/729522‑729797

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