Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,860,700:1	+C	8.3%	coding (237/654 nt)	purN ←	phosphoribosylglycinamide formyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,860,700	1	.	C	8.3%	232.8 / 14.2	60	coding (237/654 nt)	purN	phosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base . (26/29);  new base C (3/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CCCGGCGATGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAC  >  NC_002947/1860552‑1860840                                                                                                                                                      |                                                                                                                                              cccGGCGATGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATg                                                                                                                                                <  8:144523/149‑1 (MQ=255)  ccGGCGATGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATg‑c                                                                                                                                              <  1:261591/149‑1 (MQ=255)  ccGGCGATGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATg‑c                                                                                                                                              >  7:212532/1‑149 (MQ=255)         tGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCAt                                                                                                                                       <  7:238743/149‑1 (MQ=255)            gTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑aa                                                                                                                                   <  8:287800/149‑2 (MQ=255)                gCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                               >  6:56874/1‑149 (MQ=255)                 cAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                               <  5:224240/148‑1 (MQ=255)                     cccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCCGCGCTAAT‑GGcgcgcatc                                                                                                                           <  1:64689/149‑5 (MQ=255)                      ccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCg                                                                                                                         <  8:212532/149‑1 (MQ=255)                       cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATc                                                                                                                          <  8:136835/147‑1 (MQ=255)                       cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGa                                                                                                                        >  7:136835/1‑149 (MQ=255)                          gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGaa                                                                                                                       <  8:83323/147‑1 (MQ=255)                               tGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCa                                                                                                                >  8:129761/1‑149 (MQ=255)                                  gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACgg                                                                                                             <  4:278008/149‑1 (MQ=255)                                      cagcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGt                                                                                                         <  1:82154/149‑1 (MQ=255)                                         cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCg                                                                                                       >  1:175118/1‑148 (MQ=255)                                           aCGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAAc                                                                                                    <  6:276291/149‑1 (MQ=255)                                              ggTGCAGGTTCAGCAGGCGGCCCTGATAGTGGCCCACGAAGCCGCCACTGAGGAGGCGCATCAAACCGCCCAGCACCACCAGGTCCGGGGAAAAGCCGTCTATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                 <  2:175118/149‑1 (MQ=255)                                               gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCttt                                                                                                                 <  5:64671/132‑2 (MQ=37)                                               gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCttt                                                                                                                 >  6:64671/1‑131 (MQ=37)                                                   aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                            >  6:339426/1‑149 (MQ=255)                                                               cggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                       <  5:8035/142‑1 (MQ=255)                                                               cggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGa                                                                                       >  6:8035/1‑142 (MQ=255)                                                                      gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                               >  6:45909/1‑95 (MQ=255)                                                                      gATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGcgcg                                                                                                                               <  5:45909/95‑1 (MQ=255)                                                                       aTAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                          <  8:23111/147‑1 (MQ=255)                                                                          gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                         >  3:94580/1‑129 (MQ=255)                                                                          gTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                         <  4:94580/129‑1 (MQ=255)                                                                                   gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                            >  4:297269/1‑149 (MQ=255)                                                                                   gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                            >  4:258618/1‑149 (MQ=255)                                                                                        cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                 <  7:227977/107‑1 (MQ=255)                                                                                        cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                 >  8:227977/1‑107 (MQ=255)                                                                                        cgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggcg                                                       >  4:25371/1‑149 (MQ=255)                                                                                          ccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGcgcg                                                     <  6:148733/149‑1 (MQ=255)                                                                                             cTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTg                                                  >  7:63592/1‑149 (MQ=255)                                                                                                aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                <  7:129761/148‑1 (MQ=255)                                                                                                  gATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGc                                             >  3:236325/1‑149 (MQ=255)                                                                                                    tGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgtgt                                                                                             <  4:268229/99‑1 (MQ=255)                                                                                                    tGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgtgt                                                                                             >  3:268229/1‑99 (MQ=255)                                                                                                       gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATTGAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCg                                           >  1:277678/1‑147 (MQ=255)                                                                                                       gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                        <  5:56874/149‑1 (MQ=255)                                                                                                       gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGCTGTGGTCGAGCACGACCCTGTCGATTCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                        >  3:248594/1‑149 (MQ=255)                                                                                                       gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATAGAAGGCTTCAAGGCCGTCGAACTGGGTGTGGTCGAGCACCACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGGCGGAg                                        >  2:173536/1‑149 (MQ=255)                                                                                                           gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGAGCGGCATCGAACGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                    >  5:61848/1‑149 (MQ=255)                                                                                                             aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCATGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                  >  8:37886/1‑149 (MQ=255)                                                                                                             aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCATGCCGTAGGCATCg                                  >  4:310290/1‑149 (MQ=255)                                                                                                             aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCg                                  <  4:236325/149‑1 (MQ=255)                                                                                                                     aGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                          <  5:339426/149‑1 (MQ=255)                                                                                                                      gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                    >  6:148792/1‑138 (MQ=255)                                                                                                                      gCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                    <  5:148792/138‑1 (MQ=255)                                                                                                                       caccacCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                         >  7:156848/1‑148 (MQ=255)                                                                                                                         ccaccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCCTT‑AGCGCGGGATGGAAGGCTTGACGGGCGGCGAACGGGGGGGGGTCGAGGACGACGCTGTCGCTGCCGGCCGCTGCGGGGGGGTGGAGGCCGGAGGCATCGGCGCGGTTggc                       >  2:170975/1‑147 (MQ=255)                                                                                                                          caccaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGa                       >  7:99926/1‑147 (MQ=255)                                                                                                                           accaGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCAAT‑AGCGCGCCATCCAAGGCTAAACGGCCGGCCAACAGCGAGTGGACGAGCACGAAGCACTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAAAGGCCCGCTTGGAAAg                    >  2:124088/1‑149 (MQ=255)                                                                                                                                gTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcg               <  1:173536/149‑1 (MQ=255)                                                                                                                                     gggCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg          <  8:156848/149‑1 (MQ=255)                                                                                                                                     gggCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg          <  8:99926/149‑1 (MQ=255)                                                                                                                                       gCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGAt        <  3:258618/149‑1 (MQ=255)                                                                                                                                         aaaGCCGTCGATG‑CGCGCCATA‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATAc      <  1:170975/149‑1 (MQ=255)                                                                                                                                           aGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGc    <  4:248594/149‑1 (MQ=255)                                                                                                                                           aGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGc    <  3:25371/149‑1 (MQ=255)                                                                                                                                            gCCTTCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCa   >  3:22125/1‑149 (MQ=38)                                                                                                                                             ccttCGATGCCGCGCTAAT‑GGCGC‑GCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc  >  7:255816/4‑149 (MQ=38)                                                                                                                                                      |                                                                                                                                              CCCGGCGATGCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATG‑CGCGCCATT‑AGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAC  >  NC_002947/1860552‑1860840

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