Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,626,459	C→G	6.6%	G348A (GGG→GCG)	PP_4944 ←	carbamoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,626,459	0	C	G	6.6%	185.6 / 4.4	61	G348A (GGG→GCG)	PP_4944	carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (2/2);  total (29/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 = 5.14e-01

CGCCACTTCGGCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAGGTAATGGTC  >  NC_002947/5626315‑5626604                                                                                                                                                  |                                                                                                                                                  cgcCACTTCGGCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGc                                                                                                                                               <  2:17954/149‑1 (MQ=255) cgcCACTTCGGCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGc                                                                                                                                               >  3:249300/1‑149 (MQ=255)    cACTTCGGCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTc                                                                                                                                            >  6:247303/1‑149 (MQ=255)           gCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCTGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAg                                                                                                                                     <  1:47367/149‑1 (MQ=255)           gCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAg                                                                                                                                     <  3:109728/149‑1 (MQ=255)                ttCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCgg                                                                                                                                <  1:56484/149‑1 (MQ=255)                  cGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGc                                                                                                                               <  3:332936/148‑1 (MQ=255)                     gTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGa                                                                                                                           >  7:237446/1‑149 (MQ=255)                          gggCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAAc                                                                                                                      >  8:20699/1‑149 (MQ=255)                          gggCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAAc                                                                                                                      <  3:255529/149‑1 (MQ=255)                            gcgcAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGTGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAgg                                                                                                                                    <  1:245904/134‑1 (MQ=255)                            gcgcAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAgg                                                                                                                                    >  2:245904/1‑133 (MQ=255)                                  gcggcgATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCttt                                                                                                              <  4:249585/149‑1 (MQ=255)                                   cggcgATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACGGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTc                                                                                                             <  8:237446/149‑1 (MQ=255)                                         tCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATcc                                                                                                       >  3:222688/1‑149 (MQ=255)                                         tCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATcc                                                                                                       <  4:40430/149‑1 (MQ=255)                                           aCGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCg                                                                                                      >  7:76695/1‑148 (MQ=255)                                                  cGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGc                                                                                               >  3:16413/1‑148 (MQ=255)                                                    ttGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCGCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCgat                                                                                            >  4:12964/1‑149 (MQ=255)                                                    ttGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCGCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCgat                                                                                            >  4:17933/1‑149 (MQ=255)                                                    ttGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCgat                                                                                            >  3:44890/1‑149 (MQ=255)                                                    ttGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCgat                                                                                            <  4:16413/149‑1 (MQ=255)                                                          tACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGAtctt                                                                                      >  7:97919/1‑149 (MQ=255)                                                            cgacgGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGAtcttct                                                                                    <  7:20699/149‑1 (MQ=255)                                                             gacgGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTg                                                                                   >  5:5513/1‑149 (MQ=255)                                                              acgGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTgg                                                                                  <  4:44890/149‑1 (MQ=255)                                                                  gCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCGCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGa                                                                              <  3:12964/149‑1 (MQ=255)                                                                  gCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCGCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGa                                                                              <  3:17933/149‑1 (MQ=255)                                                                  gCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGa                                                                              <  4:109645/149‑1 (MQ=255)                                                                  gCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCAGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGCTGATCTTCTGGtt                                                                                >  5:197926/1‑147 (MQ=255)                                                                      aGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCtt                                                                          <  8:97919/149‑1 (MQ=255)                                                                      aGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCtt                                                                          <  5:247303/149‑1 (MQ=255)                                                                        gTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGa                                                                        >  7:190432/1‑149 (MQ=255)                                                                              gTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCa                                                                  <  6:5513/149‑1 (MQ=255)                                                                                     aTCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGcaca                                                           <  8:76695/149‑1 (MQ=255)                                                                                      tcttctAGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGCCGTCGCCGGACGCCGGCTGGACGCACAGCTCTTTCCCATCCGGACGGGCGATGATCTTCTGGATGAGCATGACGTACAGCGCACAg                                                          >  8:283374/1‑149 (MQ=255)                                                                                           tcGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCa                                                                  >  5:90179/1‑136 (MQ=255)                                                                                           tcGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCa                                                                  <  6:90179/136‑1 (MQ=255)                                                                                           tcGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAgccgcc                                                     >  5:2074/1‑149 (MQ=255)                                                                                                 tgtttcccACGCCCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCg                                                 >  6:188785/7‑147 (MQ=255)                                                                                                 ggAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGaa                                               <  8:190432/149‑1 (MQ=255)                                                                                                       ccACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGCCGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAg                                         >  3:307407/1‑149 (MQ=255)                                                                                                        cACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGc                                        >  5:84985/1‑149 (MQ=255)                                                                                                        cACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGc                                        >  5:170683/1‑149 (MQ=255)                                                                                                         acgcccGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCt                                       <  5:188785/144‑1 (MQ=255)                                                                                                         aCGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCt                                       >  1:187454/1‑149 (MQ=255)                                                                                                                  ggAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGttt                              <  6:2074/149‑1 (MQ=255)                                                                                                                  ggAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGttt                              >  5:92024/1‑149 (MQ=255)                                                                                                                    aaaCATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTcg                            >  5:361727/1‑149 (MQ=255)                                                                                                                      aCATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGCAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCACCTTGCCGGTTTcgcg                          >  7:71208/1‑149 (MQ=255)                                                                                                                           g‑gCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAgg                      <  7:39112/148‑1 (MQ=255)                                                                                                                             gCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGGCGGCTTGACGAACAGCTCTTTCACATCGGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTGCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGgt                    <  6:197926/149‑1 (MQ=255)                                                                                                                             gCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGgt                    >  6:34627/1‑149 (MQ=255)                                                                                                                                    cgcCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGcc              >  4:19951/1‑148 (MQ=255)                                                                                                                                      ccAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAg           <  8:193199/149‑1 (MQ=255)                                                                                                                                       cAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAgg          >  5:16969/1‑149 (MQ=255)                                                                                                                                          cGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAGGTaa       <  6:170683/149‑1 (MQ=255)                                                                                                                                            gCAGTCCGGGCGTCGCCGGAGGCCGGCTGGAAGTACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTGAGCGCACAGCCGCCGGGGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGGCCAGGTAATg     <  8:71208/149‑1 (MQ=255)                                                                                                                                            gCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAGGTAATg     <  2:187454/149‑1 (MQ=255)                                                                                                                                              aGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAAGGTGTCGCCCAGGTAATGGt   <  4:222688/149‑1 (MQ=255)                                                                                                                                               gTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAGGTAATGGTc  <  6:16969/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CGCCACTTCGGCTGGTTCGGGTGACGGGCGCAGGCGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAG‑GCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTCGCCCAGGTAATGGTC  >  NC_002947/5626315‑5626604

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