Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,311,314	T→A	8.0%	F12I (TTT→ATT)	PP_3784 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,311,314	0	T	A	8.0%	173.8 / 3.5	59	F12I (TTT→ATT)	PP_3784	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base A (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAG  >  NC_002947/4311169‑4311462                                                                                                                                                   |                                                                                                                                                     acaACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTg                                                                                                                                                   >  6:324192/1‑149 (MQ=255)   aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCg                                                                                                                                                 >  5:36328/1‑149 (MQ=255)   aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCg                                                                                                                                                 >  8:164908/1‑149 (MQ=255)       gcgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGc                                                                                                                                             <  8:277939/149‑1 (MQ=255)         gTGCGCTCGCAAACCACTTCTTGCAATGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGccc                                                                                                                                           <  8:241373/149‑1 (MQ=255)             gcTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCg                                                                                                                                       >  7:209773/1‑149 (MQ=255)               tCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGcc                                                                                                                                     >  7:183641/1‑149 (MQ=255)                 gCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTg                                                                                                                                              >  2:268964/1‑138 (MQ=255)                 gCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTg                                                                                                                                              <  1:268964/138‑1 (MQ=255)                            gTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                         >  3:346431/1‑148 (MQ=255)                             tGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcatc                                                                                                                        >  1:273897/1‑148 (MQ=255)                                 aaGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAt                                                                                                                                                      <  2:341834/114‑1 (MQ=255)                                 aaGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAt                                                                                                                                                      >  1:341834/1‑114 (MQ=255)                                 aaGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                                    >  5:193295/1‑148 (MQ=255)                                   ggCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcat                                                                                                                 <  1:53954/149‑1 (MQ=255)                                     cGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                               <  6:36328/149‑1 (MQ=255)                                     cGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                               <  5:324192/149‑1 (MQ=255)                                      gTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                         >  4:224840/1‑138 (MQ=255)                                      gTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                         <  3:224840/138‑1 (MQ=255)                                       tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             <  7:164908/149‑1 (MQ=255)                                         tAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGc                                                                                                           >  4:164294/1‑149 (MQ=255)                                            ccGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCAc                                                                                                        <  6:193295/149‑1 (MQ=255)                                             cGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCt                                                                                                                                          <  8:212311/114‑1 (MQ=255)                                              gAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAATTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgt                                                                                                      >  7:291916/1‑149 (MQ=255)                                                      tAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCt                                                                                              <  4:346431/149‑1 (MQ=255)                                                        aCGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa                                                                                            <  3:164294/149‑1 (MQ=255)                                                                gCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACAtt                                                                                      <  2:273897/147‑1 (MQ=255)                                                                  ccccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGTGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAg                                                                                   >  1:232857/1‑149 (MQ=255)                                                                  ccccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGc                                                                                  >  6:191281/1‑149 (MQ=255)                                                                      gAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACt                                                                              <  8:209773/149‑1 (MQ=255)                                                                                 gcgcCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCaa                                                                   <  8:183641/149‑1 (MQ=255)                                                                                 gcgcCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCa                                                                    >  7:237187/1‑148 (MQ=255)                                                                                       tCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTg                                                             >  2:120547/1‑149 (MQ=255)                                                                                       tCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTg                                                             >  3:292903/1‑149 (MQ=255)                                                                                         ccGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgcc                                                           >  5:260344/1‑149 (MQ=255)                                                                                         ccGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAATTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgcc                                                           <  8:291916/149‑1 (MQ=255)                                                                                                tGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCa                                                                                    >  3:189747/1‑117 (MQ=255)                                                                                                tGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCa                                                                                    <  4:189747/117‑1 (MQ=255)                                                                                                       gAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAAc                                              <  1:120547/149‑1 (MQ=255)                                                                                                             catcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAg                                        >  8:8246/1‑149 (MQ=255)                                                                                                               tcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAgg                                       >  1:22595/1‑148 (MQ=255)                                                                                                               tcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGc                                      <  2:232857/149‑1 (MQ=255)                                                                                                                 aaTGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGcac                                    <  2:22595/149‑1 (MQ=255)                                                                                                                 aaTGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGca                                     <  8:237187/148‑1 (MQ=255)                                                                                                                  aTGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             <  4:195987/75‑1 (MQ=255)                                                                                                                  aTGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             >  3:195987/1‑75 (MQ=255)                                                                                                                       tAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTg                              >  8:98751/1‑149 (MQ=255)                                                                                                                       tAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTg                              >  7:126008/1‑149 (MQ=255)                                                                                                                         aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGTCTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGc                            <  5:191281/149‑1 (MQ=255)                                                                                                                           aGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGt                          <  6:260344/149‑1 (MQ=255)                                                                                                                             cATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgc                        >  4:315639/1‑149 (MQ=255)                                                                                                                                 cccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATAAGTGCCACGTGTAAGCCTGACACATTAAGCGACTGCCTGTCGCAATAATTGCCGCCATCCAGCAACGAACAGGCACCACATGGCGTGGTGCg                    >  2:331356/1‑149 (MQ=255)                                                                                                                                      cGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCtt                <  7:8246/148‑1 (MQ=255)                                                                                                                                      cGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTg               <  8:126008/149‑1 (MQ=255)                                                                                                                                           gCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTg          <  1:331356/149‑1 (MQ=255)                                                                                                                                           gCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTg          >  2:187348/1‑149 (MQ=255)                                                                                                                                             aaCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg         >  6:63237/1‑148 (MQ=255)                                                                                                                                             aaCGCATTTGCGTTGCCCTTCGCCGCCTCTCGTCTTCAATTCAACAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGAGTCTTGTCGTgg         >  4:32886/1‑148 (MQ=255)                                                                                                                                              aCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAg                                                  >  6:15310/1‑106 (MQ=255)                                                                                                                                                gCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGt     >  5:152479/1‑149 (MQ=255)                                                                                                                                                  aTTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGa   <  4:292903/149‑1 (MQ=255)                                                                                                                                                   tttGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAg  >  3:319879/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAG‑GAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAG  >  NC_002947/4311169‑4311462

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