Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,508,789	A→T	6.7%	M109K (ATG→AAG)	PP_2203 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,508,789	0	A	T	6.7%	279.8 / 3.0	99	M109K (ATG→AAG)	PP_2203	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (61/29);  new base T (2/6);  total (63/36)
Fisher's exact test for biased strand distribution p-value = 2.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGAC  >  NC_002947/2508644‑2508937                                                                                                                                                  |                                                                                                                                                       cggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑c                                                                                                                                                    <  8:162118/149‑1 (MQ=255)  ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑cc                                                                                                                                                   >  5:18341/1‑149 (MQ=255)   gcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑ccc                                                                                                                                                  <  7:184308/149‑1 (MQ=255)     aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGTCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCgg                                                                                                                                                <  6:209416/149‑1 (MQ=255)        aaaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCG‑CCCGGCat                                                                                                                                             <  1:238956/149‑2 (MQ=255)         aaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCAcc                                                                                                                                            >  4:122262/1‑149 (MQ=255)              gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAg                                                                                                                                       >  7:6810/1‑149 (MQ=255)              gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAg                                                                                                                                       <  5:190012/149‑1 (MQ=255)                tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGc                                                                                                                                     >  8:162235/1‑149 (MQ=255)                      aCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGCCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGa                                                                                                                               <  2:226046/149‑1 (MQ=255)                      aCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTCTCCTGGCCATCCTTGGTCACCCCGATGACCTGGAAGGCATCCTTCTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCATGcccgt                                                                                                                                 >  8:304170/1‑145 (MQ=255)                       cAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAg                                                                                                                              <  3:122262/149‑1 (MQ=255)                       cAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGGCACCCATCTCCATGCCGGGCGAGCCCATGGTCATG‑CCCGGCACCGCCAGGCCCTTGAg                                                                                                                              <  3:54548/149‑1 (MQ=255)                         gcgcCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGt                                                                                                                            >  8:204454/1‑149 (MQ=255)                          cgcCCCTTGGCTGAGCTACTGTAGTTGGCGACCAGGGTGTCCTGGCCTTCCTGGGTCACCCCGATGCCCTGGTAGGCATCCTTGTGGTCACCCATCTCAATGCCGGGCGAGCCCTTGGGCTTG‑CCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                                           <  4:53802/149‑1 (MQ=255)                           gcCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCCCCGCCATGCCCTTGAGGTCg                                                                                                                          >  8:243604/1‑149 (MQ=255)                             cccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAg                                                                                                                              >  7:138989/1‑143 (MQ=255)                                    cTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTg                                                                                                                 <  7:207103/149‑1 (MQ=255)                                        gCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAg                                                                                                                              >  5:223609/1‑132 (MQ=255)                                         cATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGtt                                                                                                                        >  7:296442/1‑137 (MQ=255)                                             agtagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                                        >  7:52073/1‑149 (MQ=255)                                             agtagtTGGCGACCACGGTGGCCTGCCCACCCGGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCAGCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                                        <  7:102242/149‑1 (MQ=255)                                               tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                                           >  4:50806/1‑128 (MQ=255)                                               tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                                           <  3:50806/128‑1 (MQ=255)                                                   tGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagca                                                                                                          >  4:283170/1‑141 (MQ=255)                                                   tGGCGACCACGGGGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagca                                                                                                          <  3:283170/141‑1 (MQ=255)                                                   gggCGCCCCCGGGGTCCGGCCCATGCTGGGCCCCCCGTTTCACGTGGTAGGCACCATTGTGGCCACCAATCTCCATGCCGGGCGCGCCCTGGGGCTGG‑CCGGGCACCGCCAGGCCCTTGAGGTCTTTGCGTTTGGCCAGCAGGCGCAc                                                                                                   <  6:18341/147‑1 (MQ=255)                                                     gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                                  >  5:276178/1‑131 (MQ=255)                                                     gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTg                                                                                                >  1:44267/1‑149 (MQ=255)                                                     gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTg                                                                                                >  3:135195/1‑149 (MQ=255)                                                      cGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGc                                                                                               >  3:151564/1‑149 (MQ=255)                                                      cGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCGTGAGGTCGTTGCGTTTGGCCAGCAGGCGCCCCTGc                                                                                               >  6:13407/1‑149 (MQ=255)                                                       gACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCt                                                                                              >  7:298305/1‑149 (MQ=255)                                                        aCCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTc                                                                                             <  2:74120/149‑1 (MQ=255)                                                                   cTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGt                                                                                  >  7:173553/1‑149 (MQ=255)                                                                   cTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGt                                                                                  >  7:240837/1‑149 (MQ=255)                                                                    tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑g                                                                                >  2:36767/1‑149 (MQ=255)                                                                    tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑g                                                                                >  1:79683/1‑149 (MQ=255)                                                                    tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTGGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑g                                                                                >  3:26796/1‑149 (MQ=255)                                                                      gCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGc                                                                              >  8:188593/1‑149 (MQ=255)                                                                            cTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGcc                                                                                                              <  5:255476/112‑1 (MQ=255)                                                                            cTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGcc                                                                                                              >  6:255476/1‑112 (MQ=255)                                                                               ggTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                                   >  6:114797/1‑120 (MQ=255)                                                                                 tCACCCCGATGACGTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAAc                                                                   <  2:199435/149‑1 (MQ=255)                                                                                 cccccccCATCACCTGGTTGCCATCCTTTGGGCCACCCATCTCCATTCGGGGCGTGCCCTTGGGCTTC‑CCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAAc                                                                   <  7:204454/146‑1 (MQ=255)                                                                                     cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTgc                                                               >  3:310172/1‑149 (MQ=255)                                                                                     cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTgc                                                               >  5:225969/1‑149 (MQ=255)                                                                                       cGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTgccg                                                             >  4:202876/1‑149 (MQ=255)                                                                                              cTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCa                                                      >  3:169434/1‑149 (MQ=255)                                                                                              cTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGCCCCCGCCAGGCCCTAGAGGTCGTTGCGTTTGGCCACCAGGAGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCTCCGATCa                                                      >  4:236589/1‑149 (MQ=255)                                                                                               tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTcgac                                                                       >  6:77245/1‑131 (MQ=255)                                                                                               tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTcgac                                                                       <  5:77245/131‑1 (MQ=255)                                                                                               tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCAc                                                     >  5:184776/1‑149 (MQ=255)                                                                                                 gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCAccc                                                   >  5:134300/1‑149 (MQ=255)                                                                                                 gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCAccc                                                   <  7:188593/149‑1 (MQ=255)                                                                                                 gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCAccc                                                   >  6:82170/1‑149 (MQ=255)                                                                                                     gCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGgt                                               >  3:293027/1‑149 (MQ=255)                                                                                                       aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGgtgt                                             >  5:265647/1‑149 (MQ=255)                                                                                                       aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGgtgt                                             >  6:146856/1‑149 (MQ=255)                                                                                                       aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGgcgt                                             >  8:136954/1‑149 (MQ=255)                                                                                                          cTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGAACACCCCGGTGTgg                                           <  2:44267/148‑1 (MQ=255)                                                                                                             gtgGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                                        >  2:255762/1‑148 (MQ=255)                                                                                                                gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                                        <  5:150283/145‑1 (MQ=255)                                                                                                                gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                                        >  6:150283/1‑145 (MQ=255)                                                                                                                 tCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                                    >  5:73496/1‑148 (MQ=255)                                                                                                                   aCCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGc                                   >  3:116849/1‑147 (MQ=255)                                                                                                                     ccATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                               >  8:178666/1‑149 (MQ=255)                                                                                                                     ccATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                               >  7:109391/1‑149 (MQ=255)                                                                                                                       aTCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                               >  6:114368/1‑147 (MQ=255)                                                                                                                       aTCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                               <  5:114368/147‑1 (MQ=255)                                                                                                                          tcCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg                          >  2:36305/1‑149 (MQ=255)                                                                                                                           ccctGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTCGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg                          <  1:255762/146‑1 (MQ=255)                                                                                                                            cATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCCCCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCACGGCGTGgcg                        >  1:289168/1‑149 (MQ=255)                                                                                                                            cATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                                           <  1:241650/130‑1 (MQ=255)                                                                                                                            cATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                                           >  2:241650/1‑130 (MQ=255)                                                                                                                            cATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg                          <  7:178666/147‑1 (MQ=255)                                                                                                                            cATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcg                        >  1:290184/1‑149 (MQ=255)                                                                                                                             aTGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc                       >  5:205248/1‑149 (MQ=255)                                                                                                                               gCCGGGCCGCCCCATGGGGCTGCCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGcc                      <  4:151564/149‑1 (MQ=255)                                                                                                                               gCCGGGCATGCCCATGGGCTCG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTcgac                                                                       <  7:237294/99‑1 (MQ=38)                                                                                                                               gCCGGGCATGCCCATGGGCTCG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTcgac                                                                       >  8:237294/1‑99 (MQ=38)                                                                                                                                ccGGGCATGCCCATGGGCTCG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCAc                    >  1:297764/1‑149 (MQ=255)                                                                                                                                 cGGGCGCGCCCTGGGGCTGG‑CCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCAca                   <  4:310172/149‑2 (MQ=255)                                                                                                                                  gggCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc                  >  6:211327/1‑149 (MQ=255)                                                                                                                                    gCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                                  >  7:77256/1‑52 (MQ=255)                                                                                                                                    gCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                                  <  8:77256/52‑1 (MQ=255)                                                                                                                                    gCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCg                <  7:136954/149‑1 (MQ=38)                                                                                                                                     gcaGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGGCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACTAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGa               <  8:52073/147‑1 (MQ=255)                                                                                                                                     cGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCgg               <  4:169434/149‑2 (MQ=38)                                                                                                                                        gCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTACCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGc            >  8:300387/1‑149 (MQ=255)                                                                                                                                         cccATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGTACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGACACGCCGAGGc            >  1:113102/1‑148 (MQ=255)                                                                                                                                         cccATGGGAATG‑CCCGGAACCGCCAGGCCCTTGCGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCGGGTAGGT‑GGCCTTCGCCGAACTTCCCGCCTATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGc            <  2:113102/148‑1 (MQ=255)                                                                                                                                          ccATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt          >  4:167776/1‑149 (MQ=255)                                                                                                                                          ccATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt          >  3:275463/1‑149 (MQ=255)                                                                                                                                          ccATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt          >  1:71628/1‑149 (MQ=255)                                                                                                                                            aTGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTg        >  2:308884/1‑149 (MQ=255)                                                                                                                                               ggCATC‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGGTTGGCCAGCAGGCGCACCGGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGGGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGc       <  4:195240/147‑1 (MQ=255)                                                                                                                                                gCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCtt     >  7:231661/1‑148 (MQ=255)                                                                                                                                                 cATG‑CCCGGCACCGCCCGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGa   <  1:36305/149‑1 (MQ=255)                                                                                                                                                  aTG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGAc  >  1:207387/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATG‑CCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGT‑GGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGAC  >  NC_002947/2508644‑2508937

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