Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,103,972	C→T	5.5%	Q269* (CAA→TAA) ‡	PP_5353 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,103,972	0	C	T	5.5%	269.1 / 2.3	89	Q269* (CAA→TAA) ‡	PP_5353	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base T (3/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

TACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103825‑6104120                                                                                                                                                      |                                                                                                                                                     tacATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                      >  3:241805/1‑148 (MQ=255) tacATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                      >  6:310025/1‑148 (MQ=255)  acATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                      >  3:102833/1‑147 (MQ=255)    aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCaatg                                                                                                                                                  <  2:250075/149‑3 (MQ=255)    aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                                  >  5:249035/1‑149 (MQ=255)     tCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATAACGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                                  <  1:298302/149‑1 (MQ=255)       gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCga                                                                                                                                               <  4:59204/149‑1 (MQ=255)       gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGCTCGCGCCCAAGGCga                                                                                                                                               <  7:4900/149‑1 (MQ=255)            ccTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  4:59887/1‑148 (MQ=255)             cTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                           >  5:243608/1‑147 (MQ=255)                    tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  8:336343/149‑8 (MQ=255)                    tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTgg                                                                                                                                  <  4:241805/149‑1 (MQ=255)                     gCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                                  <  4:307558/148‑8 (MQ=255)                      cAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  4:354260/148‑1 (MQ=255)                      cAAATGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTTTCGAAGATT‑TCGCGCTTCGACGAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATTGTTCGCGGCCAAGGCGAGAACCAGCCTTGGc                                                                                                                                 <  3:122498/148‑1 (MQ=255)                       aaCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTGCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTg                                                                                                                               <  3:34012/149‑1 (MQ=255)                        aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  5:247633/1‑149 (MQ=255)                        aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  6:249035/149‑1 (MQ=255)                              cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCTCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCGCAAGTCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  8:149030/149‑1 (MQ=255)                              cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:181409/1‑149 (MQ=255)                               ggTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCgg                                                                                                                       <  5:251697/149‑1 (MQ=255)                                   tCGATGACCGCGCCGAGCGGACCGACGCTACCCTGACCTTCCGGGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCTAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCTCCCAAGGCGAGAACCAGCGTTGGCTGGTAGCCGGTAt                                                                                                                    <  4:324506/148‑1 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              <  2:301422/136‑1 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                              >  1:301422/1‑136 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCg                                                                                                                 >  4:185334/1‑149 (MQ=255)                                          ccGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCACGGCGAA‑GTGTTCAGCGAAACCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCAGTCaa                                                                                                             >  6:174235/1‑148 (MQ=255)                                              gcCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAg                                                                                                        <  3:132764/149‑1 (MQ=255)                                              gcCCACCGCACCTAACGTACCCTGTCCTTACGTGGGGTTTCGAATAAT‑TGGGGGTTTGCTCATGGGGAATGTCTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACa                                                                                                         <  1:117811/149‑1 (MQ=255)                                                      cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  3:296732/1‑149 (MQ=255)                                                      cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                                >  7:127943/1‑149 (MQ=255)                                                         cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:133575/1‑149 (MQ=255)                                                         cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             >  5:328896/1‑149 (MQ=255)                                                            cGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTg                                                                                          >  5:70200/1‑149 (MQ=255)                                                                 cccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  2:346478/148‑1 (MQ=255)                                                                 cccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                      <  4:287399/148‑1 (MQ=255)                                                                  ccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  2:166267/1‑148 (MQ=255)                                                                  ccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                     >  7:345296/1‑148 (MQ=255)                                                                    tGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCGTTGTCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                             <  3:351565/138‑1 (MQ=255)                                                                    tGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACAAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCgacaacgg                                                                                     >  4:351565/1‑139 (MQ=255)                                                                       ccTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  7:346079/1‑148 (MQ=255)                                                                       ccTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                                >  5:342965/1‑148 (MQ=255)                                                                         ttCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  4:276051/1‑149 (MQ=255)                                                                         ttCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                             >  3:63542/1‑149 (MQ=255)                                                                            cGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAcccc                                                                          <  6:247633/149‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAATCAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  6:243608/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  8:132991/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                       <  2:254878/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCACGGGc                                                                      >  6:264426/1‑149 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        >  3:95013/1‑94 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                        <  4:95013/94‑1 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGcc                                                                 >  1:259659/1‑149 (MQ=255)                                                                                       gaagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCg                                                               <  4:102833/149‑1 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:219508/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:190046/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCAGGTATCCGTCAAAACAGCTAATCGCAGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                              >  6:346944/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCACCGAAAGCTGGCACATGGTAAGCGCCCAAGGGGAGAACCAGCCTTGGCTGGTAGCCGCTATCCGTCAAAACAGCTAATCGCTGCCTGCTAAACACAGACCCACGGGCTTGCCCgg                                                              >  4:31648/1‑149 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  6:342965/149‑1 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             <  8:346079/149‑1 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                             >  1:144031/1‑149 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCggg                                                             <  4:296732/149‑1 (MQ=255)                                                                                          gaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                            >  2:95508/1‑149 (MQ=255)                                                                                           aaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGt                                                           <  6:70200/149‑1 (MQ=255)                                                                                                 gcgcTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                      <  1:145511/149‑5 (MQ=255)                                                                                                    cTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                                   >  5:161083/1‑149 (MQ=255)                                                                                                       cGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttt                                                  <  5:190046/147‑1 (MQ=255)                                                                                                        gACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttttg                                                <  4:181409/148‑1 (MQ=255)                                                                                                                cGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTATTCTTTTGCCAAAGTgg                                       >  7:7162/1‑149 (MQ=255)                                                                                                                  aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    >  5:94509/1‑102 (MQ=255)                                                                                                                  aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                                    <  6:94509/102‑1 (MQ=255)                                                                                                                       gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  7:17414/1‑148 (MQ=255)                                                                                                                       gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  5:213418/1‑148 (MQ=255)                                                                                                                        ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                                  >  6:14685/1‑147 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                             <  3:276051/149‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  5:219508/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  8:345296/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  6:180329/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  2:259659/148‑1 (MQ=255)                                                                                                                               aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              >  2:34280/1‑144 (MQ=255)                                                                                                                               aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGCTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                              <  1:34280/144‑1 (MQ=255)                                                                                                                                 aGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAg                       <  5:310025/149‑1 (MQ=255)                                                                                                                                  gCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAgg                      <  3:237052/149‑1 (MQ=255)                                                                                                                                   cTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                     >  8:160733/1‑149 (MQ=255)                                                                                                                                      aaggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcg                    >  2:372894/6‑147 (MQ=255)                                                                                                                                       aggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACACCGGGCTTGGCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAg                 >  5:242158/5‑149 (MQ=255)                                                                                                                                        acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                <  8:17414/149‑1 (MQ=255)                                                                                                                                        acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc                >  2:127263/1‑149 (MQ=255)                                                                                                                                           tGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt             >  8:99709/1‑149 (MQ=255)                                                                                                                                             gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  6:71724/1‑149 (MQ=255)                                                                                                                                             gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           <  5:224881/149‑1 (MQ=255)                                                                                                                                              ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt           >  7:128062/1‑148 (MQ=255)                                                                                                                                               tCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCTGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTc          <  6:328896/148‑1 (MQ=255)                                                                                                                                                 gcgcCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATc       <  6:213418/149‑1 (MQ=255)                                                                                                                                                    cccAATGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg    <  3:59887/149‑1 (MQ=255)                                                                                                                                                     ccAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATACGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCAGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAgg   >  5:298397/1‑149 (MQ=255)                                                                                                                                                      cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  8:128062/149‑1 (MQ=255)                                                                                                                                                      cAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGt  <  1:166267/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGT  >  NC_002947/6103825‑6104120

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: ‑