Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,389,501	T→C	6.8%	T414A (ACG→GCG)	PP_3865 ←	tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,389,501	0	T	C	6.8%	168.5 / 2.1	59	T414A (ACG→GCG)	PP_3865	tail protein
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base C (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 = 9.84e-01

GTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGA  >  NC_002947/4389358‑4389647                                                                                                                                                |                                                                                                                                                   gtAGGCCTTGATCAGCTCCTACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCGCGGCGATCAGGCCAATAGGCCTGGCCGCGCAGGCCCCACCGAGCACCCCCGTGGccc                                                                                                                                               >  4:25984/1‑148 (MQ=255)    ggCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCAcag                                                                                                                                            >  7:252324/1‑147 (MQ=255)    ggCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGACCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACtg                                                                                                                                            >  8:110207/1‑149 (MQ=255)       cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGGGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                         <  6:1311/149‑1 (MQ=255)       cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGCTAAGGCCAATCGGCGTGGCGGCGAAGGCCCCACCGCGCACCCGCGGGGGCCCTGTgg                                                                                                                                         >  3:262452/1‑149 (MQ=255)       cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                         >  2:268833/1‑149 (MQ=255)       cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACAGTgg                                                                                                                                         <  8:252324/149‑1 (MQ=255)           aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgc                                                                                                                                     <  4:122275/149‑1 (MQ=255)           aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcg                                                                                                                                      <  1:268833/148‑1 (MQ=255)              aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGaa                                                                                                                                  <  2:299068/149‑1 (MQ=255)              aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGAGCAGGCCAAACGGCGTGGCCGCCAAGGCCACACCCAGCACCCCCGCGGCCCCCGGGGCcccgcg                                                                                                                                  >  4:221872/1‑147 (MQ=255)              aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCCCGCTCCCCGCGGCGAGCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCCCCCCCGCCGCCACACTGGCcacgcg                                                                                                                                  >  4:217434/1‑147 (MQ=255)               gCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                 <  3:217434/149‑4 (MQ=255)               gCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                 <  3:221872/149‑4 (MQ=255)                cTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACAATCAGACCGGCCGCCCCCCCCAAACCCCCGGCGAAAAGGCCAAAACGCGTGGCCGCGAACGCCACACCGCGCACCCCCCTTGCCAAAGTGGCCGCGAaca                                                                                                                                >  8:180081/1‑149 (MQ=255)                     acaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGGCGCGTTCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGGCACTGTGGCCGCGAACAccgcc                                                                                                                           <  6:262979/149‑1 (MQ=255)                        gggCGCGGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACAATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGGGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCcgca                                                                                                                        <  7:180081/149‑1 (MQ=255)                              tGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCGCGGCGATCAGGCCAAACGGCCTGGCCGCGAAGGCCCCACCCAGCCCCCGCCTGGCCACTGTGGCCCCGAACACCGCCGCACGCAgc                                                                                                                  >  7:9985/1‑149 (MQ=255)                              tGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCCATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCCCGAACACCGCCGCACGCAgc                                                                                                                  >  7:4735/1‑149 (MQ=255)                                aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACAGCCGCACGCAgcgc                                                                                                                >  5:360166/1‑149 (MQ=255)                                       gCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                         <  3:25984/149‑1 (MQ=255)                                        cccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAgg                                                                                                        >  2:58545/1‑149 (MQ=255)                                              cTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCCCTCCCCCCGGCGCTCAGGCCAATCGGCGTGGCCGCGCAGGCCACACCGAGCACCCCCGTGGCCACAGTGGCCGCGAACACCCCCGCACGCAGCGCGGTGAAGGCCGccc                                                                                                  >  3:69933/1‑149 (MQ=255)                                                    ccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                                                            <  1:58545/149‑1 (MQ=255)                                                     cAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATACCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAt                                                                                           >  6:74682/1‑149 (MQ=255)                                                                gctCAGACCGCCCGCCGCCGCGATCCCCCCGGCGATCACGCCAATCGGCCTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCcgc                                                                                <  4:262452/147‑1 (MQ=255)                                                                gaTCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAac                                                                                                                                 <  5:294840/100‑1 (MQ=255)                                                                    aGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCa                                                                             <  3:166852/148‑1 (MQ=255)                                                                     gACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCAc                                                                            >  1:231495/1‑148 (MQ=255)                                                                             gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                   <  8:97577/149‑1 (MQ=255)                                                                             gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                   >  6:314556/1‑149 (MQ=255)                                                                             gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                   <  5:314556/149‑1 (MQ=255)                                                                                gccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                <  1:117682/149‑1 (MQ=255)                                                                                         cccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                        <  4:69933/148‑1 (MQ=255)                                                                                          tcgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatc                                                      <  7:210386/148‑1 (MQ=255)                                                                                          tcgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatc                                                      <  7:204958/148‑1 (MQ=255)                                                                                          ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGat                                                       >  4:305552/1‑148 (MQ=255)                                                                                           cgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                     <  3:305552/149‑1 (MQ=255)                                                                                              ggcgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCt                                                  <  6:360166/149‑1 (MQ=255)                                                                                                cgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgc                                                <  7:128396/149‑1 (MQ=255)                                                                                                 gATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgcg                                               <  7:110207/149‑1 (MQ=255)                                                                                                       gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                         >  6:324373/1‑149 (MQ=255)                                                                                                        ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                         >  5:67584/1‑132 (MQ=255)                                                                                                        ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                         <  6:67584/132‑1 (MQ=255)                                                                                                            tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACg                                                                                       <  8:76851/98‑1 (MQ=255)                                                                                                            tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACg                                                                                       >  7:76851/1‑98 (MQ=255)                                                                                                              ggCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                  <  5:156450/149‑1 (MQ=255)                                                                                                                cGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATcc                                <  5:74682/149‑1 (MQ=255)                                                                                                                  tGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                              <  2:231495/149‑1 (MQ=255)                                                                                                                      cgcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGcc                          >  4:353607/1‑149 (MQ=255)                                                                                                                           aGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGt                     <  8:9985/149‑1 (MQ=255)                                                                                                                           aGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGt                     <  8:4735/149‑1 (MQ=255)                                                                                                                              ccACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGc                  >  5:17666/1‑149 (MQ=255)                                                                                                                               cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                  <  4:136293/132‑1 (MQ=255)                                                                                                                                  acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCTCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                >  2:53791/1‑147 (MQ=255)                                                                                                                                    cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa              >  5:248510/1‑147 (MQ=255)                                                                                                                                    cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa              <  6:248510/147‑1 (MQ=255)                                                                                                                                           ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAt     >  1:177259/1‑149 (MQ=255)                                                                                                                                              cgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGa  >  4:205221/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   GTAGGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGA  >  NC_002947/4389358‑4389647

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