Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,278,096	G→T	6.3%	L9I (CTT→ATT) ‡	ampC ←	beta‑lactamase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,278,096	0	G	T	6.3%	203.2 / 2.9	63	L9I (CTT→ATT) ‡	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base T (2/2);  total (33/30)
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.40e-01

GCGTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGTGAT  >  NC_002947/3277948‑3278240                                                                                                                                                     |                                                                                                                                                 gcgTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAg                                                                                                                                                  <  2:413292/149‑1 (MQ=255)              cATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGCGCGGATGATCTGGTCGACCTTGGCCTGTAACAGGTCGTCGGCCTCGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGt                                                                                                                                     >  2:219715/1‑149 (MQ=255)                tAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTgg                                                                                                                                   <  5:261016/149‑1 (MQ=255)                 aGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCACTTGCAGccgccatt                                                                                                                                     <  7:335765/146‑7 (MQ=255)                 aGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCACTTGCAGccgccatt                                                                                                                                     >  8:335765/1‑140 (MQ=255)                  gATGGCCACAGCCATGCCGGTGATACCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGc                                                                                                                                 <  3:416970/149‑1 (MQ=255)                      gCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAAGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGa                                                                                                                              <  1:1309/148‑1 (MQ=255)                          caGCCTTGCCGGTGATGCCTTGTTCCTGCTTGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATg                                                                                                                         <  8:190875/149‑1 (MQ=255)                              cATGCCGGTGAATACTTGTTCTGTAATGAGCGTTCGTATTATCTGTTCTACCTTTGCCTGTAACTGGTCGTCGTCCTGGGCGCTGGCAGGGCCTAAGTCAAGGGTAATGGCGTCTGCAAGTTCAGCCAGACGTGTGATGATCATGGtta                                                                                                                     <  2:97822/149‑3 (MQ=255)                                tGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAg                                                                                                                             >  6:107953/1‑139 (MQ=255)                                tGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAg                                                                                                                             <  5:107953/139‑1 (MQ=255)                                         tGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGATGCAAGTGCAGCCAGACGTGGGCGGAGCATTTTTAGTCCTTTtgtt                                                                                                          <  6:244013/149‑2 (MQ=255)                                          gCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTtgtg                                                                                                          >  6:184607/1‑148 (MQ=255)                                           ccTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTtgtg                                                                                                          >  4:429004/1‑147 (MQ=255)                                            cTTGATCCGGCCTGAGCGGGCGGAGGAGCTGGTCGACGTTTGCCTGTTACGGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTACTGGCGGCTGCAAGTGCAGCCAGTCGTGGGCGTAGCATGGTGAGTCCTTTTGTGCAg                                                                                                       <  8:372059/149‑1 (MQ=255)                                             ttGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGc                                                                                                      >  7:160949/1‑149 (MQ=255)                                                 tCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCgtg                                                                                                   >  4:2095/1‑148 (MQ=255)                                                   cTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTcc                                                                                                >  8:99999/1‑149 (MQ=255)                                                   cTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTcc                                                                                                >  2:227072/1‑149 (MQ=255)                                                      cATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                             <  3:429004/149‑1 (MQ=255)                                                        tGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGa                                                                                           <  4:11319/149‑1 (MQ=255)                                                          aGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAg                                                                                          <  3:418245/148‑1 (MQ=255)                                                          aGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGc                                                                                         <  1:298765/149‑1 (MQ=255)                                                           gcgggcggATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGc                                                                                         >  7:253249/1‑148 (MQ=255)                                                            cgggcggATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGa                                                                                       >  2:435355/1‑149 (MQ=255)                                                                 ggATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCa                                                                                  >  3:182774/1‑149 (MQ=255)                                                                   atgatCTGGTCGACCTTGGCCTGTAACTGGACGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCTAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACt                                                                                >  1:258628/1‑149 (MQ=255)                                                                        cTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATa                                                                                     <  8:125953/139‑1 (MQ=255)                                                                        cTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATa                                                                                     >  7:125953/1‑139 (MQ=255)                                                                         tGGTCGACCTTGGCCTGTAACTGCTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGCTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGt                                                                          >  3:344713/1‑149 (MQ=255)                                                                               aCCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGc                                                                    <  8:253249/149‑1 (MQ=255)                                                                                 cTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCTAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGt                                                                  <  2:258628/149‑1 (MQ=255)                                                                                  ttGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGa                                                                                                                                        <  7:198024/78‑1 (MQ=255)                                                                                  ttGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGa                                                                                                                                        >  8:198024/1‑78 (MQ=255)                                                                                   tGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGt                                                                  >  2:139494/1‑147 (MQ=255)                                                                                      ccTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGc                                                              >  8:141755/1‑148 (MQ=255)                                                                                        tGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGAc                                                           >  3:180890/1‑149 (MQ=255)                                                                                            aCTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACggg                                                        <  8:405139/148‑1 (MQ=255)                                                                                             cTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGa                                                      <  1:227072/149‑1 (MQ=255)                                                                                                gtcgtcGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAg                                                     <  5:397030/147‑1 (MQ=255)                                                                                                gtcgtcGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGc                                                    >  2:20666/1‑148 (MQ=255)                                                                                                   gtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCa                                                  >  3:114375/1‑147 (MQ=255)                                                                                                   gtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCa                                                  >  6:298927/1‑147 (MQ=255)                                                                                                   gtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACACGGGGTTCCGGCGTTAGCGACGGGGAGCCa                                                  >  6:95404/1‑147 (MQ=255)                                                                                                           gggCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTATCCGTTGAGCGATACCACTCGGGGTTCCCGCGTAAGCGACGGGGAGCCAAAGACCTTAt                                        >  6:30750/1‑149 (MQ=255)                                                                                                              cgcTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGt                                     >  3:157624/1‑149 (MQ=255)                                                                                                                    cAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGccc                               <  1:214538/149‑1 (MQ=255)                                                                                                                    cAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGccc                               <  8:160949/149‑1 (MQ=255)                                                                                                                         ccTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGCGGTTCCGGCGAAAGCGACGGGGAGCCAAAGACCTTATCGTCTGcgcaccc                           >  3:387936/1‑148 (MQ=255)                                                                                                                          cTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCccacccac                         >  3:289261/1‑149 (MQ=255)                                                                                                                             aGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCAt                      <  1:20666/149‑1 (MQ=255)                                                                                                                               gCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGa                    >  7:257484/1‑149 (MQ=255)                                                                                                                               gCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGa                    >  2:56890/1‑149 (MQ=255)                                                                                                                               gCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCAT‑GTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGa                    >  2:197192/1‑148 (MQ=255)                                                                                                                                cAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGAc                   <  7:99999/149‑1 (MQ=255)                                                                                                                                   gggTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAg                <  4:182774/149‑1 (MQ=255)                                                                                                                                     gTAATGGCGGGTGCAAGTGCAGCCAGAGGTGGGCGGTGCATGGTGAGTCCTTTTTTGGAGCGTGTCCGTTGAGCTATACCAGTCGGGGTTCCGGCGTTAGCGTCGTGGAGCCATAGACCTTATCGTCTGCCCAGCCACGATGTCAAGAt              <  4:344713/149‑1 (MQ=255)                                                                                                                                      tAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATg             <  3:2095/149‑1 (MQ=255)                                                                                                                                        aTGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATgg            <  4:180890/148‑1 (MQ=255)                                                                                                                                           ctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATc                                       >  6:197360/3‑118 (MQ=255)                                                                                                                                           ctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATc                                       <  5:197360/116‑1 (MQ=255)                                                                                                                                              gcTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGt     >  1:262088/1‑149 (MQ=255)                                                                                                                                                  cAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGtgat  >  4:200223/1‑148 (MQ=255)                                                                                                                                                  cAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGtgat  <  4:157624/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCGTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGTGAT  >  NC_002947/3277948‑3278240

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