Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,040,747	C→G	6.5%	F41L (TTC→TTG)	pyrF →	orotidine 5'‑phosphate decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,040,747	0	C	G	6.5%	192.7 / 4.1	62	F41L (TTC→TTG)	pyrF	orotidine 5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CGACTTGCCAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCGTGTGGATGGTCAACGTGCA  >  NC_002947/2040602‑2040893                                                                                                                                                  |                                                                                                                                                   cGACTTGCCAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCAc                                                                                                                                                  >  6:21987/1‑148 (MQ=255)  gACTTGACAAGAGACCCTGCCCATGTCCTCCTGCCAGCCGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGTCTGACCAGCTTGACCCTGCCCTGTGCCGGGTGAAGGTTGGCAAGGATCTGTTCACca                                                                                                                                                <  6:212938/149‑1 (MQ=255)     ttGCCAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACcag                                                                                                                                               >  8:284599/1‑147 (MQ=255)     ttGCCAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACACGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACcag                                                                                                                                               >  4:40696/1‑147 (MQ=255)        ccAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAgcg                                                                                                                                          >  3:25008/1‑149 (MQ=255)            gagaCCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTc                                                                                                                                      >  5:149137/1‑149 (MQ=255)               aCCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACcag                                                                                                                                               >  6:322352/1‑137 (MQ=255)               aCCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACcag                                                                                                                                               <  5:322352/137‑1 (MQ=255)                 ccTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCa                                                                                                                                 >  5:292064/1‑149 (MQ=255)                      ccATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTTGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTc                                                                                                                            <  3:40696/149‑1 (MQ=255)                      ccATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTc                                                                                                                            >  3:286205/1‑149 (MQ=255)                           tCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGaaa                                                                                                                        <  4:25008/148‑1 (MQ=255)                           tCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAAc                                                                                                                       <  5:21987/149‑1 (MQ=255)                                 tGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCtgtg                                                                                                                 >  1:158701/1‑149 (MQ=255)                                       aCGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCTGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACaa                                                                                                           <  5:365274/149‑1 (MQ=255)                                         gCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAgg                                                                                                         <  6:283033/149‑1 (MQ=255)                                           cccTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGc                                                                                                       >  7:224009/1‑149 (MQ=255)                                            ccTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCt                                                                                                      >  1:333750/1‑149 (MQ=255)                                                 tcgtcgCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGaa                                                                                                 <  2:333750/149‑1 (MQ=255)                                                         cTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTGACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         >  4:134810/1‑149 (MQ=255)                                                         cTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTGACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         >  4:140466/1‑149 (MQ=255)                                                         cTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         >  3:357851/1‑149 (MQ=255)                                                         cTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         >  2:35394/1‑149 (MQ=255)                                                         cTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         <  7:284599/149‑1 (MQ=255)                                                             aTTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGAc                                                                                     >  5:147538/1‑149 (MQ=255)                                                             aTTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGAc                                                                                     >  3:311191/1‑149 (MQ=255)                                                               ttCCCTACCCGTGAGGCCGCCCAGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTAGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCt                                                                                   >  7:55429/1‑149 (MQ=255)                                                                  ccTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCAtt                                                                                                                               >  5:363240/1‑102 (MQ=255)                                                                  ccTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCAtt                                                                                                                               <  6:363240/102‑1 (MQ=255)                                                                   cTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCg                                                                                                                                     <  7:30265/95‑1 (MQ=255)                                                                   cTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCg                                                                                                                                     >  8:30265/1‑95 (MQ=255)                                                                   cTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAg                                                                               >  5:74790/1‑149 (MQ=255)                                                                   cTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAg                                                                               >  1:240956/1‑149 (MQ=255)                                                                       ccGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTc                                                                            >  4:50762/1‑148 (MQ=255)                                                                        cGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTGGAAGTGTTCCTGGACCTCAAGTTCCa                                                                          <  7:272207/149‑1 (MQ=255)                                                                        cGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCa                                                                          <  4:318008/149‑1 (MQ=255)                                                                         gTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCAc                                                                         <  6:149137/149‑1 (MQ=255)                                                                          tGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACg                                                                        <  8:224009/149‑1 (MQ=255)                                                                             ggCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACa                                                                     <  6:135178/149‑1 (MQ=255)                                                                                     tGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCaa                                                              >  8:131706/1‑148 (MQ=255)                                                                                                  ccAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATgg                                                <  4:286205/149‑1 (MQ=255)                                                                                                       ttGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         <  5:318128/103‑1 (MQ=255)                                                                                                       ttGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCt                                                                                         >  6:318128/1‑103 (MQ=255)                                                                                                          aCCCTGCGCTGTGCCGGGTGAAGGTTGTCAAGGTGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAAGAGCTCCGCCATGGCGGTAAAg                                        <  8:267664/149‑1 (MQ=255)                                                                                                          aCCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGCGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGtggtg                                        <  1:35394/149‑5 (MQ=255)                                                                                                           cccTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAgg                                       <  6:147538/149‑1 (MQ=255)                                                                                                             cTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCt                                     <  3:50762/149‑1 (MQ=255)                                                                                                                cgcTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAacac                                                           >  3:234554/1‑124 (MQ=255)                                                                                                                cgcTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAacac                                                           >  3:239262/1‑124 (MQ=255)                                                                                                                cgcTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAacac                                                           <  4:239262/124‑1 (MQ=255)                                                                                                                cgcTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAacac                                                           <  4:234554/124‑1 (MQ=255)                                                                                                                cgcTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTgc                                   <  6:74790/148‑1 (MQ=255)                                                                                                                           gTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATggg                       >  8:34114/1‑149 (MQ=255)                                                                                                                             gAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACAtt                                                                   >  4:187594/1‑103 (MQ=255)                                                                                                                             gAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACAtt                                                                   <  3:187594/103‑1 (MQ=255)                                                                                                                             gAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCg                     >  2:22858/1‑149 (MQ=255)                                                                                                                              aaGGTTGGCAAGGAGCTGTTGACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCgt                    <  3:140466/149‑1 (MQ=255)                                                                                                                               aGGTTGGCAAGGAGCTGTTGACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCgt                    <  3:134810/148‑1 (MQ=255)                                                                                                                                     gCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCGTGTGGAt              <  7:131706/148‑1 (MQ=255)                                                                                                                                     gCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCGTGTGGAt              <  4:311191/148‑1 (MQ=255)                                                                                                                                     gCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCGTGTGGAt              <  1:22858/148‑1 (MQ=255)                                                                                                                                                ttCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAACGGCTGCCGCCGAGATGGGCGTGTGGGTGGGCAACGTGCa  >  4:121986/1‑149 (MQ=255)                                                                                                                                                  cACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCg                                              >  6:110989/1‑103 (MQ=255)                                                                                                                                                  cACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCg                                              <  5:110989/103‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGACTTGCCAAGAGACCCTGCCCATGTCCGCCTGCCAGACGCCCCTGATCGTCGCCCTGGATTTCCCTACCCGTGAGGCCGCCCTGAAGCTGGCTGACCAGCTTGACCCTGCGCTGTGCCGGGTGAAGGTTGGCAAGGAGCTGTTCACCAGCAGCGCTTCGGGCATTGTCGAAACCCTGTGTGACAAGGGCTTCGAAGTGTTCCTGGACCTCAAGTTCCACGACATTCCCAACACCACCGCCATGGCGGTAAAGGCTGCCGCCGAGATGGGCGTGTGGATGGTCAACGTGCA  >  NC_002947/2040602‑2040893

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