Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,260,613	G→C	6.4%	T69T (ACC→ACG)	PP_2857 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,260,613	0	G	C	6.4%	191.6 / 4.0	63	T69T (ACC→ACG)	PP_2857	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (2/2);  total (32/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 = 3.87e-01

TGCTCAGTGCGTTGGCTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCGTTGTTGATGACCGGAACCAGG  >  NC_002947/3260465‑3260760                                                                                                                                                      |                                                                                                                                                      tGCTCAGTGCGTTGGCTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTgg                                                                                                                                                       <  2:162645/149‑1 (MQ=255)    tCAGTGCGTTGGCTTGCTGGGTGTTGGCACCTTCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTg                                                                                                                                                    >  3:156411/1‑149 (MQ=255)           gTTGGCTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCAt                                                                                                                                             >  2:161964/1‑149 (MQ=255)               gCTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGAGCCGCTCCCGTTGCTGAAAGCGATGC‑CCTCGATGGTTGAGCGGGCGTCAATCGCGGGGTCGCCGGGCGAGTCGAGCCGTGGTCTGATCTGGGTGCTTGCATTGGc                                                                                                                                         >  6:214836/1‑149 (MQ=255)                cTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCt                                                                                                                                        >  3:257791/1‑149 (MQ=255)                 ttGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCt‑t                                                                                                                                      >  1:212710/1‑149 (MQ=255)                     tGGGTGTTGGCACCTGTGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCggg                                                                                                                                  >  2:284895/1‑149 (MQ=255)                             gtcACCTGCGCTCTGGTGGCCGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGGCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTg                                                                                                                                           <  6:36178/131‑1 (MQ=255)                                      gcTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGAT‑GTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGgcgcg                                                                                                                >  1:166571/1‑149 (MQ=255)                                           ggTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCa                                                                                                            >  8:57471/1‑149 (MQ=255)                                            gTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑t                                                                                                          >  6:257536/1‑149 (MQ=255)                                                  gCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGgcgc                                                                                                                 >  6:247686/1‑137 (MQ=255)                                                  gCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGgcgc                                                                                                                 <  5:247686/137‑1 (MQ=255)                                                    ccAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGCACCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACctg                                                                                                   >  1:207267/1‑149 (MQ=255)                                                          cgccgcTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGt                                                                                            <  4:156411/149‑1 (MQ=255)                                                            ccgcTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGttt                                                                                          <  7:253330/149‑1 (MQ=255)                                                            ccgcTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGttt                                                                                          <  1:161964/149‑1 (MQ=255)                                                              gcTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGt                                                                                        >  6:189487/1‑149 (MQ=255)                                                                  ccGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCTCTCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTc                                                                                       <  1:73724/147‑1 (MQ=255)                                                                  ccGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCAcc                                                                                    >  3:197816/1‑149 (MQ=255)                                                                    gTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTg                                                                                  >  2:160487/1‑149 (MQ=255)                                                                     ttGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACctg                                                                                                   <  8:279230/131‑1 (MQ=255)                                                                     ttGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACctg                                                                                                   >  7:279230/1‑131 (MQ=255)                                                                        cTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGTCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCgg                                                                               <  4:257791/148‑1 (MQ=255)                                                                        cTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGc                                                                              <  1:160487/149‑1 (MQ=255)                                                                          gAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCt                                                                            >  8:108865/1‑149 (MQ=255)                                                                            aaGCGTTGC‑CCTGGATGGGTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTgg                                                                          <  2:207267/149‑1 (MQ=255)                                                                               cGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCATTTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGtt                                                                        <  1:215558/149‑1 (MQ=255)                                                                                 ttGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGa                                                                      >  2:62155/1‑148 (MQ=255)                                                                                 ttGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGAc                                                                     <  1:140486/149‑1 (MQ=255)                                                                                   gC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACa                                                                    <  7:108865/148‑1 (MQ=255)                                                                                   gC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACa                                                                    <  1:284895/148‑1 (MQ=255)                                                                                    gcgcTGGAGGGTTGAGCGGGCGTCCATCGCGGGGGCGGCGGGCGTGTGGCGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACag                                                                   <  2:128644/146‑1 (MQ=255)                                                                                      ccTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACa                                                                    >  6:82420/1‑146 (MQ=255)                                                                                       cTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCt                                                                                   <  7:251100/130‑1 (MQ=255)                                                                                       cTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCt                                                                                   >  8:251100/1‑130 (MQ=255)                                                                                       cTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACaga                                                                  >  7:74945/1‑147 (MQ=255)                                                                                        tGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGCGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGtt                                                                                           >  4:143058/1‑121 (MQ=255)                                                                                        tGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGCGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGtt                                                                                           >  6:327804/1‑121 (MQ=255)                                                                                        tGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGCGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGtt                                                                                           <  3:143058/121‑1 (MQ=255)                                                                                        tGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGCGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGtt                                                                                           <  5:327804/121‑1 (MQ=255)                                                                                         ggATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGACGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGAc                                                                     <  4:98152/142‑1 (MQ=255)                                                                                         ggATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGACGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGAc                                                                     >  3:98152/1‑142 (MQ=255)                                                                                              gTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGcc                                                         <  5:82420/149‑1 (MQ=255)                                                                                               ttGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGccc                                                        >  6:234684/1‑149 (MQ=255)                                                                                                  aGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGa                                                     <  1:62155/149‑1 (MQ=255)                                                                                                      ggCGTCCATCTCGGGGTAGGCGGGCGTGGCGAGCCTGTGTCTGATCTGGTTGATTGCATTGGCT‑TCGGTGCCGACGCCGATGGCGCGAGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGc                                                 <  2:212710/149‑1 (MQ=255)                                                                                                            cATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGc                                           >  5:14327/1‑149 (MQ=255)                                                                                                             aTCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGAc                                                                     >  6:266133/1‑122 (MQ=255)                                                                                                              tggcgGGGTCGGCGGGCGTGTCGAGGCGTTGTCTGTTATGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGGGATGGCGCGCGCT‑TAGAGCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGa                                                     <  3:192078/135‑1 (MQ=255)                                                                                                              tCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCg                                         >  7:3427/1‑149 (MQ=255)                                                                                                              tCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCg                                         <  7:57471/149‑1 (MQ=255)                                                                                                              tCGCGGGGTCCGCGGGCGACTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGCCA‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTTACCAGCTGCTCTTTGTCACCTGCGGCCTGGATGACAGAGAGCTTGCCCTGa                                                     >  4:192078/1‑137 (MQ=255)                                                                                                                gcgGGGGCGGCGGGCGTGTCGAGCCGTTGGCTGATCTGGGGGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGa                                       <  5:214836/149‑1 (MQ=255)                                                                                                                     gTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGt                                  <  5:189487/149‑1 (MQ=255)                                                                                                                             gcgTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCg                          <  6:14327/149‑1 (MQ=255)                                                                                                                                  tCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCgtt                     >  3:130873/1‑149 (MQ=255)                                                                                                                                  tCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCgtt                     >  5:191577/1‑149 (MQ=255)                                                                                                                                    gAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCgttgt                   >  4:50018/1‑149 (MQ=255)                                                                                                                                    gAGCCGTTGTCTGATCTGGGTGATTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGACAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCgttgt                   <  4:197816/149‑1 (MQ=255)                                                                                                                                     aGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCgttgtt                  <  5:257536/149‑1 (MQ=255)                                                                                                                                      gCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCgg                                               >  3:32582/1‑119 (MQ=255)                                                                                                                                      gCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCgg                                               <  4:32582/119‑1 (MQ=255)                                                                                                                                          ttGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCGTTGTtgatga             <  4:130873/149‑1 (MQ=255)                                                                                                                                                      ggTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCGTTGTTGATGACCGGAACCAgg  >  3:63716/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGCTCAGTGCGTTGGCTTGCTGGGTGTTGGCACCTGCGCTCTGGTTGACGCCCAGGGCGCCGCTGCCGTTGCTGAAAGCGTTGC‑CCTGGATGGTTGAGCGGGCGTCCATCGCGGGGTCGGCGGGCGTGTCGAGCCGTTGTCTGATCTGGGTGCTTGCATTGGCT‑TCGGGGCCGACGGCGATGGCGCGCGCA‑TTGACCTGCTGCTGTTTGTCACCTGCGGCCTGGTTGACAGAGAGGTTGCCCTGATAGCGGCTGCCGGAGCTGTCGATATCGGCGTTGTTGATGACCGGAACCAGG  >  NC_002947/3260465‑3260760

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑