Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,016,908	C→A	6.0%	R119S (CGC→AGC)	fliI →	flagellum‑specific ATP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,016,908	0	C	A	6.0%	242.1 / 2.8	67	R119S (CGC→AGC)	fliI	flagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base C (24/39);  new base A (2/2);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCC  >  NC_000913/2016769‑2017047                                                                                                                                            |                                                                                                                                            tAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTc                                                                                                                                             >  1:843467/1‑140 (MQ=255)       tCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTc                                                                                                                                       >  2:873239/1‑140 (MQ=255)         aaCGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTg                                                                                                                                     >  1:319739/1‑140 (MQ=255)         aaCGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTg                                                                                                                                     >  1:1166944/1‑140 (MQ=255)         aaCGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTg                                                                                                                                     >  1:1166990/1‑140 (MQ=255)           cGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGa                                                                                                                                   <  2:418299/140‑1 (MQ=255)           cGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGa                                                                                                                                   <  2:1179024/140‑1 (MQ=255)              cTGTTTTTAATGCCGCTGGTGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACgg                                                                                                                                <  1:451570/140‑1 (MQ=255)              cTGTTTCTCATGCCGCTGGAGGTAGTCGACGGTGTCCTGCCCGGCGCGCGTGGATATGCCAAAAACCTTTCGGCAGAAGGGCTGCAAAGCGGATATCAGTTGCCGCTCTTTCCGGCGTTCTTAGGTCGCGTTCTGGTCgg                                                                                                                                <  2:1089213/140‑1 (MQ=255)                           cGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCg                                                                                                                   >  2:987608/1‑140 (MQ=255)                           cGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCg                                                                                                                   >  1:943222/1‑140 (MQ=255)                             cgggaggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCt                                                                                                                 <  2:571647/138‑1 (MQ=255)                                     aGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTAGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGcc                                                                                                         >  2:443591/1‑140 (MQ=255)                                     aGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTAGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGcc                                                                                                         >  2:443201/1‑140 (MQ=255)                                         gAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGcc                                                                                                     <  1:1058354/140‑1 (MQ=255)                                         gAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGcc                                                                                                     <  2:843467/140‑1 (MQ=255)                                         gAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGcc                                                                                                     <  1:1148768/140‑1 (MQ=255)                                                tCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTccccc                                                                                              >  2:640376/1‑140 (MQ=255)                                                tCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTccccc                                                                                              >  2:830372/1‑140 (MQ=255)                                                tCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTccccc                                                                                              >  1:810531/1‑140 (MQ=255)                                                    gCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATa                                                                                          <  2:159028/140‑1 (MQ=255)                                                           gcgcgTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGa                                                                                   <  1:921245/140‑1 (MQ=255)                                                              cgTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAAc                                                                                >  2:297422/1‑140 (MQ=255)                                                                   ttATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTg                                                                           <  2:402500/140‑1 (MQ=255)                                                                               cttttAGGCAGAAGGGCTGCAATGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAAGCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTAccc                                                               <  1:610965/138‑1 (MQ=255)                                                                                ctttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTAcccc                                                              <  2:386650/139‑1 (MQ=255)                                                                                 tttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                                             <  1:408218/140‑1 (MQ=255)                                                                                 tttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                                             <  1:1013860/140‑1 (MQ=255)                                                                                 tttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                                             <  2:837694/140‑1 (MQ=255)                                                                                 tttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                                             <  1:296814/140‑1 (MQ=255)                                                                                 cttCGGCTGTCGGGCTGCACAGCGGCTAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                                             <  2:684903/139‑1 (MQ=255)                                                                                        aGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTa                                                      <  2:15351/140‑1 (MQ=255)                                                                                          aaGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAAc                                                    >  1:891507/1‑140 (MQ=255)                                                                                          aaGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTAGTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAAc                                                    >  1:955190/1‑140 (MQ=255)                                                                                               cTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                                               <  1:743000/140‑1 (MQ=255)                                                                                               cTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                                               <  1:635691/140‑1 (MQ=255)                                                                                               cTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                                               <  1:159494/140‑1 (MQ=255)                                                                                                tGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTg                                              <  1:273000/140‑1 (MQ=255)                                                                                                        ggCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTAGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTac                                      <  1:443591/140‑1 (MQ=255)                                                                                                        ggCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTAGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTac                                      <  1:443201/140‑1 (MQ=255)                                                                                                         gCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTaca                                     <  2:472800/140‑1 (MQ=255)                                                                                                         gCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTaca                                     <  2:1178330/140‑1 (MQ=255)                                                                                                          cAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTacac                                    >  2:360311/1‑140 (MQ=255)                                                                                                          cAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTacac                                    >  2:360268/1‑140 (MQ=255)                                                                                                          cAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTacac                                    >  2:1180170/1‑140 (MQ=255)                                                                                                          cAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTacac                                    <  2:810531/140‑1 (MQ=255)                                                                                                           aaGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACAcc                                   <  1:908547/140‑1 (MQ=255)                                                                                                            aGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCg                                  >  1:1160339/1‑140 (MQ=255)                                                                                                            aGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCg                                  >  1:1160506/1‑140 (MQ=255)                                                                                                              cAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGAt                                <  1:1057082/140‑1 (MQ=255)                                                                                                              cAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGAt                                <  1:792181/140‑1 (MQ=255)                                                                                                              cAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGAt                                <  1:1197792/140‑1 (MQ=255)                                                                                                                 ttGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGa                             >  1:871415/1‑140 (MQ=255)                                                                                                                 ttGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGa                             >  1:149155/1‑140 (MQ=255)                                                                                                                 ttGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGa                             >  1:122645/1‑140 (MQ=255)                                                                                                                   gCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAAc                           <  1:1081142/140‑1 (MQ=255)                                                                                                                   gCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAAc                           <  1:882850/140‑1 (MQ=255)                                                                                                                   gCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAAc                           <  1:330718/140‑1 (MQ=255)                                                                                                                   gCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAAc                           <  1:298038/140‑1 (MQ=255)                                                                                                                     cGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACAt                         <  1:1032115/140‑1 (MQ=255)                                                                                                                       cTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACAtgt                       >  2:1064362/1‑140 (MQ=255)                                                                                                                               ggCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGaca               <  1:902568/140‑1 (MQ=255)                                                                                                                                gCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGacac              <  2:1098328/140‑1 (MQ=255)                                                                                                                                 cGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACAcc             <  1:167948/140‑1 (MQ=255)                                                                                                                                     tttaGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCg         <  1:1217092/139‑1 (MQ=255)                                                                                                                                      tgaGGTCGCGGTCTGGACGGCAGCGGTAAACCGCTCGATGGCCAGCCCTCCCCCGTTACGACGGAAACCGGTGCGCTGATAACCCCGCCAATAAACCCTTAGCAACGATCACCGATGGAACAAGTGCTGGACACCGGCGt        >  2:1227724/3‑140 (MQ=255)                                                                                                                                          gacgcgTTCTGGACGGCAGCGCTAAACCGCACGATGGCCTGCCCTCCCCCGATACGACTGAACCCGGTGCGCTGATTACCCCGCCATGTAACCCGTTGCAACGTACACCGCTTGAACATGTGCTGGACACCGGCGTgcgc    <  2:94065/138‑1 (MQ=255)                                                                                                                                           tCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGTTACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACGGATTGAACATGTGCTGGACACCGGCGTGCGcc   >  2:799012/1‑140 (MQ=255)                                                                                                                                           tCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGcc   >  2:503622/1‑140 (MQ=255)                                                                                                                                            cgcgTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGccc  <  2:573108/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            TAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCC  >  NC_000913/2016769‑2017047

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