Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,501	C→G	8.6%	A306G (GCC→GGC)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,501	0	C	G	8.6%	174.3 / 7.4	58	A306G (GCC→GGC)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/3);  total (27/31)
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.22e-01

CGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTC  >  NC_002947/4605354‑4605643                                                                                                                                                     |                                                                                                                                               cGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑Gcc                                                                                                                                               >  7:124229/1‑149 (MQ=255)              cccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGcc                                                                                                                                   >  6:297540/1‑148 (MQ=255)                 acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTAGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAgg                                                                                                                               <  6:264436/149‑1 (MQ=255)                  caGGTGGGCATCGAGTTCCATTTCACGCGTCACCTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                              <  4:234132/149‑1 (MQ=255)                         gCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                        >  6:259021/1‑148 (MQ=255)                            tCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAg                                                                                                                    <  8:124229/149‑1 (MQ=255)                              gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                              <  4:224288/137‑1 (MQ=255)                              gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GGCAGCCTCGCTGGCCAGGt                                                                                                                              >  3:224288/1‑137 (MQ=255)                              gAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGcc                                                                                                                  <  5:259021/149‑1 (MQ=255)                                    cATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCAGAGAAGATTTCAGGCTGTTCTGCAGCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGTCCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                             <  1:55966/148‑1 (MQ=255)                                      tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCg                                                                                                          <  2:262042/149‑1 (MQ=255)                                         cACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAg                                                                                                        <  5:105633/148‑1 (MQ=255)                                                       cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                         <  3:122129/149‑1 (MQ=255)                                                       cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                         <  5:297540/149‑1 (MQ=255)                                                        aGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGACGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGTCGCAGTCCGCTAAc                                                                                         >  8:263412/1‑148 (MQ=255)                                                         gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                            <  2:237030/144‑1 (MQ=255)                                                         gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCt                                                                                            >  1:237030/1‑144 (MQ=255)                                                         gACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACg                                                                                        >  8:24035/1‑148 (MQ=255)                                                          acccccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaa                                                                                      <  8:278372/146‑1 (MQ=255)                                                           cacCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGaaa                                                                                     >  5:88636/1‑149 (MQ=255)                                                              ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                  >  1:240711/1‑149 (MQ=255)                                                              ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCt                                                                                  >  5:306820/1‑149 (MQ=255)                                                                        ggCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                       >  8:289127/1‑102 (MQ=255)                                                                        ggCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTg                                                                                                                       <  7:289127/102‑1 (MQ=255)                                                                         gCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                       <  4:157111/149‑1 (MQ=255)                                                                           acaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                             >  5:277060/1‑109 (MQ=255)                                                                           acaGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGaa                                                                                                             <  6:277060/109‑1 (MQ=255)                                                                             agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGa                                                                                                                                          >  4:144118/1‑78 (MQ=255)                                                                             agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGa                                                                                                                                          <  3:144118/78‑1 (MQ=255)                                                                             agaagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAg                                                                   >  6:99951/1‑149 (MQ=255)                                                                               aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAGGt                                                                  >  8:98419/1‑148 (MQ=255)                                                                               aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtc                                                                 <  6:88636/149‑1 (MQ=255)                                                                                  aTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                              <  1:251464/149‑1 (MQ=255)                                                                                   tttCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatca                                                             <  5:99951/149‑1 (MQ=255)                                                                                      cAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAGGTCATCAGc                                                           <  7:98419/148‑1 (MQ=255)                                                                                         ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                       >  4:242076/1‑149 (MQ=255)                                                                                         ccTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAccc                                                       >  5:113088/1‑149 (MQ=255)                                                                                             ttCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                     <  3:242076/147‑1 (MQ=255)                                                                                                 gCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                >  1:264789/1‑148 (MQ=255)                                                                                                   aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAc                                                         <  8:36987/137‑1 (MQ=255)                                                                                                   aCCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCAc                                                         >  7:36987/1‑137 (MQ=255)                                                                                                     cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAg                                                            <  1:252020/132‑1 (MQ=255)                                                                                                     cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAg                                                            >  2:252020/1‑132 (MQ=255)                                                                                                             tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTTTCGATGAAGTCATAAGCACCCGAACAACCACTGACGGaa                                     <  7:263412/147‑1 (MQ=255)                                                                                                             tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGaa                                     <  2:264789/147‑1 (MQ=255)                                                                                                                 ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                               >  1:117718/1‑149 (MQ=255)                                                                                                                 ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                               >  2:165167/1‑149 (MQ=255)                                                                                                                 ccTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCg                               >  3:151470/1‑149 (MQ=255)                                                                                                                    gTTCAAGCACAGTGCCGAGCCCATCGACCTGCGCCAGCGAGGCTGCCCAGGTCGAGATGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGt                             >  1:188828/1‑149 (MQ=255)                                                                                                                                    gAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGATCAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACttt            >  6:8708/1‑149 (MQ=255)                                                                                                                                      gCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTc          <  6:306820/149‑1 (MQ=255)                                                                                                                                        ccATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTAAGTGTCGATGAAGTCATCCGACCCCGAACAACCACGGACGGAAGTACCGGGGCGCACCGTCGGACATTgcgc        >  4:228774/1‑149 (MQ=255)                                                                                                                                         cATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGcc       <  2:240711/149‑1 (MQ=255)                                                                                                                                          aTCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGAc                                         <  7:245857/114‑1 (MQ=255)                                                                                                                                          aTCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGAc                                         >  8:245857/1‑114 (MQ=255)                                                                                                                                             gACCTG‑GCCAGCGAGGCTGCCCAGTTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt   <  3:228774/149‑1 (MQ=255)                                                                                                                                             gACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt   <  7:24035/149‑1 (MQ=255)                                                                                                                                              aCCTG‑GGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGCGCACCTGCGGACTTTTCGCCCTTTc  >  3:44188/1‑149 (MQ=38)                                                                                                                                              aCCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc  <  6:113088/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTG‑GCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTC  >  NC_002947/4605354‑4605643

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