Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,917,209	T→A	7.6%	F331I (TTC→ATC) ‡	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,209	0	T	A	7.6%	186.2 / 5.1	66	F331I (TTC→ATC) ‡	mexB	multidrug resistance protein MexB
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

GGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCA  >  NC_002947/3917062‑3917350                                                                                                                                                     |                                                                                                                                              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                                                                                                 >  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: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               >  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)                                                                                                                                                     |                                                                                                                                              GGCAAACCGGCCACGGCCTTCAGCGTGCAGTTGTCACCCGGCGCCAACGCCATGGAAACCGCTACCCTGGTACGGGCGAAGATGCAGGACCTGGCGCGCTACTTCC‑CGGAAGGCGTCAAGTACGACATCCCCTACGACACCTCGCCGTTCGTCAAGGTGTCGATCGAGCAGGTCATCAACACCCTGTTCGAGGCCATGCTGCTGGTGTTCGCAGTCATGTTCCTGTTCCTGCAGAACCTGCGCTACACGCTGATCCCGACCCTGGTCGTGCCAGTGGCATTGATGGGCA  >  NC_002947/3917062‑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: ‑