Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,209,214	C→T	6.9%	Y389Y (TAC→TAT)	PP_3692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,209,214	0	C	T	6.9%	172.9 / 2.2	58	Y389Y (TAC→TAT)	PP_3692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/21);  new base T (2/2);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

GAATGATGCTGCCGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTG  >  NC_002947/4209067‑4209359                                                                                                                                                    |                                                                                                                                                   gAATGATGCTGCCGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACg                                                                                                                                                   <  4:169787/149‑1 (MQ=255)  aaTGATGCAGCCGCACTTGAGCTAGTCAACTAGACTGAAAAAGCAGCGAGAACACGCAAAGGCAGTGAAGCGGCTCAGGCCGCTATTTTCGGTAACGCAGAGCTGGAGGGCGTTGGCTCATGGACGTGGCAGCAGATCAGGGAGTACGc                                                                                                                                                  >  6:98896/1‑149 (MQ=255)            ccGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑c                                                                                                                                       >  8:115477/1‑149 (MQ=255)            ccGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑c                                                                                                                                       >  2:269829/1‑149 (MQ=255)             cGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑Ct                                                                                                                                      >  6:127214/1‑149 (MQ=255)                    gAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTccgccgc                                                                                                                               >  1:222886/1‑149 (MQ=255)                          gTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTGCGCCGCCGTGGc                                                                                                                         >  1:129072/1‑149 (MQ=255)                                ttGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGCAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCATGGCGtctac                                                                                                                   >  6:284512/1‑145 (MQ=255)                                     tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTATGCGCGAGCGTA‑CTcccaca                                                                                                                                <  1:228147/131‑5 (MQ=255)                                     tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTATGCGCGAGCGTA‑CTcccaca                                                                                                                                >  2:228147/1‑127 (MQ=255)                                     tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGt                                                                                                              >  2:82401/1‑149 (MQ=255)                                        aaaaGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCAc                                                                                                           <  1:269829/149‑1 (MQ=255)                                          aaGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCAccc                                                                                                         >  5:280683/1‑149 (MQ=255)                                            gCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCtt                                                                                                       <  1:135375/149‑1 (MQ=255)                                                        cGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATAGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCCCCCTTCCACATCgtgt                                                                                            >  4:55017/1‑148 (MQ=255)                                                         gCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTcc                                                                                                     <  8:88499/138‑1 (MQ=255)                                                         gCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTcc                                                                                                     >  7:88499/1‑138 (MQ=255)                                                          cAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCATAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATACCTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCa                                                                                          <  1:82401/149‑1 (MQ=255)                                                           aaaGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGc                                                                                        <  2:222886/149‑1 (MQ=255)                                                                 aGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTGACCCTTCCCCAACGTGTCAGCTGAGTc                                                                                  <  6:280683/149‑1 (MQ=255)                                                                 aGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGt                                                                                   >  3:269857/1‑148 (MQ=255)                                                                            tGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGc                                                                       <  3:55017/149‑1 (MQ=255)                                                                                  cGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGc                                                                  >  4:302422/1‑148 (MQ=255)                                                                                       ttttCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAg                                                            >  3:21840/1‑149 (MQ=255)                                                                                       ttttCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAg                                                            <  1:165054/149‑1 (MQ=255)                                                                                            ggTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtct                                                       >  7:138395/1‑149 (MQ=255)                                                                                             gTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGCGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtctc                                                      >  4:310985/1‑149 (MQ=255)                                                                                             gTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtctc                                                      >  5:315190/1‑149 (MQ=255)                                                                                              tAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGTGGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCATGGCGTATCCTGTGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcg                                                     <  5:284512/149‑1 (MQ=255)                                                                                               aaGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcga                                                    >  6:116337/1‑149 (MQ=255)                                                                                                  gCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcgacga                                                 >  5:94260/1‑149 (MQ=255)                                                                                                   cAGAGCTGGGGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCa                                                             <  8:268655/136‑1 (MQ=255)                                                                                                   cAGAGCTGGGGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCa                                                             >  7:268655/1‑136 (MQ=255)                                                                                                        cTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                           <  7:115477/149‑1 (MQ=255)                                                                                                           gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCttt                                                                    >  4:225040/1‑121 (MQ=255)                                                                                                           gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAt                                                         >  5:272064/1‑132 (MQ=255)                                                                                                           gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAt                                                         <  6:272064/132‑1 (MQ=255)                                                                                                            aGGGCTTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCaaaa                                       <  7:202836/149‑1 (MQ=255)                                                                                                                 gTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTCGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACg                                  >  6:11738/1‑149 (MQ=255)                                                                                                                         aTCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGcc                           <  2:129072/148‑1 (MQ=255)                                                                                                                         aTCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGcc                           >  6:217665/1‑148 (MQ=255)                                                                                                                            gACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGaa                        <  5:217665/148‑1 (MQ=255)                                                                                                                            gACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAAt                       >  5:143122/1‑149 (MQ=255)                                                                                                                              cGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAAttt                     <  7:201087/149‑1 (MQ=255)                                                                                                                                 ggCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTg                                                                   >  8:276373/1‑100 (MQ=255)                                                                                                                                    agcagATGTGGGAGTACGCTCTCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAg               >  6:127703/1‑149 (MQ=255)                                                                                                                                    agcagATGTGGGAGTACGCTCGCGCATT‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAg               >  8:228339/1‑149 (MQ=255)                                                                                                                                       agATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGATTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCAGCAAAAGCACGCCTGGCCGAATTTGAAAAGttt            <  5:127214/149‑1 (MQ=255)                                                                                                                                         aTGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCACCGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGt          <  4:117394/149‑1 (MQ=255)                                                                                                                                         aTGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTTTCAGCTGAGTCACGATGCGTGCTTTTCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTg           >  7:231076/1‑148 (MQ=255)                                                                                                                                           gtgGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTc         >  7:180950/1‑148 (MQ=255)                                                                                                                                             gggAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                           <  1:221138/112‑1 (MQ=255)                                                                                                                                             gggAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                           >  2:221138/1‑112 (MQ=255)                                                                                                                                             cggAGTATGCGCGAGCGTA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGt                                                                         <  6:28322/81‑1 (MQ=38)                                                                                                                                             cggAGTATGCGCGAGCGTA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGt                                                                         >  5:28322/2‑82 (MQ=38)                                                                                                                                               gAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgc    <  4:269857/149‑1 (MQ=255)                                                                                                                                                 gTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGCTCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgctg  >  6:282167/1‑149 (MQ=255)                                                                                                                                                 gTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgctg  >  4:118847/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAATGATGCTGCCGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTG  >  NC_002947/4209067‑4209359

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