Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,122,760	A→G	8.9%	intergenic (+31/+116)	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,760	0	A	G	8.9%	183.4 / 7.8	67	intergenic (+31/+116)	PP_0983/PP_0984	membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (3/3);  total (36/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 = 5.17e-01

CCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122612‑1122905                                                                                                                                                     |                                                                                                                                                   ccAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCa                                                                                                                                                    <  4:120487/149‑1 (MQ=255)  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                                                                       >  1:108962/1‑148 (MQ=255)                                                                              gtgCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAgg                                                                       >  3:218576/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                                            >  4:291129/1‑135 (MQ=255)                                                                                                                      aTAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTc                                            <  3:291129/135‑1 (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              >  3:157624/1‑149 (MQ=255)                                                                                                                                      aaGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCgg              <  2:202186/149‑1 (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)                                                                                                                                                     |                                                                                                                                                   CCAGGACCTGCTTGGGCCGTCGAGCCAGGTATTCGGCTTCCCGCCACTGCTGGCGGTGGTGATACCAGCGGGCATCTGTGCCTTGGCCGGTGTGTGGTTGTTGCGCCGGGCAGGGTAATAGCCTGAGCTGCACAAGAAAGCCGACCTCAGGGTC‑GGCTTTTTCGTGGGTGTTGATCTGGCCCGGCCTCTTCGCGGGTAAGCCCGCGCCTACAGGTACCGCGCAGGCATCAGGAAAGGGCCCGAAGGCAACTCCAATCCTCTGGCTTATTTCTTCTGCTTCGGCACCCGCACCAG  >  NC_002947/1122612‑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: ‑