Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,973,377	C→T	A650T (GCG→ACG)  A650T (GCG→ACG)	cheA ← cheA ←	chemotaxis protein CheA chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,377	0	C	T	98.3%	190.6 / ‑4.2	59	A650T (GCG→ACG)  A650T (GCG→ACG)	cheA cheA	chemotaxis protein CheA chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/31);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGC  >  NC_000913/1973132‑1973614                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              caacaaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATAc                                                                                                                                                                                                                                          <  1:245040/251‑1 (MQ=255)   acaaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACgt                                                                                                                                                                                                                                        >  1:295228/1‑251 (MQ=255)    caaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACgtt                                                                                                                                                                                                                                       <  2:247341/251‑1 (MQ=255)    caaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACgtt                                                                                                                                                                                                                                       <  1:246133/251‑1 (MQ=255)    caaTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACgtt                                                                                                                                                                                                                                       <  1:312583/251‑1 (MQ=255)      aTCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACgttgt                                                                                                                                                                                                                                     <  1:327953/251‑1 (MQ=255)         ggCACAATCACGCCGCGCAGATCCGCGACGCCTTTGATAAACGCTGGCGTGTTCGCACTCCGTGTTACCTGATCGTAGCCACGGATCTCCTGGACTTTCAGGATATAAATAAAGTACTATTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGAAGGCTCGCTGGCCCGCTTTGTTACATTCGTCATACCGGTCATATAGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTcg                                                                                                                                                                                                                                  >  2:422026/1‑251 (MQ=255)             caATAACGCCGCGCAGAGTCGTGACGACAGTGATAAACGCTGGCGTGTTCGCAAGCCGTGTTACCTGATCGTAGCCACGGATCTCCTTCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATGCCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTTCATTCGTCATACCGGTCATATTGTTGCCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGt                                                                                                                                                                                                                              <  1:470616/251‑1 (MQ=255)               aTAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCATGCGGCGGTGTTCGTCATACGTTGTTCGCGGttt                                                                                                                                                                                                                            <  1:552778/251‑1 (MQ=255)                 aaagccgcGCAGATTCGTGGCGACTTTGATAAACGCTGGCGTTTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACGCTTCATCGCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGTCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTTTTCGCGGTTTa                                                                                                                                                                                                                           <  2:74596/247‑1 (MQ=255)                     ccgcGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGcc                                                                                                                                                                                                                      >  1:69956/1‑251 (MQ=255)                       gcgcAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTg                                                                                                                                                                                                                    <  1:262762/251‑1 (MQ=255)                          cAGATTAGTGACGCATTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCAAGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATAGCCTGCa                                                                                                                                                                                                                  >  2:417948/1‑250 (MQ=255)                                                  cGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGc                                                                                                                                                                                         >  2:23101/1‑251 (MQ=255)                                                         gtgatCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTACTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGAGGAGACATCAACACTCAGAGCTACGCTg                                                                                                                                                                                  >  2:637329/5‑251 (MQ=255)                                                               gCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTcc                                                                                                                                                                            >  2:316194/1‑250 (MQ=255)                                                                 aaTCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCGACGCTGCCGTCGcc                                                                                                                                                                          <  1:170837/251‑1 (MQ=255)                                                                      gtgtTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGa                                                                                                                                                                      >  1:104627/1‑250 (MQ=255)                                                                        gtTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATg                                                                                                                                                                   <  2:168170/251‑1 (MQ=255)                                                                          tACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATgg                                                                                                                                                                  <  1:491429/250‑1 (MQ=255)                                                                                aTCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCACGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATCCCGGTCCTATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCgcagc                                                                                                                                                           >  2:397721/1‑251 (MQ=255)                                                                                  cGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCgcagca                                                                                                                                                          <  1:457059/250‑1 (MQ=255)                                                                                      gCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAAt                                                                                                                                                     >  2:367276/1‑251 (MQ=255)                                                                                        cACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGGCTGACGTCTCGCTGGCCATCTTTGTTACATTCGTCATGCCTGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGCCGCCAAGAGTGGTCGCAGCAGAAATGc                                                                                                                                                   <  2:55363/251‑1 (MQ=255)                                                                                         aCGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGcc                                                                                                                                                  <  1:290286/251‑1 (MQ=255)                                                                                           ggATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCgg                                                                                                                                                >  1:587253/1‑251 (MQ=255)                                                                                              tctcCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCgggg                                                                                                                                              >  1:49213/1‑250 (MQ=255)                                                                                                 ccTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCggg                                                                                                                                               <  2:549430/246‑1 (MQ=255)                                                                                                  cTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACttt                                                                                                                                         >  1:107103/1‑251 (MQ=255)                                                                                                    gcaagttCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACTTTCGTCATACCGGTCATATTGTTACCTTTTCACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAATGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGCACTTTg                                                                                                                                        <  2:405358/245‑1 (MQ=255)                                                                                                     cACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCg                                                                                                                                      >  1:370465/1‑251 (MQ=255)                                                                                                         ttCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTTGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATa                                                                                                                                   <  1:646256/250‑1 (MQ=255)                                                                                                            aGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTa                                                                                                                               <  1:409310/251‑1 (MQ=255)                                                                                                             ggATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTAc                                                                                                                              <  1:367276/251‑1 (MQ=255)                                                                                                                       aCCGTACTCTTCATCACCAAGGGTAAATACCATAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGtt                                                                                                                    <  1:19098/251‑1 (MQ=255)                                                                                                                          gTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGttttt                                                                                                                 >  1:421404/1‑251 (MQ=255)                                                                                                                                  catcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACcaca                                                                                                          <  2:331047/250‑1 (MQ=255)                                                                                                                                   atcaCCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAAc                                                                                                        >  1:580681/1‑251 (MQ=255)                                                                                                                                         aaGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATCGTTACTTTCAAGGTTTTTAAc                                                                                                              >  2:433143/1‑239 (MQ=255)                                                                                                                                         aaGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAAc                                                                                                              <  1:433143/239‑1 (MQ=255)                                                                                                                                          aGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTa                                                                                                                >  1:463645/1‑236 (MQ=255)                                                                                                                                          aGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTa                                                                                                                <  2:463645/236‑1 (MQ=255)                                                                                                                                                aaTACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCa                                                                                           <  2:69956/251‑1 (MQ=255)                                                                                                                                                 aTACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCa                                                                                           >  2:512918/1‑250 (MQ=255)                                                                                                                                                        aaaTTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCaattaattt                                                                                   >  2:168117/1‑250 (MQ=255)                                                                                                                                                         aaTTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGa                                                                                  <  1:515950/251‑1 (MQ=255)                                                                                                                                                                              cttgCCAGATTTGCTACATTCGTCATCCCGGTCATATCGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCGGTGCCACGCTGCCGTCGCCAAGACTGGTCGCAGCAGAAATGCAGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCACTTAATTGATACACCAGCAAGGCGTAncgg                                                             >  2:609917/4‑251 (MQ=255)                                                                                                                                                                                ggCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcggcg                                                           <  2:295228/251‑1 (MQ=255)                                                                                                                                                                                ggCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcggcg                                                           <  2:370465/251‑1 (MQ=255)                                                                                                                                                                                ggCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcggcg                                                           >  1:176721/1‑251 (MQ=255)                                                                                                                                                                                ggCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcggcg                                                           <  2:107103/251‑1 (MQ=255)                                                                                                                                                                                ggCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAgcggc                                                            <  1:409277/250‑1 (MQ=255)                                                                                                                                                                                    aGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCg                                                       >  1:637869/1‑251 (MQ=255)                                                                                                                                                                                                           ggTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAAt                                >  2:128026/1‑251 (MQ=255)                                                                                                                                                                                                              cATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTcc                             >  2:235563/1‑251 (MQ=255)                                                                                                                                                                                                                  ttGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACcaca                                                                                                          <  1:232317/170‑1 (MQ=255)                                                                                                                                                                                                                  ttGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACcaca                                                                                                          >  2:232317/1‑170 (MQ=255)                                                                                                                                                                                                                                   aTTCAGGCGGCGGTGCTCGTCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGCCGCCAAGACTGGTCGCCGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCCAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCACTGGGTGGCTTCGGTTTTctc        >  2:377132/1‑249 (MQ=255)                                                                                                                                                                                                                                         gcggcgGTGTTCGCCCTACGTTGTTCGCGGTTTATCGCCCGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCCCCCAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGCTAGTTACTTTCAAGGGTTTTAACCCAAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGGTTTCGCGCCcgc  >  2:401483/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              CAACAATCGGCACAATAACGCCGCGCAGATTCGTGACGCCTTTGATAAACGCTGGCGTGTTCGCAATCCGTGTTACCTGATCGTAGCCACGGATCTCCTGCACTTTCAGGATATCAATACCGTACTCTTCATCACCAAGGGTAAATACCAGAAATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGC  >  NC_000913/1973132‑1973614

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