Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	401,371	G→T	6.0%	R319S (AGG→AGT)	PP_0334 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	401,371	0	G	T	6.0%	198.5 / 3.3	67	R319S (AGG→AGT)	PP_0334	transposase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (2/2);  total (35/32)
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.58e-01

TGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCC  >  NC_002947/401253‑401515                                                                                                                       |                                                                                                                                                 tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                                                                                    >  1:239414/1‑149 (MQ=32)   atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGc                                                                                                                  >  2:387117/1‑149 (MQ=32)          aGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGtt                                                                                                            >  5:241953/1‑148 (MQ=32)          aGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                                                                                           >  1:9287/1‑149 (MQ=32)           gggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTAc                                                                                                          <  7:298032/149‑1 (MQ=32)             gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTa                                                                                                           >  5:295642/1‑146 (MQ=32)              cgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGTTCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTTATCGAAATCTGCGGATTATCTGGGCCCTGTTACaga                                                                                                       <  1:11981/149‑1 (MQ=17)               gcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACagag                                                                                                      <  6:295642/149‑1 (MQ=32)                 gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGc                                                                                                     >  5:88951/1‑148 (MQ=32)                   aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCt                                                                                                               <  8:136609/136‑1 (MQ=32)                   aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCt                                                                                                               >  7:136609/1‑136 (MQ=32)                    ttCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGACGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGAt                                                                                                  >  6:364612/1‑148 (MQ=17)                       cGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAg                                                                                               <  1:410135/148‑1 (MQ=32)                        gTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGt                                                                                             <  5:284858/149‑1 (MQ=32)                           ttGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTa                                                                                          >  2:93651/1‑149 (MQ=32)                             gAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCGTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTAc                                                                                         <  2:9287/148‑1 (MQ=17)                             gAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGttt                                                                                           <  6:253573/146‑1 (MQ=32)                                 cAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                                                                                    <  7:227874/149‑1 (MQ=32)                                  aGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                                                                                    <  6:206803/116‑1 (MQ=32)                                     tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGc                                                                                                                  >  2:349637/1‑115 (MQ=32)                                     tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGc                                                                                                                  <  1:349637/115‑1 (MQ=32)                                     tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCAAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCaa                                                                                 >  8:344244/1‑148 (MQ=17)                                          ccTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                                                                           <  6:241953/149‑1 (MQ=32)                                              aTTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacga                                                                        >  1:128617/1‑148 (MQ=32)                                                    aTCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                 <  1:343441/149‑1 (MQ=32)                                                      cAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                 >  7:223600/1‑147 (MQ=32)                                                        gCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAgg                                                             <  2:128617/149‑1 (MQ=32)                                                           aGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                          >  8:187251/1‑149 (MQ=32)                                                             atgatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcggc                                                         >  2:269648/1‑148 (MQ=32)                                                                ATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTAt                                                      >  6:88951‑M1/1‑148 (MQ=255)                                                                  CGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAg                                                   <  8:157085‑M1/149‑1 (MQ=255)                                                                   GCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGc                                                  >  1:269648‑M1/1‑149 (MQ=255)                                                                     cgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCct                                                >  3:167089‑M2/50‑149 (MQ=255)                                                                     cgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTGATCGGAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctc                                               >  8:69241‑M2/50‑149 (MQ=255)                                                                     CGCCACTTGCGCTGGCTATACCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCGGGGTAATAGAAATCTGCGGATTATCTGGGCCCt                                                                                                               <  2:250708‑M1/86‑1 (MQ=255)                                                                      gcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctc                                               >  5:262437‑M2/49‑149 (MQ=255)                                                                       ctccttgcgcccaaggcgggatagccagcgcaaggaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGCTTTACCAACGCAAACCGAACGAAAAACCGGACGCGGCTATGACCctct                                              >  4:140166‑M2/48‑149 (MQ=255)                                                                       cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTa                                                       <  5:263122‑M2/93‑1 (MQ=255)                                                                       cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTa                                                       >  6:263122‑M2/48‑140 (MQ=255)                                                                       CCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctct                                              <  4:324791‑M1/149‑1 (MQ=255)                                                                       CCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCctct                                              <  3:408959‑M1/149‑1 (MQ=255)                                                                        ctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTGATCGGAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                            <  7:69241‑M2/103‑1 (MQ=255)                                                                              gcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTg                                             <  2:100315‑M2/102‑1 (MQ=255)                                                                              gcgcccaaggcgggatagccagcgcaagaggcggcgatca‑TCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTg                                             >  1:100315‑M2/41‑142 (MQ=255)                                                                               gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                               >  3:235642/1‑124 (MQ=32)                                                                               gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                               <  4:235642/124‑1 (MQ=32)                                                                               gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTg                                             <  4:163942/142‑1 (MQ=32)                                                                               gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTg                                             >  3:163942/1‑142 (MQ=32)                                                                               gcTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAAcc                                      >  3:162776/1‑149 (MQ=32)                                                                                   gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCggg                                   <  6:99304/148‑1 (MQ=32)                                                                                   gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTGTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCggg                                   >  2:254815/1‑148 (MQ=17)                                                                                    cTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGt                                  >  1:12640/1‑148 (MQ=32)                                                                                     tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTaa                                >  2:3586/1‑149 (MQ=32)                                                                                         ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTa                                 <  5:144928/144‑1 (MQ=32)                                                                                         ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTa                                 >  6:144928/1‑144 (MQ=32)                                                                                          cGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCa                           >  4:143884/1‑149 (MQ=32)                                                                                               tGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCggg                      >  8:374682/1‑149 (MQ=32)                                                                                                gggCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGTCACCggg                      <  5:364612/148‑1 (MQ=17)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                                                                               >  2:207650/1‑85 (MQ=33)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                                                                               <  1:207650/85‑1 (MQ=33)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATg                                                     <  3:125262/111‑1 (MQ=32)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATg                                                     >  4:125262/1‑111 (MQ=32)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAg               <  7:157085/149‑1 (MQ=32)                                                                                                      aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAg               <  2:12640/149‑1 (MQ=32)                                                                                                               gcgatcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCaa      >  7:78061‑M2/8‑149 (MQ=255)                                                                                                                   tcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAcccc  <  7:344244‑M2/146‑1 (MQ=255)                                                                                                                   tcaGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAcccc  <  8:223600‑M2/146‑1 (MQ=255)                                                                                                                       |                                                                                                                                                 TGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCC  >  NC_002947/401253‑401515

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