Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,069,392	C→T	6.0%	V923I (GTC→ATC)	mdtB ←	multidrug transporter membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,069,392	0	C	T	6.0%	197.2 / 2.3	66	V923I (GTC→ATC)	mdtB	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base T (2/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTACAGCACCC  >  NC_002947/4069249‑4069538                                                                                                                                                |                                                                                                                                                   accAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATAc                                                                                                                                               <  2:62469/149‑1 (MQ=255)  ccAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATAcc                                                                                                                                              >  6:92724/1‑149 (MQ=255)   cAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATAcc                                                                                                                                              >  8:78378/1‑148 (MQ=255)            cgGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATcagcag                                                                                                                                     >  2:35501/1‑148 (MQ=255)            cgGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAgg                                                                                                                                    <  1:35501/149‑1 (MQ=255)             ggCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGga                                                                                                                                   >  3:296907/1‑149 (MQ=255)                cAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGgatga                                                                                                                                >  5:330735/1‑149 (MQ=255)                   ggtggtCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGatgatg                                                                                                                              >  7:76357/1‑148 (MQ=255)                   ggtggtCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGatgatg                                                                                                                              <  8:76357/148‑1 (MQ=255)                       gtCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCga                                                                                                                          <  8:40497/148‑1 (MQ=255)                            caGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGATGATACCGATCAGCAGGATGATGCCGATGATgg                                                                                                                    >  7:110509/1‑149 (MQ=255)                              ggATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGc                                                                                                                   >  2:1440/1‑148 (MQ=255)                               gATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGa                                                                                                                 >  2:336005/1‑149 (MQ=255)                               gATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGa                                                                                                                 >  7:182814/1‑149 (MQ=255)                                aTCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGAt                                                                                                                >  2:166620/1‑149 (MQ=255)                                      cGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGcc                                                                                                          >  1:93330/1‑149 (MQ=255)                                       ggAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCGGATGATGGCGATCATGcc                                                                                                          >  6:271706/1‑148 (MQ=255)                                             gcagcaACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCg                                                                                                   >  7:151819/1‑149 (MQ=255)                                              cagcaACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGt                                                                                                  <  7:23399/149‑1 (MQ=255)                                              cagcaACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATAGCCTGCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGt                                                                                                  >  7:282080/1‑149 (MQ=255)                                                gcaACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGgat                                                                                                                                  >  3:211796/1‑115 (MQ=255)                                                gcaACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGgat                                                                                                                                  <  4:211796/115‑1 (MQ=255)                                                     gccgccTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCt                                                                                           >  2:188457/1‑149 (MQ=255)                                                           tGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAg                                                                                     >  7:159749/1‑149 (MQ=255)                                                              tAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCaa                                                                                  <  8:151819/149‑1 (MQ=255)                                                                 aTGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGc                                                                               >  7:111771/1‑149 (MQ=255)                                                                   ggCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGGGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGc                                                                               >  7:221399/1‑147 (MQ=255)                                                                         gcgcgGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAg                                                                        <  3:199071/148‑1 (MQ=255)                                                                           gcgGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCa                                                                     <  8:159749/149‑1 (MQ=255)                                                                              ggCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGATGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCaccg                                                                  <  8:110509/149‑1 (MQ=255)                                                                                       ccTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCt                                                                                           <  5:165567/115‑1 (MQ=255)                                                                                       ccTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCt                                                                                           >  6:165567/1‑115 (MQ=255)                                                                                       ccTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGAccgccgcc                                                         >  6:171124/1‑149 (MQ=255)                                                                                       ccTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGAccgccgcc                                                         >  4:221907/1‑149 (MQ=255)                                                                                            ttGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACgg                                                    >  5:270430/1‑149 (MQ=255)                                                                                               cgcTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAg                                                 >  8:117365/1‑149 (MQ=255)                                                                                                gcTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAgg                                                <  8:221399/149‑1 (MQ=255)                                                                                                    ggCCTCCAGGGCGAAGTCGATAATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTc                                             >  8:194380/1‑148 (MQ=255)                                                                                                      ccTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGATGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTc                                             >  1:217492/1‑146 (MQ=255)                                                                                                      ccTCCAGGGCGAACTCGATCATCATGATCGCGTTCTTCTTGACCATACCGATAAGCAGGATGATGCCGATCAAGCCGATCATGCCCAGGTAGTTAACGCTGATGAGCAATGACAGCAACGCACCGCACCCCGCCCAAGCAAGGGGCgag                                          >  8:292121/1‑149 (MQ=255)                                                                                                       cTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGATGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCgaga                                         <  2:217492/149‑1 (MQ=255)                                                                                                       cTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCgaga                                         <  1:336005/149‑1 (MQ=255)                                                                                                       cTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCgaga                                         <  5:92724/149‑1 (MQ=255)                                                                                                        tCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGagag                                        >  2:156076/1‑149 (MQ=255)                                                                                                            gggCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATg                                    <  7:194380/149‑1 (MQ=255)                                                                                                             ggCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATgg                                   <  7:117365/149‑1 (MQ=255)                                                                                                             ggCGAAGTCGAACATCATGATCGCGTTCTTCTAGACGATACCGATCAGCACGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATgg                                   <  5:271706/149‑1 (MQ=255)                                                                                                                 aaGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGAt                               >  1:226632/1‑149 (MQ=255)                                                                                                                 aaGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGAt                               <  2:93330/149‑1 (MQ=255)                                                                                                                   gTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCg                             <  8:111771/149‑1 (MQ=255)                                                                                                                      gATCATCATGATCGCGTTCTTCTTGACGCTACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCGAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCggg                           <  7:216419/148‑1 (MQ=255)                                                                                                                      gATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGt                          <  1:1440/149‑1 (MQ=255)                                                                                                                         catcatGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGa                       >  3:65743/1‑149 (MQ=255)                                                                                                                           tcatGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATg                     <  3:221907/149‑1 (MQ=255)                                                                                                                             atGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTAACGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATg                     >  1:309540/1‑147 (MQ=255)                                                                                                                              tGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCTTGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAg                  <  1:268055/149‑1 (MQ=255)                                                                                                                              tGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTa                   <  4:65743/148‑1 (MQ=255)                                                                                                                                     gTTCTTCTTGACCATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTAGTTACCGCAGATGAGCAAGCCCAGCAAAGCACAGACAGCCGCAGACGGCAGCGTCGAGAGGATGGTGATCGCGTGGATGTAGCTTTCg            >  5:157968/1‑148 (MQ=255)                                                                                                                                       tcttctTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCtt               >  5:51064/1‑143 (MQ=255)                                                                                                                                       tcttctTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCtt               <  6:51064/143‑1 (MQ=255)                                                                                                                                        cttcttGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTa          >  3:326163/1‑147 (MQ=255)                                                                                                                                        cttcttGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTAc         <  4:296907/148‑1 (MQ=255)                                                                                                                                           cttGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTACa        <  5:171124/146‑1 (MQ=255)                                                                                                                                            ttGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCtt               >  5:94440/1‑138 (MQ=255)                                                                                                                                            ttGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCtt               <  6:94440/138‑1 (MQ=255)                                                                                                                                              gACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTACAGCAccc  >  2:218095/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   ACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGGCATCGCGCGGGCTCATGCCCTGATTGCGCTCGGCCTCCAGGGCGAAGTCGATCATCATGATCGCGTTCTTCTTGACGATACCGATCAGCAGGATGATGCCGATGATGGCGATCATGCCCAGGTCGTTACCGCTGATGAGCAAGGCCAGCAAGGCACCGACCGCCGCCGACGGCAGGGTCGAGAGGATGGTGATCGGGTGGATGTAGCTTTCGTACAGCACCC  >  NC_002947/4069249‑4069538

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