Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,765,101	G→T	6.1%	intergenic (+1621/+54)	PP_3326 → / ← PP_3328	hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,765,101	0	G	T	6.1%	209.4 / 2.4	65	intergenic (+1621/+54)	PP_3326/PP_3328	hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base T (2/2);  total (32/33)
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.06e-01

ACAAGACCTGCAATAATTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCGCTG  >  NC_002947/3764953‑3765241                                                                                                                                                     |                                                                                                                                             acaAGACCTGCAATAATTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACg                                                                                                                                              >  1:26873/1‑149 (MQ=255) acaAGACCTGCAATAATTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACg                                                                                                                                              >  3:1701/1‑149 (MQ=255)   aaGACCTGCAATAATTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGcc                                                                                                                                            <  3:166899/149‑1 (MQ=255)            aataatTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAtct                                                                                                                                    >  3:85296/1‑148 (MQ=255)                atTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATAGTTGGAGCATGTTATTTTTATTGAGTGGCCACTCCCCGATACACGCCTCTTATCTCACg                                                                                                                                >  3:252396/1‑148 (MQ=255)                    tAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTAACCGATACACGCCTCTTATCTCACGCAGa                                                                                                                            >  3:219439/1‑148 (MQ=255)                      gTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGa                                                                                                                            <  2:26873/146‑1 (MQ=255)                          ccAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCg                                                                                                                     <  1:258436/149‑1 (MQ=255)                           cAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCgg                                                                                                                    >  7:305510/1‑149 (MQ=255)                             gACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGa                                                                                                                  <  4:219439/149‑1 (MQ=255)                              aCTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAtctc                                                                                                                                   >  6:39291/1‑131 (MQ=255)                              aCTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAtctc                                                                                                                                   <  5:39291/131‑1 (MQ=255)                                  gTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGACACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGgcgc                                                                                                             >  1:331438/1‑149 (MQ=255)                                   tCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCGTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGcgcg                                                                                                            >  8:285235/1‑149 (MQ=255)                                   tCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACTCCTCTTATCTCACGCAGACAAAGCGGGAGGcgcg                                                                                                            >  5:113288/1‑149 (MQ=255)                                   tCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGGTATTTTTATTGAGTGGCCACTCACCGATACACTCCTCTTATCTCACGCAGACAAAGCGGGAGGcgcg                                                                                                            >  7:313669/1‑149 (MQ=255)                                     tttGCCATACGTCCTCCAGTTATTATATTTTGGACTGTTCATTGGCCGGCAAAAATGCCAAACACATCGTT‑GAGCGTGTTATTTTTATTGAGTGGCAACTCACCGATACAGGCCTCTAATCTCACGCTGACAAAGCGGGCGGCGCGTcc                                                                                                         >  4:202014/1‑149 (MQ=255)                                     tttGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTc                                                                                                          >  3:50076/1‑149 (MQ=255)                                     tttGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTc                                                                                                          >  3:303731/1‑149 (MQ=255)                                     tttGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTc                                                                                                          >  8:152730/1‑149 (MQ=255)                                          cATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGaa                                                                                                     <  3:166794/149‑1 (MQ=255)                                           aTACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAg                                                                                                    >  7:169832/1‑149 (MQ=255)                                           aTACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAg                                                                                                    >  3:173158/1‑149 (MQ=255)                                                tcctccAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTa                                                                                               <  1:314882/149‑1 (MQ=255)                                                 cctccAGTTATGATACTTTGGACTGTACATTGGCCGGCAAAACTGCCAAACACAGCGTTGGAGCCTGTTATTTTTATTGATTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTAt                                                                                              <  4:252396/149‑1 (MQ=255)                                                  ctccAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCTCTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTAt                                                                                              <  7:333407/148‑1 (MQ=255)                                                    ccAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGc                                                                                            >  2:304504/1‑148 (MQ=255)                                                       gTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACTAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAAt                                                                                        >  1:344750/1‑149 (MQ=255)                                                       gTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAAt                                                                                        >  7:222823/1‑149 (MQ=255)                                                       gTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAAt                                                                                        <  4:50076/149‑1 (MQ=255)                                                           ttatAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAgg                                                                                    >  4:165440/1‑149 (MQ=255)                                                              taAATTGGACTGTAACTTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTg                                                                                  <  4:278029/148‑1 (MQ=255)                                                                 tttGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCa                                                                               >  5:39550/1‑148 (MQ=255)                                                                     gACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATac                                                                          <  8:222823/149‑1 (MQ=255)                                                                     gACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATa                                                                           <  7:281301/148‑1 (MQ=255)                                                                        tGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAACTCACGCAGACAAAGCg                                                                                                                     <  7:32378/103‑1 (MQ=255)                                                                        tGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAACTCACGCAGACAAAGCg                                                                                                                     >  8:32378/1‑103 (MQ=255)                                                                        tGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAACTCACGCAGACAAAGCg                                                                                                                     <  7:38166/103‑1 (MQ=255)                                                                        tGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTAACTCACGCAGACAAAGCg                                                                                                                     >  8:38166/1‑103 (MQ=255)                                                                          tACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCtt                                                                     >  8:229682/1‑149 (MQ=255)                                                                              ttGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCaaa                                                                 <  2:344750/149‑1 (MQ=255)                                                                                    ggCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCtt                                                           <  4:303731/149‑1 (MQ=255)                                                                                        aaaaTGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAAc                                                       <  8:305510/149‑1 (MQ=255)                                                                                            tGCCAAACACCTCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCtt                                                   <  4:85296/149‑1 (MQ=255)                                                                                              ccAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCttt                                                  >  5:341688/1‑148 (MQ=255)                                                                                                  acacATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCtt                                                           >  7:95005/1‑135 (MQ=255)                                                                                                  acacATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCtt                                                           <  8:95005/135‑1 (MQ=255)                                                                                                   cacaTGGTTGGAGCCTGTTATTTTTATTGAGTGGACACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGAtctc                                            <  2:331438/149‑1 (MQ=255)                                                                                                         gTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCgg                                      >  4:264612/1‑149 (MQ=255)                                                                                                          ttGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCggg                                     <  3:264612/149‑1 (MQ=255)                                                                                                                ccTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCa                               >  1:84468/1‑149 (MQ=255)                                                                                                                      aTTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTggg                          <  4:173158/148‑1 (MQ=255)                                                                                                                      aTTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTggg                          <  3:165440/148‑1 (MQ=255)                                                                                                                       tttttATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGt                        >  4:301037/1‑149 (MQ=255)                                                                                                                           tATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGTTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCa                     <  8:169832/148‑1 (MQ=255)                                                                                                                                aGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGaaa                <  7:285235/148‑1 (MQ=255)                                                                                                                                     ccACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATc          <  7:229682/149‑1 (MQ=255)                                                                                                                                        cTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGtt       <  7:152730/149‑1 (MQ=255)                                                                                                                                         tCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTc      >  2:366178/1‑149 (MQ=255)                                                                                                                                         tCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTc      >  3:152561/1‑149 (MQ=255)                                                                                                                                           aCCGATACACGCCTCTTATCTCACGCAGACACAGCGGGAGGCGCGTCCCGACGGTTTATGCCACTCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCgc    <  6:201519/149‑1 (MQ=255)                                                                                                                                           aCCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCgc    <  6:39550/149‑1 (MQ=255)                                                                                                                                            ccGATACACTCCTCTTATCTCACGCATACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCgct   <  6:113288/149‑1 (MQ=255)                                                                                                                                            ccGATACACTCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCgct   <  8:313669/149‑1 (MQ=255)                                                                                                                                             cGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCgctg  >  8:195445/1‑149 (MQ=255)                                                                                                                                                 acacGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCt                                                            <  5:63274/87‑1 (MQ=255)                                                                                                                                                 acacGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCt                                                            >  6:63274/1‑87 (MQ=255)                                                                                                                                                     |                                                                                                                                             ACAAGACCTGCAATAATTCTAGTCACCAGACTAGTCTTTGCCATACGTCCTCCAGTTATTATAATTTGGACTGTACATTGGCCGGCAAAAATGCCAAACACATCGTTGGAGCCTGTTATTTTTATTGAGTGGCCACTCACCGATACACGCCTCTTATCTCACGCAGACAAAGCGGGAGGCGCGTCCCGAAGGTTTATGCCAATCAGGTGGCAAATACACCTTCAAAGTCCTTGAACCCTTTGATCTCGATCGGGTTGCCACTGGGGTCCAGGAAAAACATCGTTCGCTG  >  NC_002947/3764953‑3765241

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