Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,122,761	G→A	9.2%	intergenic (+32/+115)	PP_0983 → / ← PP_0984	membrane protein/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,122,761	0	G	A	9.2%	179.9 / 6.7	65	intergenic (+32/+115)	PP_0983/PP_0984	membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base A (3/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122613‑1122905                                                                                                                                                     |                                                                                                                                                  cAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTACGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAg                                                                                                                                                   >  8:91754/1‑149 (MQ=255)  aGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAgg                                                                                                                                                  <  7:36539/149‑1 (MQ=255)    gACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGt                                                                                                                                                >  6:112553/1‑149 (MQ=255)          cTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCt                                                                                                                                          <  8:119864/148‑1 (MQ=255)              ggCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTC‑GGCtttct                                                                                                                                      <  8:290920/148‑3 (MQ=255)                 cGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCTGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTC‑GGCTTTCTtgtg                                                                                                                                  <  1:246754/149‑1 (MQ=255)                        ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGtt                                                                                                                            >  2:113723/1‑148 (MQ=255)                        ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGtt                                                                                                                            >  5:109395/1‑148 (MQ=255)                        ccAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGtt                                                                                                                            >  4:134298/1‑148 (MQ=255)                           ggTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATc                                                                                                                        >  8:124702/1‑149 (MQ=255)                            gTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCt                                                                                                                       <  8:42755/149‑1 (MQ=255)                            gTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCt                                                                                                                       <  3:134298/149‑1 (MQ=255)                                 cGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGccc                                                                                                                  >  1:209050/1‑149 (MQ=255)                                       cccGCCACTGCTGGCTGTGGTGATACCAGCGGGCATGTGTGCCTTGGCAGGTGTGTGGTTGGTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGTCCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCct                                                                                                             <  8:53680/148‑1 (MQ=255)                                        ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCctct                                                                                                           <  4:162928/149‑1 (MQ=255)                                        ccGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCctct                                                                                                           >  5:219191/1‑149 (MQ=255)                                          gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCt                                                                                                                                          >  4:262817/1‑116 (MQ=255)                                          gCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCt                                                                                                                                          <  3:262817/116‑1 (MQ=255)                                          gCCACTGCTGGCGGTGGTGATACCAGCAGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTc                                                                                                         >  6:229363/1‑149 (MQ=255)                                              ctgctgGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCgg                                                                                                     <  5:112553/149‑1 (MQ=255)                                                       gtggtgATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcg                                                                                            >  3:141321/1‑149 (MQ=255)                                                            gATACCAGGGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTC‑GGCtttcttgt                                                                                                                                   <  4:112874/105‑6 (MQ=38)                                                            gATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCCTGAGGTC‑GGCtttcttgt                                                                                                                                   >  3:112874/1‑100 (MQ=38)                                                               aCCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGt                                                                                                                                   >  4:296610/1‑102 (MQ=255)                                                               aCCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGt                                                                                                                                   <  3:296610/102‑1 (MQ=255)                                                                  aGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAgg                                                                                  >  4:46621/1‑148 (MQ=255)                                                                             gtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAgg                                                                       >  3:218576/1‑148 (MQ=255)                                                                             gtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAgg                                                                       >  1:108962/1‑148 (MQ=255)                                                                              tgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         <  1:113723/145‑1 (MQ=255)                                                                              tgCCTTGGCCGGTGTGTGGTTGTTACGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAg                                                                        <  7:91754/146‑1 (MQ=255)                                                                                ccTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATc                                                                   <  2:209050/149‑1 (MQ=255)                                                                                  ttGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAg                                                                 >  3:232477/1‑149 (MQ=255)                                                                                    ggCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAgg                                                                <  5:229363/148‑1 (MQ=255)                                                                                     gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         >  7:20556/1‑138 (MQ=255)                                                                                     gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCa                                                                         <  8:20556/138‑1 (MQ=255)                                                                                     gCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGa                                                               <  7:124702/148‑1 (MQ=255)                                                                                              tgGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCg                                                      >  2:180795/1‑148 (MQ=255)                                                                                                    ttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCcgc                                                                                           <  1:138969/105‑1 (MQ=255)                                                                                                    ttgCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGACCGGCCTCTTCGCGGGTAAGCCcgc                                                                                           >  2:138969/1‑105 (MQ=255)                                                                                                       cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACt                                             <  2:108962/148‑1 (MQ=255)                                                                                                       cgcCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            >  5:198643/1‑149 (MQ=255)                                                                                                          cGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCTTCAGGAAAGGGCCCGAAGGCAACTCCaa                                         <  1:180795/149‑1 (MQ=255)                                                                                                            ggCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATc                                       <  4:218576/149‑1 (MQ=255)                                                                                                               aGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                  <  2:243271/119‑1 (MQ=255)                                                                                                               aGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCa                                                                  >  1:243271/1‑119 (MQ=255)                                                                                                                 ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAgg                                                           >  5:205824/1‑124 (MQ=255)                                                                                                                 ggTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAgg                                                           <  6:205824/124‑1 (MQ=255)                                                                                                                  gTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTgg                                 <  4:232477/149‑1 (MQ=255)                                                                                                                     aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            <  3:291129/135‑1 (MQ=255)                                                                                                                     aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            >  4:291129/1‑135 (MQ=255)                                                                                                                       aGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAt                            >  1:63531/1‑149 (MQ=255)                                                                                                                         ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTCAGGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAtt                           >  1:202186/1‑149 (MQ=255)                                                                                                                         ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAttt                          >  8:86178/1‑149 (MQ=255)                                                                                                                         ccTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTAtt                           >  7:204821/1‑148 (MQ=255)                                                                                                                           tGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttct                        >  1:83825/1‑149 (MQ=255)                                                                                                                             aGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttcttc                      <  3:46621/149‑1 (MQ=255)                                                                                                                                  cacgaaaaaGCCGACCCTGAGGTC‑GGCTTTCTTGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTATAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCtt                 >  6:160389/7‑149 (MQ=38)                                                                                                                                  cacaAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGTAAAGGGCCCGAAGGCACCTCCAATCCTCTGGCTTATTTCTTCTGCtt                 >  2:226479/1‑149 (MQ=255)                                                                                                                                     aaGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg              <  2:202186/149‑1 (MQ=255)                                                                                                                                     aaGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg              >  3:157624/1‑149 (MQ=255)                                                                                                                                       gAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCAGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCa            >  7:5886/1‑149 (MQ=255)                                                                                                                                       aaaaGCCGACCCTGAGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCa            >  4:202621/2‑149 (MQ=255)                                                                                                                                        aaaGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCAc           <  2:83825/149‑1 (MQ=255)                                                                                                                                           gCCGGCCTCAGGGTC‑GGCTTTTTCGTGGGAGTTGATCTGGCCCGGCCTCTTCGCGCGTAACCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGCGCCCGAAGGCACCTCCAATCCTCTGGCGAATTTCTTATGCTTCGGGACCCg        >  7:120158/1‑149 (MQ=255)                                                                                                                                                 cTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                       <  1:68360/128‑1 (MQ=255)                                                                                                                                                 cTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATttctt                       >  2:68360/1‑128 (MQ=255)                                                                                                                                                 cTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAg  <  8:204821/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122613‑1122905

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