Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,799,826	G→T	6.7%	G143C (GGC→TGC)	PP_5664 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,799,826	0	G	T	6.7%	187.7 / 3.6	60	G143C (GGC→TGC)	PP_5664	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base T (2/2);  total (29/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.10e-01

TACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCGAGCACTTAA  >  NC_002947/4799682‑4799964                                                                                                                                                  |                                                                                                                                           tACGGAGAATATTCCGGTTTTTCATTAGACATCGTAGCCAACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCTA‑GGCTTCTTTCTATTGATCGAACGGTAGTCCAGCGACTTTgg                                                                                                                                           <  1:124676/146‑1 (MQ=255) tACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGc                                                                                                                                          <  2:185187/147‑1 (MQ=255)    ggAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGAtt                                                                                                                                      >  7:306875/1‑148 (MQ=255)     gagaAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTa                                                                                                                                    <  5:16770/149‑1 (MQ=255)     gagaAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTa                                                                                                                                    >  1:225487/1‑149 (MQ=255)      agaAAATTACGGTTTTTCATAAGACATCGTAGCCAACAGTTCAAAATACTTCGGCATAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGAttt                                                                                                                                     <  8:306875/147‑1 (MQ=255)          aaTTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGa                                                                                                                               >  2:167199/1‑149 (MQ=255)            ttCGGGTTTTTCATTAGACATCGTATCCAACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATc                                                                                                                             <  6:328874/149‑1 (MQ=255)               cGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTg                                                                                                                          <  2:225487/149‑1 (MQ=255)                ggTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTc                                                                                                                                 >  8:256974/1‑141 (MQ=255)                ggTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTc                                                                                                                                 <  7:256974/141‑1 (MQ=255)                ggTTTTTCATTAGACATCGTAGCCAACACTTCACAATACTGCGGCAGAGTGAATGAAATTGATTCAGGCAACATTTACAAGGAGGTGAAA‑GGCTTCTTTCTAATGATCTAACGGGAGTCCAGCGACTTTGGAGATTTATTAGATCCTg                                                                                                                          >  5:72948/1‑148 (MQ=255)                 gTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCa                                                                                                                        <  1:167199/149‑1 (MQ=255)                         ttAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACa                                                                                                                 >  2:109528/1‑148 (MQ=255)                           aGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAAGATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACaaa                                                                                                               <  7:181828/148‑1 (MQ=255)                            gACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAAt                                                                                                              >  7:36365/1‑148 (MQ=255)                               aTCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAgg                                                                                                          >  7:51201/1‑149 (MQ=255)                               aTCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAgg                                                                                                          >  1:285615/1‑149 (MQ=255)                                 cGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATc                                                                                                                             >  4:47244/1‑128 (MQ=255)                                 cGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATc                                                                                                                             <  3:47244/128‑1 (MQ=255)                                     gCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACt                                                                                                    <  8:51201/149‑1 (MQ=255)                                       ccACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGc                                                                                                                         <  3:107457/126‑1 (MQ=255)                                       ccACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGc                                                                                                                         >  4:107457/1‑126 (MQ=255)                                       ccACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTc                                                                                                   >  4:186299/1‑148 (MQ=255)                                        cacaCTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCaa                                                                                                 >  6:313717/1‑149 (MQ=255)                                         acacTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCaa                                                                                                 >  6:315524/1‑148 (MQ=255)                                         acacTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTCCATATAACAAATTAGGCAAACTCAAt                                                                                                <  1:177721/149‑1 (MQ=255)                                            cTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCat                                                                                             >  6:327839/1‑149 (MQ=255)                                             ttCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCatg                                                                                            <  3:186299/149‑1 (MQ=255)                                               cacaATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCatgat                                                                                          >  3:53019/1‑149 (MQ=255)                                                acaATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCatgat                                                                                          <  2:285615/148‑1 (MQ=255)                                                  aaTACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAAGGGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATACTGCATATAACAAATTAGGCAAACTCAATCAtgatga                                                                                        >  1:91192/1‑149 (MQ=255)                                                        gCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTAt                                                                                 >  2:239389/1‑149 (MQ=255)                                                         cGGCAGAGTGAATGAACTTGATGCAGGCAACATTGAGAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTAt                                                                                 <  2:139519/148‑1 (MQ=255)                                                                   aatgaaCTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTc                                                                       >  4:219123/1‑148 (MQ=255)                                                                    atgaaCTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGgttcgttc                                                                       >  4:1249/1‑147 (MQ=255)                                                                    atgaaCTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCa                                                                     <  4:53019/149‑1 (MQ=255)                                                                    atgaaCTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCa                                                                     <  5:327839/149‑1 (MQ=255)                                                                          ttGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTg                                                               <  5:315524/149‑1 (MQ=255)                                                                            gATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTgg                                                              <  1:109528/148‑1 (MQ=255)                                                                                            aCAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTg                                              <  3:316720/148‑1 (MQ=255)                                                                                            aCAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGc                                             <  8:36365/149‑1 (MQ=255)                                                                                                      gCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCatgat                                                                                          <  8:264962/94‑1 (MQ=255)                                                                                                      gCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCatgat                                                                                          >  7:264962/1‑94 (MQ=255)                                                                                                        aa‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCgg                                                       >  3:268234/1‑127 (MQ=255)                                                                                                        aa‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCgg                                                       <  4:268234/127‑1 (MQ=255)                                                                                                           ggCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCAtgatga                                                                                        >  6:11417/1‑92 (MQ=255)                                                                                                           ggCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCAtgatga                                                                                        <  5:11417/92‑1 (MQ=255)                                                                                                                     tATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGt                      >  2:316480/1‑148 (MQ=255)                                                                                                                     tATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGt                      <  5:313717/148‑1 (MQ=255)                                                                                                                       ttGATCGAACGGGAGTCCAGCGACTTTTGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCa                    <  3:199103/148‑1 (MQ=255)                                                                                                                       ttGATCGAACGGGAGTCCAGCGACTTTTGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCa                    <  3:195127/148‑1 (MQ=255)                                                                                                                       ttGATCGAACGGGAGTCCAGCGACTTTTGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCa                    >  4:195127/1‑148 (MQ=255)                                                                                                                       ttGATCGAACGGGAGTCCAGCGACTTTTGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCa                    >  4:199103/1‑148 (MQ=255)                                                                                                                         gATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGt                                                                                   <  5:167540/83‑1 (MQ=255)                                                                                                                         gATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGt                                                                                   >  6:167540/1‑83 (MQ=255)                                                                                                                              aaCGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGACCACCTAGAAAAAGTCAAAGCTGTc            >  6:117303/1‑149 (MQ=255)                                                                                                                                 gggAGTCCAGAGACTTTGGCGATTTATTCGATCCTTCATATAACAAATTAGGCAAACTCTATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCTAg         <  6:72948/149‑1 (MQ=255)                                                                                                                                        cAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCGAGCACTTaa  >  1:305712/1‑149 (MQ=255)                                                                                                                                        cAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCGAGCACTTaa  <  2:305712/149‑1 (MQ=255)                                                                                                                                        cAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCGAGCACTTaa  <  1:127862/149‑1 (MQ=255)                                                                                                                                                  ggCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGc                                     >  7:208862/1‑104 (MQ=255)                                                                                                                                                  ggCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGc                                     <  8:208862/104‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           TACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAATACTGCGGCAGAGTGAATGAACTTGATGCAGGCAACATTGACAAGGAGGTGCAA‑GGCTTCTTTCTATTGATCGAACGGGAGTCCAGCGACTTTGGCGATTTATTCGATCCTGCATATAACAAATTAGGCAAACTCAATCATGATGAGATGTATGGCTTCGTTCCAGCCTTGGCATTCGGTGGCCCCTGCGATCTTGCCCACCTAGAAAAAGTCAAAGCTGTCGAGCACTTAA  >  NC_002947/4799682‑4799964

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