Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,053,516:1	(G)6→7	100%	coding (2663/2745 nt)	torS ←	hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,510	1	.	G	91.3%	253.8 / 18.3	80	coding (2669/2745 nt)	torS	hybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base G (36/37);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606                                                                                                      |                                                                                                  gttCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑aa                                                                                                    >  2:840040/1‑100 (MQ=255)   tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gg                                                                                                >  1:865859/1‑100 (MQ=255)   tCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gg                                                                                                >  1:1403405/1‑100 (MQ=255)     ttAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggg                                                                                              <  2:134660/100‑1 (MQ=255)       aaCTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTTTTATGCAACCATGCTTCCAGATCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggg                                                                                            <  2:1513268/100‑1 (MQ=255)         cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggg                                                                                            <  1:2345289/98‑1 (MQ=255)         cTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTGATGCAACCATGCTTCCAGTGCGGCAACACTGCGTGTAATTTCTTCGTGGGG‑AA‑‑ggggg                                                                                             <  1:2317452/97‑1 (MQ=255)           ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggggcag                                                                                        >  2:1028453/1‑96 (MQ=255)           ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggggcag                                                                                        <  2:433012/100‑5 (MQ=255)           ggCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggggcag                                                                                        >  2:2477416/1‑96 (MQ=255)            gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTCTGCAACCATGCTTCCAGAGCGGCAACCCTGCGTGTAATTTCTTCGTGTGG‑AA‑‑gggggg                                                                                            >  1:464583/1‑95 (MQ=255)            gCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGAAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑ggggggggca                                                                                        <  1:140079/99‑5 (MQ=255)             cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGAt                                                                                       <  1:593498/100‑1 (MQ=255)             cTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGAt                                                                                       >  2:2012036/1‑100 (MQ=255)                 tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATg                                                                                   <  1:1450215/100‑1 (MQ=255)                 tatTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATg                                                                                   <  1:1699167/100‑1 (MQ=255)                 tatTGATTCAAAGCGCGTTCAGGTCTTTCTGATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAAtn                                                                                   <  1:1411636/100‑2 (MQ=255)                  atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATgg                                                                                  <  1:2508158/100‑1 (MQ=255)                  atTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATgg                                                                                  <  2:384966/100‑1 (MQ=255)                   ttGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgc                                                                                 <  1:1206591/100‑1 (MQ=255)                    tgcttCAAAGCGCGTTTAGGTCTTTCGGATGGAACCATGCTTCCAGAGCGGCAGCACTGCGTGTAAGTTCTTCGTGTGG‑GAG‑GGGGGGCAGATAATGgct                                                                                <  1:1176646/97‑1 (MQ=255)                    tGATTCAAATCGCGTTCAGGTCTTTTTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgcn                                                                                <  1:1366600/100‑2 (MQ=255)                    tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgct                                                                                >  1:1352030/1‑100 (MQ=255)                    tGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgct                                                                                <  1:2316757/100‑1 (MQ=255)                     gATTCAAGTGGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGGTCGGCAACACTGCGGGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgttg                                                                               >  1:361523/1‑100 (MQ=255)                     gATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgctg                                                                               >  1:1348865/1‑100 (MQ=255)                      aTTCGAATCGCGGTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgctgc                                                                              <  1:2289425/100‑1 (MQ=255)                         cAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgctgctgc                                                                           <  1:78067/100‑1 (MQ=255)                          aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgctgctgct                                                                          >  2:2189283/1‑100 (MQ=255)                          aaaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGgctgctgct                                                                          >  1:2096989/1‑100 (MQ=255)                           aaTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTc                                                                         >  2:2236840/1‑100 (MQ=255)                            aTCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTcc                                                                        <  1:554087/100‑1 (MQ=255)                             tCGCGTTCAGGTCTTTGTTATGCAACCATGCTGCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCa                                                                       <  1:1516216/100‑1 (MQ=255)                             tCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCa                                                                       <  2:2260431/100‑1 (MQ=255)                                 gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                   <  2:797256/100‑1 (MQ=255)                                 gTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTtg                                                                   <  2:796713/100‑1 (MQ=255)                                  ttCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTtgt                                                                  >  2:1997069/1‑100 (MQ=255)                                   tCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTtgtg                                                                 <  1:909387/100‑1 (MQ=255)                                     aGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcg                                                               >  2:2128226/1‑100 (MQ=255)                                      ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgcgc                                                              >  2:2571907/1‑100 (MQ=255)                                      ggTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTgc                                                                >  1:1559207/1‑98 (MQ=255)                                       gTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCa                                                             >  2:1335183/1‑100 (MQ=255)                                        tctttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCAt                                                            <  2:9139/100‑1 (MQ=255)                                          tttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgc                                                         <  1:1625774/100‑1 (MQ=255)                                           ttcttATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGCGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAg                                                          >  1:1811721/1‑99 (MQ=255)                                               tATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTCCGTGTGGCAG‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggc                                                     >  1:1186315/1‑100 (MQ=255)                                                aTGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAgctggct                                                    >  1:1316441/1‑100 (MQ=255)                                                  gCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  2:1887376/1‑100 (MQ=255)                                                  gCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  2:1402134/1‑100 (MQ=255)                                                   cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGctggctgg                                                  >  1:1846634/1‑99 (MQ=255)                                                   cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 >  1:1979522/1‑100 (MQ=255)                                                   cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 >  1:1978507/1‑100 (MQ=255)                                                   cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1176318/100‑1 (MQ=255)                                                   cAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCtggctggc                                                 <  2:1098476/100‑1 (MQ=255)                                                     aCCATGCTTCCAGAGCGGCCACACTGCGTGTAATTTCTCCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTACTCCAGTTGTGCGCATAGCtagctggc                                                 >  1:544750/1‑98 (MQ=255)                                                     aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                               >  2:289854/1‑100 (MQ=255)                                                     aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGa                                               >  1:2453804/1‑100 (MQ=255)                                                     aCCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              >  2:113773/1‑100 (MQ=255)                                                      ccATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGat                                              <  2:233019/100‑1 (MQ=255)                                                        aTGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGatat                                            <  2:1842875/100‑1 (MQ=255)                                                          gCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAGGGGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATg                                           <  1:268271/100‑1 (MQ=255)                                                           cTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCAt                                        >  1:1725364/1‑100 (MQ=255)                                                             tCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTc                                      >  2:871461/1‑100 (MQ=255)                                                              ccAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTc                                      <  2:1885322/100‑1 (MQ=255)                                                               cAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTcc                                     >  2:1381359/1‑100 (MQ=255)                                                                 gagCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCa                                   <  2:1715226/100‑1 (MQ=255)                                                                    cGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGAc                                <  2:1316448/100‑1 (MQ=255)                                                                     ggCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:119984/1‑100 (MQ=255)                                                                     ggCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACt                               >  1:1824715/1‑100 (MQ=255)                                                                      gCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTc                              >  1:2104178/1‑100 (MQ=255)                                                                      gCAACACTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGACTGGCGATATGCATTCCCAGACTc                              >  1:254785/1‑100 (MQ=255)                                                                        aaCACTGCGTGTAATTTCCTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGa                            <  2:1261430/100‑1 (MQ=255)                                                                           acTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         >  2:366084/1‑100 (MQ=255)                                                                           acTGCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCa                         <  2:1358905/100‑1 (MQ=255)                                                                              gCGTGTAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACt                      <  2:202094/100‑1 (MQ=255)                                                                                gtgtAATTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                    >  2:645492/1‑100 (MQ=255)                                                                                   tAATTTATTCGTGTGG‑AAG‑GGGGGGCAGACAATGGCCGCTGCTCAAGTTGTCCGCATGGCTGTCTGGCGATATGCATTCCCAGACTCGAGCAACTGc                    <  1:747877/97‑1 (MQ=255)                                                                                    aaTTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCtttt                <  1:81976/100‑1 (MQ=255)                                                                                     aTTTCTTCGTGTGG‑AAG‑GGGGGGCAGATGATGGCTGCTGCTCCAGTTGTGCGCATGGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  1:2100110/100‑1 (MQ=255)                                                                                     aTTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  2:2360054/1‑100 (MQ=255)                                                                                     aTTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               <  2:2158262/100‑1 (MQ=255)                                                                                     aTTTCTTCGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCttttt               >  2:674487/1‑100 (MQ=255)                                                                                         cttcGTGTGGAAA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           >  2:1698206/1‑100 (MQ=255)                                                                                         cttcGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCt           <  2:984474/100‑1 (MQ=255)                                                                                          ttcGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTg          >  1:134813/1‑100 (MQ=255)                                                                                           tcGTGTGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTACTTTTTAGCTGa         <  1:2110070/100‑1 (MQ=255)                                                                                             gtgtGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg       <  1:2045738/100‑1 (MQ=255)                                                                                             gtgtGG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATg       >  2:2091767/1‑100 (MQ=255)                                                                                               gtgG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg     >  2:446022/1‑100 (MQ=255)                                                                                               gtgG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcg     >  1:212650/1‑100 (MQ=255)                                                                                                tgG‑AAG‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATgcgc    >  2:700406/1‑100 (MQ=255)                                                                                                 gg‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCgg  >  2:1294919/1‑100 (MQ=255)                                                                                                      |                                                                                                  GTTCTTAACTGGCTGATATTGATTCAAATCGCGTTCAGGTCTTTCTTATGCAACCATGCTTCCAGAGCGGCAACACTGCGTGTAATTTCTTCGTGTGG‑AA‑‑GGGGGGCAGATAATGGCTGCTGCTCCAGTTGTGCGCATAGCTGGCTGGCGATATGCATTCCCAGACTCGAGCAACTGCTTTTTAGCTGATGCGCGG  >  NC_000913/1053411‑1053606

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