Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,278,096	G→T	9.7%	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	9.7%	185.7 / 9.7	62	L9I (CTT→ATT) ‡	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base T (3/3);  total (30/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 = 7.56e-01

GTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGTGA  >  NC_002947/3277950‑3278239                                                                                                                                                   |                                                                                                                                                gTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCTTGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTg                                                                                                                                               <  5:16656/149‑1 (MQ=255)       gcgTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCACTTGCAGccg                                                                                                                                         <  2:272307/149‑2 (MQ=255)             aTAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTg                                                                                                                                   >  6:163638/1‑149 (MQ=255)             aTAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTg                                                                                                                                   >  8:162624/1‑149 (MQ=255)              tAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTgg                                                                                                                                  <  5:163638/149‑1 (MQ=255)                 aTGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCg                                                                                                                               <  7:162624/149‑1 (MQ=255)                    gCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGa                                                                                                                             <  6:222379/148‑1 (MQ=255)                     ccACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAg                                                                                                                            >  1:271260/1‑148 (MQ=255)                       acaGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCAt                                                                                                                         >  1:78620/1‑149 (MQ=255)                       acaGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCAt                                                                                                                         >  3:235026/1‑149 (MQ=255)                       acaGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCAt                                                                                                                         <  2:8506/149‑1 (MQ=255)                          gCCATGCCGGTGATGCCTTGTTCCTGCTTGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGt                                                                                                                      <  8:51084/149‑1 (MQ=255)                          gCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGt                                                                                                                      <  8:234965/149‑1 (MQ=255)                                    tGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCtttt                                                                                                            <  6:189268/149‑1 (MQ=255)                                      aTGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTtgt                                                                                                          <  6:243653/149‑1 (MQ=255)                                         ccTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTgg                                                                                                        <  3:57648/148‑1 (MQ=255)                                                  tGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCg                                                                                              >  4:203454/1‑149 (MQ=255)                                                    cATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCg                                                                                              <  5:69000/147‑1 (MQ=255)                                                    cATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                            >  2:273078/1‑149 (MQ=255)                                                    cATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                            <  4:215272/149‑1 (MQ=255)                                                     aTGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                            >  5:241206/1‑148 (MQ=255)                                                       gAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                            >  2:168529/1‑146 (MQ=255)                                                       gAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGtt                                                                                            <  1:168529/146‑1 (MQ=255)                                                             gcggATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATAc                                                                                   <  2:271260/149‑1 (MQ=255)                                                               ggATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCCGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCa                                                                                 >  3:13842/1‑149 (MQ=255)                                                                     tCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCggg                                                                           >  5:85993/1‑149 (MQ=255)                                                                     tCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCggg                                                                           >  1:162592/1‑149 (MQ=255)                                                                      cTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCgggg                                                                          >  5:44181/1‑149 (MQ=255)                                                                      cTGGTCGACCTTGGCCTGTAACTGGGCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGGGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCggg                                                                           <  6:241206/148‑1 (MQ=255)                                                                            gACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCgg                                                                    >  7:243714/1‑149 (MQ=255)                                                                              ccTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCg                                                                  >  2:159811/1‑149 (MQ=255)                                                                                  ggCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGtt                                                                <  6:44181/147‑1 (MQ=255)                                                                                     cTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGa                                                           >  1:277890/1‑149 (MQ=255)                                                                                        tAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACgg                                                        >  3:237081/1‑149 (MQ=255)                                                                                          aCTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACgggg                                                      <  2:78620/149‑1 (MQ=255)                                                                                          aCTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACggg                                                       <  4:237081/148‑1 (MQ=255)                                                                                              gtcgtcGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCgg                                                                            <  4:122104/123‑1 (MQ=255)                                                                                              gtcgtcGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCgg                                                                            >  3:122104/1‑123 (MQ=255)                                                                                               tcgtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGa                                                                                      >  6:71394/1‑112 (MQ=255)                                                                                               tcgtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGa                                                                                      <  5:71394/112‑1 (MQ=255)                                                                                               tcgtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACgg                                                        >  5:333236/1‑142 (MQ=255)                                                                                               tcgtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACgg                                                        <  6:333236/142‑1 (MQ=255)                                                                                                 gtcgGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCCGGGTTCCGGCGTTAGCGACGGGGAGCCa                                                 >  8:10066/1‑147 (MQ=255)                                                                                                    ggCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGAc                                            <  2:162592/149‑1 (MQ=255)                                                                                                          ggCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTAt                                       <  6:85993/148‑1 (MQ=255)                                                                                                          ggCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATc                                      <  2:277890/149‑1 (MQ=255)                                                                                                               tGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTcc                                                                      >  1:10052/1‑112 (MQ=255)                                                                                                               tGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTcc                                                                      <  2:10052/112‑1 (MQ=255)                                                                                                                 gCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTg                                 <  3:203454/147‑1 (MQ=255)                                                                                                                       ccTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGcccaccca                         >  1:175320/1‑149 (MQ=255)                                                                                                                              cAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGa                   >  8:206026/1‑148 (MQ=255)                                                                                                                               aaGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGAc                  >  3:97655/1‑148 (MQ=255)                                                                                                                                 gggTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGAc                  >  1:191357/1‑146 (MQ=255)                                                                                                                                 gggTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAg               <  1:159811/149‑1 (MQ=255)                                                                                                                                        tctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGAt             >  4:220411/4‑144 (MQ=255)                                                                                                                                        tctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGAt             <  3:220411/141‑1 (MQ=255)                                                                                                                                         gcggcTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATa        >  7:57099/1‑148 (MQ=255)                                                                                                                                         ctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCgg                                                                            <  4:314107/78‑1 (MQ=255)                                                                                                                                         ctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCgg                                                                            >  3:314107/3‑80 (MQ=255)                                                                                                                                         ctggcTGCACTTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATa        >  8:145239/3‑148 (MQ=255)                                                                                                                                             cTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGtg   <  7:206026/149‑1 (MQ=255)                                                                                                                                                cAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGGCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGtga  <  7:10066/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GTGGTCGCGTGCATAGATGGCCACAGCCATGCCGGTGATGCCTTGTTCCTGCATGAGCGGGCGGATGATCTGGTCGACCTTGGCCTGTAACTGGTCGTCGGCCTGGGCGCTGGCAGGGCCTAAGGCAAGGGTAATGGCGGCTGCAAGTGCAGCCAGACGTGGGCGGAGCATGGTGAGTCCTTTTGTGGAGCGTGTCCGTTGAGCGATACCACTCGGGGTTCCGGCGTTAGCGACGGGGAGCCAAAGACCTTATCGTCTGCCCACCCACCATGACAAGATGGATACAGTGA  >  NC_002947/3277950‑3278239

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