Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,389,792	C→T	6.8%	A234T (GCC→ACC)	PP_2992 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,389,792	0	C	T	6.8%	168.1 / 2.2	59	A234T (GCC→ACC)	PP_2992	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base T (2/2);  total (26/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 = 7.94e-01

GTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTTGATGAT  >  NC_002947/3389645‑3389939                                                                                                                                                    |                                                                                                                                                     gTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGATCAGCGAGCTACGTTGTCGCCAGAGAGCAACGCAAACACCGCTTCAGCAACGGCTAAGGGAGATTGATGGTTTAGGCACGAGACAAGCCTACCATACTAGATGGc                                                                                                                                                      >  7:155996/1‑148 (MQ=255)        tCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCaaa                                                                                                                                                   >  1:248812/1‑144 (MQ=255)        tCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCaaa                                                                                                                                                   <  2:248812/144‑1 (MQ=255)          ggTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTTGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGt                                                                                                                                             <  6:378712/148‑1 (MQ=255)             cTTTCTCCGCGGGTGGGAACAACTTCAGGGTTCAGCGAGCTACGTTGACGCCAGAGAGAACCGCCATAACCGCTTCACCTACGGCTTAGCGTGATTGAGCGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGAAAATGACGTGAAt                                                                                                                                         >  5:213970/1‑149 (MQ=255)                tctcCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAAt                                                                                                                                         <  7:132750/146‑1 (MQ=255)                tctcCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAAt                                                                                                                                         <  1:274282/146‑1 (MQ=255)                         gTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTg                                                                                                                                            <  3:131620/134‑1 (MQ=255)                         gTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTg                                                                                                                                            >  4:131620/1‑134 (MQ=255)                            ggAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTg                                                                                                                          <  1:111699/149‑1 (MQ=255)                             gAACAACTTGAGGGTGCAGCGAGCTACGTAGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTGTGGCCTGGGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGc                                                                                                                         <  6:57151/149‑1 (MQ=255)                              aacaacTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCg                                                                                                                        >  2:143938/1‑149 (MQ=255)                                   cTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTAcc                                                                                                                   >  1:65111/1‑149 (MQ=255)                                      gAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAg                                                                                                                >  6:71158/1‑149 (MQ=255)                                        gggTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCa                                                                                                              <  6:162521/149‑1 (MQ=255)                                                  aGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCa                                                                                                    >  7:7163/1‑149 (MQ=255)                                                      aCGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCa                                                                                                 >  7:25850/1‑148 (MQ=255)                                                       cGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCaa                                                                                                <  8:314590/148‑1 (MQ=255)                                                        gTTGTCGCCAGAGAGGACCGCGATCACCGGTTCAGCTTCGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCGTTCCATCCGTGATGGCAATTGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAAc                                                                                               <  1:243263/148‑1 (MQ=255)                                                        gTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAACACCCTGACTCACCAACt                                                                                              >  6:103429/1‑149 (MQ=255)                                                         ttGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACtt                                                                                             >  2:4578/1‑149 (MQ=255)                                                               ccAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGt                                                                                                                                             >  6:111502/1‑95 (MQ=255)                                                               ccAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGt                                                                                                                                             <  5:111502/95‑1 (MQ=255)                                                                cAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCa                                                                                      <  1:318216/149‑1 (MQ=255)                                                                 agagagCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCCGCCTTTGCGGTTCCGAGCAGCCTGACTCCCCACCTTGTCCTc                                                                                       >  4:249278/1‑147 (MQ=255)                                                                       cACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGgaag                                                                                >  2:334373/1‑148 (MQ=255)                                                                        aCCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGgaagaa                                                                              <  8:308613/149‑1 (MQ=255)                                                                         ccGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGaagaag                                                                             <  1:334373/149‑1 (MQ=255)                                                                            ccATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAggg                                                                           >  1:138581/1‑148 (MQ=255)                                                                              aTCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCAc                                                                        >  6:357025/1‑149 (MQ=255)                                                                                 aCCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCt                                                                      <  3:79503/148‑1 (MQ=255)                                                                                  ccGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGTAAATGACGTGAATGGCAGCCAGCCTTTGc                                                                                                                         >  8:302145/1‑96 (MQ=255)                                                                                  ccGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGTAAATGACGTGAATGGCAGCCAGCCTTTGc                                                                                                                         <  7:251927/96‑1 (MQ=255)                                                                                  ccGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGTAAATGACGTGAATGGCAGCCAGCCTTTGc                                                                                                                         <  7:302145/96‑1 (MQ=255)                                                                                  ccGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGTAAATGACGTGAATGGCAGCCAGCCTTTGc                                                                                                                         >  8:251927/1‑96 (MQ=255)                                                                                            tACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCATGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGc                                                          >  7:342216/1‑149 (MQ=255)                                                                                              cGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGc                                                          <  8:25850/147‑1 (MQ=255)                                                                                              cGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCg                                                        <  1:143938/149‑1 (MQ=255)                                                                                              cGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCg                                                        <  4:338745/149‑1 (MQ=255)                                                                                                  tttGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGgag                                                     >  3:370324/1‑148 (MQ=255)                                                                                                  tttGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGgag                                                     >  3:364340/1‑148 (MQ=255)                                                                                                  tttGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATCCCGgag                                                     >  6:203795/1‑148 (MQ=255)                                                                                                     gggTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGaa                                                   >  6:257230/1‑147 (MQ=255)                                                                                                     ggcggATTCAGCGCTTGGGCCGGAGACAAGCCTGCCATCCGTGATGGCAAATGACCTGAATGGCAGCCAGCCTTTGCGGTAGCGAGCACCCTGACTCACCAACTGGGCCCCAAGGAAGAACTGCACCTGGCTTTTCATGCCGGAGAAc                                                  <  6:200075/144‑1 (MQ=255)                                                                                                             gAGGGTTGTCTCCGTTCACAAGCCATCCAACCGAGATGGCAAATGACGTGAATGGTAGCCTCCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAATGGCACCTGGCTTTTCATGCCGGAGAACAG‑ttcttc                                         <  8:155996/148‑1 (MQ=255)                                                                                                              aGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTcgc                                       >  1:107414/1‑149 (MQ=255)                                                                                                                ggTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGa                                     <  4:184709/149‑1 (MQ=255)                                                                                                                       gCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAgaga                              <  4:188752/149‑1 (MQ=255)                                                                                                                              aCAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGt                       <  3:53030/149‑1 (MQ=255)                                                                                                                                aaGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGATGa                     >  3:143872/1‑149 (MQ=255)                                                                                                                                        cATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAGTTTCTTCGCGATCAGAGAAGCCGGTGATGTTCtt              <  1:237054/149‑1 (MQ=255)                                                                                                                                         aTCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGc            >  4:200229/1‑149 (MQ=255)                                                                                                                                          tCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATTAGAGAAGCCGGTGATGTTCTTGcc           <  4:370324/149‑1 (MQ=255)                                                                                                                                          tCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATTAGAGAAGCCGGTGATGTTCTTGcc           <  4:364340/149‑1 (MQ=255)                                                                                                                                          tCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGcc           >  6:31214/1‑149 (MQ=255)                                                                                                                                              tGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTtg       <  5:71158/149‑1 (MQ=255)                                                                                                                                                aTGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTtgat     <  2:107414/149‑1 (MQ=255)                                                                                                                                                aTGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTtga      <  1:4578/148‑1 (MQ=255)                                                                                                                                                aTGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TGCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTtga      <  5:203795/148‑1 (MQ=255)                                                                                                                                                   gCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTTgatgat  <  2:65111/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGCTTTGGGTGATTGAGGGTTTTGGCCGGTGACAAGCCTTCCATCCGTGATGGCAAATGACGTGAATGGCAGCCAGCCTTTGCGGTACCGAGCACCCTGACTCACCAACTTGTCCTCAAGGAAGAAGGGCACCTGGCTTTTCATGCCGGAGAACAG‑TTCTTCGCGATCAGAGAAGCCGGTGATGTTCTTGCCCTTGATGAT  >  NC_002947/3389645‑3389939

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