Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,917,215	A→G	7.7%	K333E (AAG→GAG)	mexB →	multidrug resistance protein MexB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,917,215	0	A	G	7.7%	184.0 / 4.8	65	K333E (AAG→GAG)	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

ACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGT  >  NC_002947/3917074‑3917360                                                                                                                                               |                                                                                                                                                  aCGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTTGACGAACGGCGAGGTGtc                                                                                                                                            <  4:168934/149‑1 (MQ=255) aCGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACAGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGtc                                                                                                                                            <  6:127968/149‑1 (MQ=255)     ccgtccgCGTGCAGTTGTCACCCGGCGCCAACGCCATTGAAACCGCTACCCTGGTACGTGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCTTCAAGGTGtcgatc                                                                                                                                        <  8:251235/143‑1 (MQ=255)     ccTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTTGACGAACGGCGAGGTGtcgtag                                                                                                                                        <  2:217380/149‑4 (MQ=255)      cTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGtcgatc                                                                                                                                        >  6:97430/1‑148 (MQ=255)       ttCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGCCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGtcgatcga                                                                                                                                      >  1:96906/1‑149 (MQ=255)          aGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCa                                                                                                                                   <  5:97430/149‑1 (MQ=255)          aGCGTGCAGGTGGCACCCGGGGCCAACGCCATGGAAACCGGTTCCCTGGTTCGGGCGAAGATGCAGGCCCTGGCGGGCTACTTCCGCGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAg                                                                                                                                     <  1:114139/148‑1 (MQ=255)           gCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAg                                                                                                                                  <  6:266206/149‑1 (MQ=255)           gCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAg                                                                                                                                  <  8:238062/149‑1 (MQ=255)              tGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtc                                                                                                                               <  2:359712/149‑1 (MQ=255)              tGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtc                                                                                                                               <  6:27819/149‑1 (MQ=255)               gCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTCGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTCCGACAACCCCAACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtca                                                                                                                              >  6:17087/1‑149 (MQ=255)                cAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CTGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcat                                                                                                                             <  8:197952/149‑1 (MQ=255)                            ggCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGtt                                                                                                                 <  3:227087/149‑1 (MQ=255)                            ggCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGCCGCGCTACTTCC‑CGGAAGCCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAACGTGACGATCGAGCAGGTCATCAACACCCTGta                                                                                                                 >  4:317735/1‑148 (MQ=255)                                ccAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                                                                                                             <  4:152760/149‑1 (MQ=255)                                ccAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                                                                                                             >  2:37968/1‑149 (MQ=255)                                  aaCGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGc                                                                                                           >  3:269832/1‑149 (MQ=255)                                       cATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCAtgc                                                                                                      <  6:320615/149‑1 (MQ=255)                                           gAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCAtgctgc                                                                                                   <  5:204974/148‑1 (MQ=255)                                           gAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATgctgct                                                                                                  >  8:95319/1‑149 (MQ=255)                                                    aCCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAg                                                                                                                                  >  4:163711/1‑108 (MQ=255)                                                    aCCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAg                                                                                                                                  <  3:163711/108‑1 (MQ=255)                                                     cccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGGAGGtcatc                                                                                                                            <  1:90036/113‑1 (MQ=255)                                                     cccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtcatc                                                                                                                            >  2:90036/1‑113 (MQ=255)                                                      ccTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAg                                                                                       <  5:213069/149‑1 (MQ=255)                                                           tACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGc                                                                                                           >  6:309895/1‑124 (MQ=255)                                                           tACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGc                                                                                                           <  5:309895/124‑1 (MQ=255)                                                              gggCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttc                                                                               <  4:269832/149‑1 (MQ=255)                                                                gCGTAGATGCAGGACCTGGCGCGCTACTTGC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCTGTCAtgttcct                                                                             <  3:317735/149‑1 (MQ=255)                                                                  gAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttcctgt                                                                           >  6:214734/1‑149 (MQ=255)                                                                  gAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttcctgt                                                                           <  6:321668/149‑1 (MQ=255)                                                                     gATGCAGGACCTGGCGCGCTACTTCC‑CGGTAGGCGTCAAGTACGACATCCGCTACGACAGCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttcctgttcc                                                                        <  5:17087/149‑1 (MQ=255)                                                                            gACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGATACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGaa                                                                 >  3:245011/1‑149 (MQ=255)                                                                                  gcgcgcTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGATACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTgc                                                            <  4:245011/148‑1 (MQ=255)                                                                                  gcgcgcTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTgcg                                                           >  5:146419/1‑149 (MQ=255)                                                                                   cgcgcTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTgcgc                                                          <  5:109560/149‑1 (MQ=255)                                                                                             cc‑cGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTg                                                 <  5:214734/148‑1 (MQ=255)                                                                                                cGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATc                                              >  6:175410/1‑149 (MQ=255)                                                                                                    aGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGa                                          <  1:37968/149‑1 (MQ=255)                                                                                                      gCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACCCCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGAcc                                        >  5:355309/1‑149 (MQ=255)                                                                                                        gTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCt                                      >  5:277809/1‑149 (MQ=255)                                                                                                         tCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTg                                     >  1:142204/1‑149 (MQ=255)                                                                                                            aGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCTTGTTTGATTCAATGATGCTGGTTTTATCAGTAATTTTCCTTTTCCTTAAGAAACTGCCCTACACGCTGATCCCGACCCTGGTc                                  <  6:207510/149‑1 (MQ=255)                                                                                                              tACGACATCCCCTACGACACATCGCCGTTCGTCAACGTGTAGATCTAGCACGTCATCAACACCCTGTTCGAGGCCATGCTCCTGGTGTTCGCAGTCATCTTCCTGTTCCTGCAGAACCTGCGCTACACGATGATCCCGAACCTGGTCg                                 >  2:179361/1‑148 (MQ=255)                                                                                                                    aTCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGt                          >  8:262251/1‑149 (MQ=255)                                                                                                                      ccccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTgg                        >  5:219249/1‑149 (MQ=255)                                                                                                                       cccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                                    >  8:193239/1‑104 (MQ=255)                                                                                                                       cccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                                    <  7:193239/104‑1 (MQ=255)                                                                                                                       cccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCATGTTAGAGGCCATGCTGCTGGTGTTCGCAGTCATCTTCCTGTTCATGCAGAACCTCCGCTACACGCTGATCCCGACCCTTGTCGTGCAATTgg                        >  3:294655/1‑148 (MQ=255)                                                                                                                         cTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAt                     <  7:95319/149‑1 (MQ=255)                                                                                                                         cTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAt                     <  7:6394/149‑1 (MQ=255)                                                                                                                           tcGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAtt                    >  6:239978/2‑148 (MQ=37)                                                                                                                           tcGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAtt                    >  4:363356/2‑148 (MQ=37)                                                                                                                           tcGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTCTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAtt                    >  1:287597/2‑148 (MQ=37)                                                                                                                           aCGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTg                         >  4:27332/1‑143 (MQ=255)                                                                                                                           aCGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTg                         <  3:27332/143‑1 (MQ=255)                                                                                                                             gACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGAt                 >  1:145303/1‑149 (MQ=255)                                                                                                                              acacCTCGCCGTTCGTCAAGGTGTCGATCGATCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATg                <  4:210847/149‑1 (MQ=255)                                                                                                                              acacCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATg                >  5:382896/1‑149 (MQ=255)                                                                                                                                  cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCa            >  7:294725/1‑149 (MQ=255)                                                                                                                                  cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCa            >  3:142532/1‑149 (MQ=255)                                                                                                                                            gTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGt  >  2:55433/1‑149 (MQ=255)                                                                                                                                            gTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGt  >  3:123819/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                  ACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGT  >  NC_002947/3917074‑3917360

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