Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	901,419	A→T	10.4%	intergenic (‑138/‑193)	cysZ ← / → trxB	sulfate transport protein CysZ/thioredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	901,419	0	A	T	10.4%	148.1 / 8.4	58	intergenic (‑138/‑193)	cysZ/trxB	sulfate transport protein CysZ/thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base T (4/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

GCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGC  >  NC_002947/901274‑901565                                                                                                                                                        |                                                                                                                                                   gCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGc                                                                                                                                                 >  6:232284/1‑149 (MQ=255)    tGCAGCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTa                                                                                                                                              <  4:21746/149‑1 (MQ=255)       aTCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa                                                                                                                                           >  4:278650/1‑149 (MQ=255)         cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAgg                                                                                                                                         >  2:63860/1‑149 (MQ=255)           agagACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTg                                                                                                                                       >  8:50086/1‑149 (MQ=255)            gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTgg                                                                                                                                      <  3:278650/149‑1 (MQ=255)            gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTg                                                                                                                                       >  5:47264/1‑148 (MQ=255)             agaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa                                                                                                                                     >  1:343978/1‑149 (MQ=255)                  ccGCAGAGAAGGAACACGCGCTGACCTTACCGAGTTTGCCGCGCCGACCAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTc                                                                                                                                <  7:66549/149‑1 (MQ=255)                   ctcAGAGAATGAAATCTCGCTGTCCTTAACGAGTTTGCCGCTCCGTAGAAATGTCACTAGCCGTTGGCCTACAGTAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt                                                                                                                                        <  1:279795/144‑1 (MQ=255)                   cGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGATTGTTTATAGATGAAATCGCAGCATAGGCGATGCATATGAGGTGGATTGTc                                                                                                                                >  8:68847/1‑148 (MQ=255)                     cAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCt                                                                                                                               >  1:203907/1‑147 (MQ=255)                       gagaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgctgc                                                                                                                                        >  3:312269/1‑132 (MQ=255)                       gagaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgctgc                                                                                                                                        <  4:312269/136‑5 (MQ=255)                          aaGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACAGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGata                                                                                                                        >  1:276574/1‑149 (MQ=255)                                    cgcTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAAt                                                                                                              >  7:194223/1‑149 (MQ=255)                                         aCCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTTTAAGGATATTTCCTTAATCttt                                                                                                          <  7:50086/148‑1 (MQ=255)                                          ccTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTg                                                                                                         <  4:32671/148‑1 (MQ=255)                                          ccTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc                                                                                                        <  1:63860/149‑1 (MQ=255)                                               ccGAGTTTGCCGTGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCt                                                                                                   >  2:328370/1‑149 (MQ=255)                                                                acgaAA‑GCCAC‑‑GCGGTGGGC‑‑‑GGGAAAGGGGATGAAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCaa                                                                                  <  2:343978/149‑1 (MQ=255)                                                                  gaAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt                                                                                >  7:29031/1‑149 (MQ=255)                                                                  gaAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt                                                                                >  6:95648/1‑149 (MQ=255)                                                                    aa‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt                                                                                          >  1:167712/1‑137 (MQ=255)                                                                    aa‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt                                                                                          <  2:167712/137‑1 (MQ=255)                                                                     a‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                          >  1:241470/1‑120 (MQ=255)                                                                     a‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt                                                                                                          <  2:241470/120‑1 (MQ=255)                                                                        gCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCg                                                                           >  7:285201/1‑149 (MQ=255)                                                                          aC‑‑GGCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCt                                                                         <  5:27056/149‑1 (MQ=255)                                                                                  tGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                   <  7:68847/149‑1 (MQ=255)                                                                                  tGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg                                                                   <  2:217308/149‑1 (MQ=255)                                                                                          gggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTg                                                              >  8:196809/1‑149 (MQ=255)                                                                                                  ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTAAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGtttt                                                      <  2:276574/149‑1 (MQ=255)                                                                                                   gATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTAAATCTTTGCGACCTCGCTACTGTGCGCTCACCTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                     <  6:47264/149‑1 (MQ=255)                                                                                                        aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc                                                                                                        >  8:138002/1‑93 (MQ=255)                                                                                                        aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc                                                                                                        <  7:138002/93‑1 (MQ=255)                                                                                                         aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                     >  6:135787/1‑143 (MQ=255)                                                                                                         aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa                                                     <  5:135787/143‑1 (MQ=255)                                                                                                           aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTccc                                             >  6:267729/1‑149 (MQ=255)                                                                                                              aGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa                                          <  2:203907/149‑1 (MQ=255)                                                                                                                         cTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg                               <  5:267729/149‑1 (MQ=255)                                                                                                                                   gATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCttttt                      >  2:58635/1‑148 (MQ=255)                                                                                                                                      gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                               <  5:78106/120‑1 (MQ=255)                                                                                                                                      gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc                                               >  6:78106/1‑120 (MQ=255)                                                                                                                                        aaTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTa                     <  4:319152/144‑1 (MQ=255)                                                                                                                                        aaTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTa                     >  3:319152/1‑144 (MQ=255)                                                                                                                                            ggagctTACGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTACTTAATCTTTGCGATCTCTCTACAGTGCGCTGAACTTTTCTCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCTAGTGCTGCGTGGGTCCTTTTTATTTTCTtg             <  7:196809/146‑1 (MQ=255)                                                                                                                                            gcagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtgg            <  5:232284/149‑1 (MQ=255)                                                                                                                                              ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAttt                 >  5:240557/3‑142 (MQ=255)                                                                                                                                              ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAttt                 <  6:240557/140‑1 (MQ=255)                                                                                                                                              ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct          >  7:89424/3‑149 (MQ=255)                                                                                                                                              ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct          >  4:56377/3‑149 (MQ=255)                                                                                                                                               gcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg                               >  3:331082/1‑127 (MQ=255)                                                                                                                                               gcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg                               <  4:331082/127‑1 (MQ=255)                                                                                                                                                caTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggctgg        <  1:79070/149‑1 (MQ=255)                                                                                                                                                     gCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAg   >  6:246212/1‑149 (MQ=255)                                                                                                                                                     gCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAg   >  8:101177/1‑149 (MQ=255)                                                                                                                                                      cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGc  <  8:194223/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   GCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGC  >  NC_002947/901274‑901565

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