Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,917,204	C→T	7.2%	S329L (TCG→TTG)	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,204	0	C	T	7.2%	202.6 / 4.0	69	S329L (TCG→TTG)	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (3/2);  total (35/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.43e-01

GAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCA  >  NC_002947/3917058‑3917350                                                                                                                                                    |                                                                                                                                                   gAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGc                                                                                                                                                  >  5:321668/1‑149 (MQ=255)  aaCGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGcc                                                                                                                                                 >  3:171308/1‑149 (MQ=255)     ggCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGtt                                                                                                                                              >  3:210847/1‑149 (MQ=255)        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:321668/149‑1 (MQ=255)                                                                                  gAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCAtgttcctgt                                                                 >  6:214734/1‑149 (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                                                          <  7:193239/104‑1 (MQ=255)                                                                                                                                       cccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCa                                                          >  8:193239/1‑104 (MQ=255)                                                                                                                                       cccTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCATGTTAGAGGCCATGCTGCTGGTGTTCGCAGTCATCTTCCTGTTCATGCAGAACCTCCGCTACACGCTGATCCCGACCCTTGTCGTGCAATTgg              >  3:294655/1‑148 (MQ=255)                                                                                                                                         cTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAt           <  7:6394/149‑1 (MQ=255)                                                                                                                                         cTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCAt           <  7:95319/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               <  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)                                                                                                                                                    |                                                                                                                                                   GAACGGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCA  >  NC_002947/3917058‑3917350

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