Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,633,397	A→T	6.3%	F572I (TTC→ATC)	uvrC ←	UvrABC system protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,633,397	0	A	T	6.3%	183.1 / 2.9	64	F572I (TTC→ATC)	uvrC	UvrABC system protein C
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/2);  total (32/32)
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.45e-01

CGGTACGACTGGCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAGGT  >  NC_002947/4633254‑4633540                                                                                                                                                |                                                                                                                                                 cGGTACGACTGGCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGt                                                                                                                                              <  3:252434/148‑1 (MQ=255)       gACTGGCTGCAAGTTGAGCCC‑GCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGt                                                                                                                                              >  8:73027/1‑141 (MQ=255)       gACTGGCTGCAAGTTGAGCCC‑GCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGt                                                                                                                                              <  7:73027/141‑1 (MQ=255)         cTGGCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAAt                                                                                                                                                <  4:300718/138‑1 (MQ=255)         cTGGCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAAt                                                                                                                                                >  3:300718/1‑138 (MQ=255)            gCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGt                                                                                                                                    <  1:281083/147‑1 (MQ=255)            gCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTcg                                                                                                                                  >  1:5038/1‑149 (MQ=255)            gCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTcg                                                                                                                                  >  3:275467/1‑149 (MQ=255)             cTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGTCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTcgc                                                                                                                                 >  2:5074/1‑149 (MQ=255)              tccAACTTGAGCCAGGAATTCAACTCCCTATGCAGGCTGGCATAAATCGACTCGCCAAGCTTTTTAATGATCCCGGGGGCATAGGCGATTTCATCGATACTGGCGCGGTTGAGCTCATGCAGGACGCCGAAATGTTTAAGCACGTcgcg                                                                                                                                >  4:320195/3‑149 (MQ=255)              tGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTcgcg                                                                                                                                >  6:244807/1‑149 (MQ=255)                   gTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCgcggcggc                                                                                                                           >  5:75544/1‑149 (MQ=255)                   gTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCgcggcggc                                                                                                                           <  5:214516/149‑1 (MQ=255)                    ttGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCgcggcggcg                                                                                                                          >  3:302889/1‑149 (MQ=255)                          cctgCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTTCTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTgg                                                                                                                     <  3:320195/145‑1 (MQ=255)                          ccGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTggg                                                                                                                    <  8:354348/149‑1 (MQ=255)                          ccGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTggg                                                                                                                    >  7:354348/1‑149 (MQ=255)                                  cTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGaaa                                                                                                                                                 <  4:78138/112‑1 (MQ=255)                                  cTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGaaa                                                                                                                                                 >  3:78138/1‑112 (MQ=255)                                   tACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTAcccc                                                                                                           >  2:84835/1‑149 (MQ=255)                                      tCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGc                                                                                                        >  7:50073/1‑149 (MQ=255)                                       cGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGc                                                                                                        >  7:291108/1‑148 (MQ=255)                                            tGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACAt                                                                                                   <  4:86684/148‑1 (MQ=255)                                               aGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTc                                                                                               <  6:103898/149‑1 (MQ=255)                                                     gcgtAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCt                                                                                         >  8:171292/4‑149 (MQ=255)                                                       aTAAATCGACTCGGCAATCTTTTTACTGATTCCGTGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGATCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCACGGCCACATCTTAAAGGCTgg                                                                                       <  2:78278/149‑1 (MQ=255)                                                       aTAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGAGATTTCATCGATACTGGCGCTGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCTCGGCGGCGCTTGGGACCTACACCGTCCCCATCTTCAAGGCTg                                                                                        >  1:78278/1‑148 (MQ=255)                                                        tAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGga                                                                                      <  8:291108/149‑1 (MQ=255)                                                           aTCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGAc                                                                                   >  2:36349/1‑149 (MQ=255)                                                           aTCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGAc                                                                                   >  1:59314/1‑149 (MQ=255)                                                           aTCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGAc                                                                                   <  5:201724/149‑1 (MQ=255)                                                           aTCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGAc                                                                                   >  6:201724/1‑149 (MQ=255)                                                              gACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTc                                                                                 <  5:244807/148‑1 (MQ=255)                                                                   ggCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAg                                                                           >  3:179708/1‑149 (MQ=255)                                                                   ggCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAg                                                                           >  5:161297/1‑149 (MQ=255)                                                                    gCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGc                                                                          <  2:5038/149‑1 (MQ=255)                                                                          tttttACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGATATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑c                                                                    >  3:239916/1‑148 (MQ=255)                                                                          tttttACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGATATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑c                                                                    >  3:244568/1‑148 (MQ=255)                                                                            tttACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑cc                                                                   <  6:75544/147‑1 (MQ=255)                                                                                 tGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACgg                                                            >  8:373/1‑149 (MQ=255)                                                                                      ccGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCg                                                       <  4:302889/149‑1 (MQ=255)                                                                                               ttGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGt                                              >  6:140604/1‑149 (MQ=255)                                                                                                 ggCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGc                                                                                          <  2:220116/104‑1 (MQ=255)                                                                                                 ggCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGc                                                                                          >  1:220116/1‑104 (MQ=255)                                                                                                    gATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAg                                         >  1:209379/1‑149 (MQ=255)                                                                                                     aTTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTGCCCACGACGGGCCCGGTGGCCGGTGATAg                                         <  1:5074/149‑1 (MQ=255)                                                                                                     aTTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGc                                        <  6:161297/149‑1 (MQ=255)                                                                                                        tCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCg                                       >  1:145349/1‑147 (MQ=255)                                                                                                         cATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGaa                                     <  8:50073/148‑1 (MQ=255)                                                                                                          aTCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCgg                                               <  8:133715/137‑1 (MQ=255)                                                                                                           tCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGCCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCg                                  <  7:171292/149‑1 (MQ=255)                                                                                                             gATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGt                                >  6:162214/1‑149 (MQ=255)                                                                                                              aTACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTg                               <  2:209379/149‑1 (MQ=255)                                                                                                                  tGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCt                           >  5:298347/1‑149 (MQ=255)                                                                                                                    gcgcGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCtt                          <  4:179708/148‑1 (MQ=255)                                                                                                                         gTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTcgc                    <  7:373/149‑1 (MQ=255)                                                                                                                          ttGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTcgcg                   <  4:275467/149‑1 (MQ=255)                                                                                                                                tCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGc                                        >  5:167087/1‑122 (MQ=255)                                                                                                                                tCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGc                                        <  6:167087/122‑1 (MQ=255)                                                                                                                                        gccgccGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTcgc                    <  1:89902/134‑1 (MQ=255)                                                                                                                                        gccgccGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTcgc                    >  2:89902/1‑134 (MQ=255)                                                                                                                                         ccgccgAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAg    <  1:267028/149‑1 (MQ=255)                                                                                                                                           gccgATATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCTCATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAGGt  <  4:239916/149‑1 (MQ=255)                                                                                                                                           gccgATATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAGGt  <  4:244568/149‑1 (MQ=255)                                                                                                                                           gccgAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAgg   >  2:259751/1‑148 (MQ=255)                                                                                                                                           gccgAAATGTTTAAGCAGGGCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAgg   <  1:259751/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CGGTACGACTGGCTGCAAGTTGAGCCCGGCATTCTACTCGCTATGCAGGCTGGCATAAATCGACTCGGCAAGCTTTTTACTGATGCCGGGGGCTTTGGCGATTTCATCGATACTGGCGCGGTTGAGCTCCTGCAGGCCGCCGAAATGTTTAAGCAGGTCGCGGCGGCGCTTGGGCCCTACCCCGGCCACATCTTCAAGGCTGGAGACTCGACGAGCCTTG‑CCACGACGGGCCCGGTGGCCGGTGATAGCGAAGCGGTGGGCTTCGTCGCGGATCTGCTGGATCAGGT  >  NC_002947/4633254‑4633540

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