Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,318,924	G→A	5.7%	A381T (GCC→ACC)	PP_3789 →	efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,318,924	0	G	A	5.7%	205.5 / 2.4	70	A381T (GCC→ACC)	PP_3789	efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base A (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

ATTGACTCGATCAGCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGC  >  NC_002947/4318780‑4319071                                                                                                                                                 |                                                                                                                                                    aTTGACTCGATCAGCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGc                                                                                                                                                <  1:488198/150‑1 (MQ=255)  ttGACTCGATCAGCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCg                                                                                                                                               >  2:165383/1‑150 (MQ=255)       tCGATCAGCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCTATCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCCGGGTTTGTGCTTCTCACTGTTCGCCTGCGTGCTg                                                                                                                                          >  2:1241242/1‑150 (MQ=255)              gCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCGTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTTTCACTGATCGCCTGCGTGCGGGGCGGTa                                                                                                                                   >  2:398846/1‑150 (MQ=255)                ccgccgCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTa                                                                                                                                   >  1:1062448/1‑148 (MQ=255)                ccgccgCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTa                                                                                                                                   <  2:1062448/148‑1 (MQ=255)                 cgccgcCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTagca                                                                                                                                <  2:123554/150‑1 (MQ=255)                       cccTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGCTTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCACCTGCGTGCTGGCGGGTAGCAGCCttt                                                                                                                          >  2:943570/1‑150 (MQ=255)                       cccTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCACCTGCGTGCTGGCGGGTAGCAGCCttt                                                                                                                          >  2:943535/1‑150 (MQ=255)                        ccTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGCGGCAGCCTTtg                                                                                                                         >  1:1051816/1‑150 (MQ=255)                          tCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTtgtg                                                                                                                       >  2:1182569/1‑150 (MQ=255)                          tCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTtgtg                                                                                                                       >  1:627147/1‑150 (MQ=255)                           cAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGATTGCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGc                                                                                                                      <  1:1005639/150‑1 (MQ=255)                              aGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCAt                                                                                                                   <  2:241488/150‑1 (MQ=255)                                  aGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATgcgc                                                                                                               >  1:1615739/1‑150 (MQ=255)                                  aGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATgcgc                                                                                                               >  1:1278640/1‑150 (MQ=255)                                     tATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCt                                                                                                            <  1:1159098/150‑1 (MQ=255)                                              gcACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGAcc                                                                                                   >  2:1726815/1‑150 (MQ=255)                                              gcACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGAcc                                                                                                   >  1:394663/1‑150 (MQ=255)                                                 cacaATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCg                                                                                                >  1:842359/1‑150 (MQ=255)                                                        tGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACt                                                                                         >  2:1902398/1‑150 (MQ=255)                                                            gAGCTTGGGGGGGGTAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGTGCCCCCAGGGTTTGTCCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGCCCCCGGGCAACTTTgc                                                                                     <  1:398846/150‑1 (MQ=255)                                                            gAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTgc                                                                                     >  1:1911799/1‑150 (MQ=255)                                                            gAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTgc                                                                                     <  2:1278640/150‑1 (MQ=255)                                                            gAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTgc                                                                                     >  2:246064/1‑150 (MQ=255)                                                            gAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTgc                                                                                     <  2:627147/150‑1 (MQ=255)                                                               cTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGcc                                                                                                                     <  1:1810160/115‑1 (MQ=255)                                                               cTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGcc                                                                                                                     >  2:1810160/1‑115 (MQ=255)                                                                 tGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                <  2:156949/150‑1 (MQ=255)                                                                 tGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                <  2:1615739/150‑1 (MQ=255)                                                                 tGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCggg                                                                                <  2:1054765/150‑1 (MQ=255)                                                                        ggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGc                                                                                                                                                    >  1:1667509/1‑75 (MQ=255)                                                                        ggcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGc                                                                                                                                                    <  2:1667509/75‑1 (MQ=255)                                                                         gcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCACCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCtt                                                                        <  1:943535/150‑1 (MQ=255)                                                                         gcAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCACCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCtt                                                                        <  1:943570/150‑1 (MQ=255)                                                                             tCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAAt                                                                    <  1:1182569/150‑1 (MQ=255)                                                                               aGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATAc                                                                  <  1:246064/150‑1 (MQ=255)                                                                                 cccGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACtt                                                                <  2:842359/150‑1 (MQ=255)                                                                                         ggtgTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTAcgc                                                        <  1:165383/150‑1 (MQ=255)                                                                                            gtgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCa                                                     >  2:112943/1‑150 (MQ=255)                                                                                             tgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTtgt                                                                                                                        >  2:1581131/1‑82 (MQ=255)                                                                                             tgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTtgt                                                                                                                        <  1:1581131/82‑1 (MQ=255)                                                                                             tgGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCAt                                                    <  2:1911799/150‑1 (MQ=255)                                                                                               gATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTg                                                  <  1:1902398/150‑1 (MQ=255)                                                                                                  tCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTCATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGgtc                                               >  1:1265919/1‑150 (MQ=255)                                                                                                    aTACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGgtcgt                                             <  2:1829184/150‑1 (MQ=255)                                                                                                    aTACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGgtcgt                                             >  1:1367141/1‑150 (MQ=255)                                                                                                        tATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAAt                                         >  2:1617168/1‑150 (MQ=255)                                                                                                         aTCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATg                                        >  1:89398/1‑150 (MQ=255)                                                                                                          tcgctcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCgg                                                                                               >  1:30776/1‑94 (MQ=255)                                                                                                          tcgctcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCgg                                                                                               <  2:30776/94‑1 (MQ=255)                                                                                                          tcgctcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGt                                       >  2:1178257/1‑150 (MQ=255)                                                                                                             ctcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCa                                                              <  1:399694/124‑1 (MQ=255)                                                                                                             ctcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCa                                                              >  2:399694/1‑124 (MQ=255)                                                                                                              tcgcGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCa                                   >  1:75574/1‑150 (MQ=255)                                                                                                               cgcgGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCAc                                  <  1:112943/150‑1 (MQ=255)                                                                                                                    gCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCa                             <  1:1178257/150‑1 (MQ=255)                                                                                                                       cccAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGCCCCCGGGCAACTTTGCGCGGGGAGGTCTTTAATACTTCAGTACGCGCATTgg                                                 >  1:692242/1‑127 (MQ=255)                                                                                                                       cccAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGGTCTTTAATACTTCAGTACGCGCATTgg                                                 <  2:692242/127‑1 (MQ=255)                                                                                                                       cccAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACtt                          >  1:1535598/1‑150 (MQ=255)                                                                                                                        ccAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTAGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTCCGCGGGGAGATCTTTAATACTTCAGTAAGCGCATTGGTCGTCAATGTTTCACCCTAAATtc                         >  2:268172/1‑150 (MQ=255)                                                                                                                         cAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCt                                                                                                                            <  1:1388294/50‑1 (MQ=255)                                                                                                                         cAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCt                                                                                                                            >  2:1388294/1‑50 (MQ=255)                                                                                                                         cAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTgg                                                 >  1:246272/1‑125 (MQ=255)                                                                                                                         cAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTgg                                                 <  2:246272/125‑1 (MQ=255)                                                                                                                              tttGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTctctacc                   >  1:439539/1‑147 (MQ=255)                                                                                                                                     ttCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTg            <  1:1241242/150‑1 (MQ=255)                                                                                                                                      tctcACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCCCACTTCTCTGATCATTCTgg           <  1:1726815/150‑1 (MQ=255)                                                                                                                                       ctcACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGt          >  2:523768/1‑150 (MQ=255)                                                                                                                                               tCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGc  >  2:1626257/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ATTGACTCGATCAGCCCGCCGCCCCTCAAAGGCAGTTATTTCGCCGCACACAATCTGGCGAGCTTGGGGGCGGCAGTCAGCCCGGCGTGGTGTGGATTCATACTATCGCTCGCGGGCCCCCAGGGTTTGTGCTTCTCACTGATCGCCTGCGTGCTGGCGGGTAGCAGCCTTTGTGCCATGCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGC  >  NC_002947/4318780‑4319071

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑