Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,181,395	C→A	6.2%	P198T (CCG→ACG)	PP_2791 →	aminoglycoside phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,181,395	0	C	A	6.2%	201.7 / 3.6	64	P198T (CCG→ACG)	PP_2791	aminoglycoside phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base A (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

TGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGTTGCACCGT  >  NC_002947/3181252‑3181541                                                                                                                                                    |                                                                                                                                                   tgCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGa                                                                                                                                                <  3:216454/148‑1 (MQ=255) tgCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGAc                                                                                                                                               <  5:366106/149‑1 (MQ=255)     gggTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGGAACGGACAgga                                                                                                                                           <  1:61937/149‑3 (MQ=255)          gcgcTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCg                                                                                                                                      <  5:163459/149‑1 (MQ=255)            gcTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCgg                                                                                                                                     >  1:243735/1‑148 (MQ=255)            gcTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCgg                                                                                                                                     <  2:243735/148‑1 (MQ=255)            gcTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCgg                                                                                                                                     <  2:249072/148‑1 (MQ=255)            gcTGCCGGCCGCGCGCTCAGCACCGGCGCCGGGCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCgg                                                                                                                                     >  1:249072/1‑148 (MQ=255)               gccggccgCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCt                                                                                                                                 >  6:97535/1‑149 (MQ=255)                cgggcggcGCGCTCTGAACGGGGGCGGGTCCTCTGGGCTGGGCCTAACCCACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCgg                                                                                                                                     <  1:316187/146‑1 (MQ=255)                 cggccgCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTgg                                                                                                                               <  5:51500/149‑1 (MQ=255)                    ccgCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCC‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCGGTGCTTCATGTCAGTTCCGGACCGCTCGACGGCCCTGGCgc                                                                                                                            >  4:230648/1‑148 (MQ=255)                     cgcgcgCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGAc                                                                                                                           >  4:185261/1‑149 (MQ=255)                     cgcgcgCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGAc                                                                                                                           <  5:97535/149‑1 (MQ=255)                        gcgcTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCAt                                                                                                                        <  7:293038/149‑1 (MQ=255)                           cTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGcc                                                                                                                                   <  3:129909/135‑1 (MQ=255)                           cTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGcc                                                                                                                                   >  4:129909/1‑135 (MQ=255)                                cACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGcc                                                                                                                 >  1:205161/1‑148 (MQ=255)                                 aCCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGGAACGGACAg                                                                                                                                             <  2:309638/119‑1 (MQ=255)                                 aCCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTCTGCTTCCTGTCCGGAACGGACAg                                                                                                                                             >  1:309638/1‑119 (MQ=255)                                         cGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACcgccgc                                                                                                        >  5:90516/1‑148 (MQ=255)                                          ggTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACACCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTa                                                                                                      >  2:20639/1‑149 (MQ=255)                                                cAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGATTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGACGTACGCCACGCTCTGCTTCCTGTCACTTCAGGACAGCTCGCACGCCCTGGCGAACATCGACGCCTACAGCCCCTACCTCGa                                                                                                >  2:32205/1‑149 (MQ=255)                                                 aGGGCCTGTCTCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGACCGGCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTTGCGACCTTCGACGCCTACCGCCGCTACCTCGac                                                                                               <  1:32205/149‑1 (MQ=255)                                                 aGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCg                                                                                                 <  3:129715/147‑1 (MQ=255)                                                 aGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGa                                                                                                >  3:287421/1‑148 (MQ=255)                                                   ggCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGacac                                                                                             >  2:305854/1‑149 (MQ=255)                                                   ggCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGacac                                                                                             <  3:189320/149‑1 (MQ=255)                                                             a‑‑‑‑aCTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGCCACCCTCGCCGAt                                                                                   >  5:306387/1‑149 (MQ=255)                                                                        gACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATc                                                                            >  3:125114/1‑149 (MQ=255)                                                                          cGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCg                                                                          <  1:305854/149‑1 (MQ=255)                                                                           gCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCt                                                                                                   <  7:344754/123‑1 (MQ=255)                                                                           gCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCt                                                                                                   >  8:344754/1‑123 (MQ=255)                                                                           gCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCgg                                                                         >  7:137069/1‑149 (MQ=255)                                                                             aGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTg                                                                       <  2:322179/149‑1 (MQ=255)                                                                                  ggCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTg                                                                    <  6:23083/147‑1 (MQ=255)                                                                                   gcgcCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGaa                                                                 >  8:157782/1‑149 (MQ=255)                                                                                     gcCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATg                                                               >  3:358717/1‑149 (MQ=255)                                                                                         aCCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGc                                                           <  2:205161/149‑1 (MQ=255)                                                                                         aCCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGc                                                           <  3:185261/149‑1 (MQ=255)                                                                                                   tCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGt                                                 <  1:20639/149‑1 (MQ=255)                                                                                                       cATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGa                                             <  4:287421/149‑1 (MQ=255)                                                                                                        aTCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAg                                            >  6:150225/1‑149 (MQ=255)                                                                                                            ggTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTg                                         >  1:244434/1‑148 (MQ=255)                                                                                                              tCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATggg                                                             >  6:369208/1‑126 (MQ=255)                                                                                                              tCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATggg                                                             <  5:369208/126‑1 (MQ=255)                                                                                                               ccgcccgcCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCt                                                                               >  1:106388/1‑107 (MQ=255)                                                                                                               ccgcccgcCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCt                                                                               <  2:106388/107‑1 (MQ=255)                                                                                                                cgcccgccGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCg                                                                                        <  3:9894/97‑1 (MQ=255)                                                                                                                cgcccgccGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCg                                                                                        >  4:9894/1‑97 (MQ=255)                                                                                                                  cccgccGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATgcgc                                  <  5:145992/149‑1 (MQ=255)                                                                                                                           cgcCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCgcagca                         <  8:137069/149‑1 (MQ=255)                                                                                                                             ccACGCTGTGCATCCTGTCCCTTCCGGACAGCTCGCCGGCCCTCGCGACCATCGACGCCTACAGCCCCTACCTATACACCCTCGCCGATCCCTATCCGGTGCGGGAATGGGGCCTGCGCTCGTGCGAGCTCCAAGAGCCGCgcagcaga                       >  4:55327/1‑148 (MQ=255)                                                                                                                             ccACGCTGAGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCgcagca                         <  3:55327/147‑1 (MQ=255)                                                                                                                              cACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTg                                                   <  6:181603/120‑1 (MQ=255)                                                                                                                              cACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTg                                                   >  5:181603/1‑120 (MQ=255)                                                                                                                                  cTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCt                  >  5:345362/1‑149 (MQ=255)                                                                                                                                    gtgCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCtgt                <  4:125114/149‑1 (MQ=255)                                                                                                                                          gcTGTCCGGAACGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGt          >  8:184833/2‑149 (MQ=255)                                                                                                                                           cTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGtt         <  5:150225/149‑1 (MQ=255)                                                                                                                                           cTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGtt         >  8:228984/1‑149 (MQ=255)                                                                                                                                           cTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGt          <  7:157782/148‑1 (MQ=255)                                                                                                                                              tccggtccggACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCTCCGATGCCTAGCCGGGGCGGGACTGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGTTg        <  3:230648/147‑1 (MQ=255)                                                                                                                                              tCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCt                  <  5:271171/137‑1 (MQ=255)                                                                                                                                                  ttCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGTTGCACCGt  >  3:184074/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGCCGGGTAGCGCTGCCGGCCGCGCGCTCAGCACCGGCGCCGGTCCTCAGGGCCGGGCGCA‑‑‑‑ACTGGCGACGCAGCTTGGCGCCAACCTGGCGCGTCTGCATCAGGTCCGCCCGCCGTGCGCCACGCTGTGCTTCCTGTCCGTTCCGGACAGCTCGCCGGCCCTGGCGACCATCGACGCCTACCGCCGCTACCTCGACACCCTCGCCGATGCCTATCCGGTGCTGGAATGGGGCCTGCGCTGGTGCGAGCTGCATGCGCCGCGCAGCAGCACCCTGTGCCTGTTGCACCGT  >  NC_002947/3181252‑3181541

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