Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,625	G→A	S70F (TCC→TTC)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,625	0	G	A	98.4%	217.4 / NA	64	S70F (TCC→TTC)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/31);  minor base C (1/0);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGG  >  NC_000913/3176378‑3176869                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ggCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCGCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAtt                                                                                                                                                                                                                                                   <  1:39720/251‑1 (MQ=255) ggCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAtt                                                                                                                                                                                                                                                   >  1:119652/1‑251 (MQ=255)  gCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGAAGAGAATTg                                                                                                                                                                                                                                                  >  1:273876/1‑251 (MQ=255)         tACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTg                                                                                                                                                                                                                                            <  1:314193/250‑1 (MQ=255)              tCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTa                                                                                                                                                                                                                                       <  2:101265/250‑1 (MQ=255)                         aaGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGCGAATTGATCCTGCGCTAAATCACCTGTcg                                                                                                                                                                                                                           >  2:71438/1‑251 (MQ=255)                         aaGCGGAAACTCGCTCATCTCACCGTGCGGCACGCCAAACCCCTGGCTATCCAGCAACAGGATTTGCCATTGCTCCCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTca                                                                                                                                                                                                                           <  2:215130/251‑2 (MQ=255)                             tGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGAc                                                                                                                                                                                                                       >  1:71075/1‑251 (MQ=255)                                 aCTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGAGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATc                                                                                                                                                                                                                   >  1:195291/1‑251 (MQ=255)                                        cagctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCGCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTc                                                                                                                                                                                                              <  2:234006/249‑1 (MQ=255)                                             caCCGTGCGGCACGCCAAACACCGGGCTATCCAGCAACAGGATTTGCCATTGCTCAACAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCg                                                                                                                                                                                                       <  2:289469/251‑1 (MQ=255)                                                       cACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCTGCATCCTGTATCGCGTTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGTCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCATCGAAATGCTGATAGGTCGCAGAGAATTGATCCTTCGATAAATCACCTGTCGCGACAATCAGGTCTAATTCGTGCTtgtgtg                                                                                                                                                                                             <  2:388379/251‑1 (MQ=255)                                                          gCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGcc                                                                                                                                                                                          >  1:230904/1‑251 (MQ=255)                                                          gCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCGGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGcc                                                                                                                                                                                          <  1:396965/251‑1 (MQ=255)                                                             aaaCACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTACCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAACTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGCTCGAATTCGTGCTGGTGTGGCCGa                                                                                                                                                                                        >  1:27238/1‑250 (MQ=255)                                                              aaCACCTGGCTAGCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCGAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAAt                                                                                                                                                                                      <  1:48257/251‑1 (MQ=255)                                                                        tATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATCGGCCGCAGAGAATTGATCCTGCGCTAAATCAACTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGc                                                                                                                                                                            >  2:117517/1‑251 (MQ=255)                                                                        tATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGc                                                                                                                                                                            >  1:217380/1‑251 (MQ=255)                                                                                 aCAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTgg                                                                                                                                                                   >  1:303474/1‑251 (MQ=255)                                                                                  cAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGt                                                                                                                                                                  <  1:363331/251‑1 (MQ=255)                                                                                      aTTTGCCATTGCTCCCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCTCCTCCAGCACCGCCTGGTAAct                                                                                                                                                              <  2:282820/251‑1 (MQ=255)                                                                                                               cgcgCTTCGCCGGGGAGATACCCGCAGCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAgg                                                                                                                                     <  2:250763/251‑1 (MQ=255)                                                                                                                     gcGCCGGGGAGATCCCCGCATCCTGTAACGCGCTGTACATCGAGGGCTGGAAATCGGGGTTGCCCGGCAGACAGACGCAGGGCGCACGAAAACTGGCGATCCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCTAATTCGTGCTGGTGTGGCCGAAACGCTTCCAGCACCGCCTGATAACTCTCCGAGGTGTTTACCCCTAACAGGGCtt                                                                                                                                 <  2:67022/248‑1 (MQ=255)                                                                                                                      cGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAACTCGTGGTTGCCCGGCCGCCAGACGCAGCGCGCACGAAAACTTGCGTTGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGCTACTGCGCTAAATAACCTGTCGCGACAATCAGGTCGAATTAGTGCTGGTGTGGACGAATCGCCTCCAGCACCGACTGGTAACTCTCCCAGGTGTTTAACCCTAACAGGGCTTCGt                                                                                                                              >  1:391586/1‑251 (MQ=255)                                                                                                                        ccGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGc                                                                                                                            <  1:307117/251‑1 (MQ=255)                                                                                                                         cGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCTGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCGGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTTGTAACTCGCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCt                                                                                                                           <  2:289471/251‑1 (MQ=255)                                                                                                                          ggggAGACACCCGCATCATGTCACGCGCAGTACATCGCGGGCTGGAAATCGAAGTTGCCCGGCAGCCAGACGACGGGCGCACGAAAACTTGCGATCCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGACTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGGGTTTACCCCTACCAGGGCTTCGTGCtt                                                                                                                          >  2:390705/1‑251 (MQ=255)                                                                                                                              agaTACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                                      <  2:230904/251‑1 (MQ=255)                                                                                                                               gaTACCCGCATCCGTTAACGCGCGGTACATGGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGATGCAGGGCGCACGAAAACTTGCGATGCCGTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCCCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgtg                                                                                                                     <  2:325921/251‑1 (MQ=255)                                                                                                                                aTACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAACTCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                                    >  1:118325/1‑251 (MQ=255)                                                                                                                                        aTCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                                            >  2:110874/1‑251 (MQ=255)                                                                                                                                         gccTGTAACGCGCAGGACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCTTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGGCGAATTCGTGCTGGTGTGGCCGAAGCGCCTCCAGCACCGCCTTGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                                           <  2:322772/250‑1 (MQ=255)                                                                                                                                          ccTGTAACGCGCGGTAAATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAGCACCTGTCGCGTCAATAAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGCTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGTTGa                                                                                                          <  2:159235/251‑1 (MQ=255)                                                                                                                                           cTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                                           >  1:133792/1‑249 (MQ=255)                                                                                                                                             ggaaCGCGCGGTACACCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCGTCCGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCTAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAGCTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTTTGCAAACAGGTGAg                                                                                                         <  2:237099/247‑1 (MQ=255)                                                                                                                                                   cgcTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGCTGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                                  >  1:185296/1‑250 (MQ=255)                                                                                                                                                     ctgtacagtgagggCTGGTAATGGTGGGTGCCCGGAAGCCCGACGCAGAGAGCAAGCAAACTTACGCTTCCATCCGCGAAATGACGATAGGACGCAGAGCATTGATCCTGCGATAAATCAACTGTCGCGACAATCAGGTCGAATTCGTGCTCGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCGCAGGTGTTTTACCCTCACAGGGCTTCGTGCTTTTGTGCACACAGGTGAGTGACGGAAAt                                                                                               >  2:148678/12‑251 (MQ=255)                                                                                                                                                         aCATCGCGGTCTGGAAAGCGTGGTTGCCCGGCAGCAAGACGCAGGGCGCACGAAAACTTGCGATGCCTTAAGCAAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTAGCGATAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTTACTCTGCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                                                           <  2:298190/251‑1 (MQ=255)                                                                                                                                                             cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaaa                                                                                       <  1:234060/251‑1 (MQ=255)                                                                                                                                                             cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                        >  1:106325/1‑250 (MQ=255)                                                                                                                                                             cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                        >  2:97595/1‑250 (MQ=255)                                                                                                                                                             cgcgGGCTGGAAATCGTGGTGGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCATCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCTCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGAGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaaa                                                                                       <  2:161818/251‑1 (MQ=255)                                                                                                                                                                     ggAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAct                                                                               <  1:331008/251‑1 (MQ=255)                                                                                                                                                                     ggAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAc                                                                                <  2:106325/250‑1 (MQ=255)                                                                                                                                                                       aaaTCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAAAAGGTGAGTGTCGGTAATTTGTAAAATCCTGActct                                                                             >  2:240544/1‑251 (MQ=255)                                                                                                                                                                            gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCt                                                                        >  1:108141/1‑251 (MQ=255)                                                                                                                                                                                ttGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGAAAGGCCGCAGAGAATTGATCCTGCGCTAAATCAACTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGAGTTTACCCCTAACCGGGCTTCGTGCTTTTGTGCAAACAGGTGTGTGTCGGTAATTTGTAAAATCCTGACTATGGCCTCAcc                                                                    >  2:241461/1‑251 (MQ=255)                                                                                                                                                                                    ccGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccagcc                                                                <  1:110874/251‑1 (MQ=255)                                                                                                                                                                                             aGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAgg                                                       >  2:132960/1‑251 (MQ=255)                                                                                                                                                                                                 gCAGGGCGCACGAAAACTTTCGATGCCTTCAGCGAAATGGTGATAGGCCGCAGAGAATTGATCCTGCGGTAAATCACCTGTCGCGACAAACAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATGTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTa                                                   <  2:331035/251‑1 (MQ=255)                                                                                                                                                                                                   aGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTACATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAAc                                                 >  1:236914/1‑251 (MQ=255)                                                                                                                                                                                                      gcgcACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGGGAATTGATCCTGCGCTAAATCACCGGGCGCGACAATCAGGTCGAATTCTTGCGGGGGTGGCCGAATCGCCTCTAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACTGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTGGTAAAATCCTGTCTCTGGCCTCACTAGCCAGAGTAAGGGTTAACAgg                                              <  2:369735/251‑1 (MQ=255)                                                                                                                                                                                                         cACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGACTTGATCCTGCGCTAAATCAACTGTCGCGACAATCAGGTCGAATTCGTGCTGTTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCtt                                           >  1:249131/1‑251 (MQ=255)                                                                                                                                                                                                                aCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTGTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGTAAGGGTTAACAGGCTTTCCAAAt                                    <  2:71075/251‑1 (MQ=255)                                                                                                                                                                                                                  ttGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATgg                                  >  1:154077/1‑251 (MQ=255)                                                                                                                                                                                                                         gCCTTCAGCGAAATGCTGATAGGGCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtt                           <  2:108141/251‑1 (MQ=255)                                                                                                                                                                                                                               aGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGttt                     >  2:344930/1‑251 (MQ=255)                                                                                                                                                                                                                               aGCGAAATGCTGATAGGCCGCAGAGAATAGCTCCTGCACCACATCACCTGTCGCGACAAACAGGTCGAATGCGTGCTGGAGTGGCCGAAACGCATCCAGCACCGCCTGGAAACTCTCCCAGGTGTTTACACCTAACAGGGCTTCGAGCTTTTGTGCAACCAGGTGCGTGTCGGAAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACGGGCTTTCCAAATGGTGTCCTTAGGttc                     >  2:328433/1‑250 (MQ=255)                                                                                                                                                                                                                                   aaaTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacg                 >  1:24117/1‑251 (MQ=255)                                                                                                                                                                                                                                        cTGATAGGCCGCAGAGAATTGATCATGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATAGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGCAGGGTTAACAGGATTTCCAAATGGGGTCCTTAGGTTTCACGACGCt            >  2:34781/1‑251 (MQ=255)                                                                                                                                                                                                                                              ggCCGCAGAGAATTGAACCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGACTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTCCGTGCTTTTGTGCAAGCAGGTGAGTGTCGGTAATTTGTAAGATCCTGACTCTGGCCTCACCCGCCAGAGGAAGGGATAACAGGCTTTCCAAATGGAGTCCTTAGGTTTCACGACGCtaataa       >  2:129952/1‑250 (MQ=255)                                                                                                                                                                                                                                              ggCCGCAGAGAAGTGATTCTGCGCTAAATCTCCGGTCGCGACACTCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCGCCAGCACCGCCTGGTAACTCTCCCAGGTGGTTACCCCGAACATGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGTCTCTGGTCTCACCAGCCAGAGGAAGGGTTAACATGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATaaa      <  1:19192/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAAcc    >  1:221843/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTACTAAACCgg  >  2:6523/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGG  >  NC_000913/3176378‑3176869

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