Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,917,213	T→G	7.7%	V332G (GTC→GGC)	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,213	0	T	G	7.7%	185.9 / 7.3	65	V332G (GTC→GGC)	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base T (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 = 3.96e-01

AAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGT  >  NC_002947/3917065‑3917360                                                                                                                                                      |                                                                                                                                                    aaaCCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGt                                                                                                                                                     <  5:289126/149‑1 (MQ=255)  aaCCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGt                                                                                                                                                     >  5:320615/1‑148 (MQ=255)  aaCCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGt                                                                                                                                                     >  6:213069/1‑148 (MQ=255)          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                                                                                                                               <  6:27819/149‑1 (MQ=255)                       tGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGtc                                                                                                                               <  2:359712/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                                                                                                             >  2:37968/1‑149 (MQ=255)                                         ccAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAg                                                                                                             <  4:152760/149‑1 (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                                                                                                                                  <  3:163711/108‑1 (MQ=255)                                                             aCCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAg                                                                                                                                  >  4:163711/1‑108 (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                                                                                                           <  5:309895/124‑1 (MQ=255)                                                                    tACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGc                                                                                                           >  6:309895/1‑124 (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                    >  4:363356/2‑148 (MQ=37)                                                                                                                                    tcGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAtt                    >  6:239978/2‑148 (MQ=37)                                                                                                                                    tcGACACCTTGACGAACGGCGAGGTGTCGATCGAGCAGGTCATCAACACCCTCTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAtt                    >  1:287597/2‑148 (MQ=37)                                                                                                                                    aCGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTg                         <  3:27332/143‑1 (MQ=255)                                                                                                                                    aCGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTg                         >  4:27332/1‑143 (MQ=255)                                                                                                                                      gACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGAt                 >  1:145303/1‑149 (MQ=255)                                                                                                                                       acacCTCGCCGTTCGTCAAGGTGTCGATCGATCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATg                <  4:210847/149‑1 (MQ=255)                                                                                                                                       acacCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATg                >  5:382896/1‑149 (MQ=255)                                                                                                                                           cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCa            >  3:142532/1‑149 (MQ=255)                                                                                                                                           cTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCa            >  7:294725/1‑149 (MQ=255)                                                                                                                                                     gTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGt  >  2:55433/1‑149 (MQ=255)                                                                                                                                                     gTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGt  >  3:123819/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCACCTTCGCCGT  >  NC_002947/3917065‑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: ‑