Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,703,478	A→T	6.4%	Q438L (CAG→CTG)	phaC‑II →	poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,703,478	0	A	T	6.4%	187.9 / 2.8	62	Q438L (CAG→CTG)	phaC‑II	poly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base T (2/2);  total (26/36)
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.97e-01

ACCGCCGGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATC  >  NC_002947/5703332‑5703620                                                                                                                                                   |                                                                                                                                                acCGCCGGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCTCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTgcagc                                                                                                                                               <  6:5276/149‑1 (MQ=255) acCGCCGGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCATCCCCGCCGCGCTTCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCAGGGCTTGAGGTGTGCTGTACGCCGAACGACATgcag                                                                                                                                                >  6:258375/1‑148 (MQ=255)     ccgGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGt                                                                                                                                           <  4:20369/149‑1 (MQ=255)     ccgGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGt                                                                                                                                           <  4:24926/149‑1 (MQ=255)          gTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTCCGGTACGCCGATCGACCTGCAGCAGGTCAata                                                                                                                                      >  3:106414/1‑149 (MQ=255)             cGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATc                                                                                                                                    <  4:106414/148‑1 (MQ=255)             cGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATc                                                                                                                                    <  3:61843/148‑1 (MQ=255)             cGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATc                                                                                                                                    >  7:46940/1‑148 (MQ=255)                cATCCTGTACTGGAATGCCGACAGCACGCGCCTTCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGa                                                                                                                                  >  8:153128/1‑147 (MQ=255)                  tCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGc                                                                                                                              <  5:146086/149‑1 (MQ=255)                         cTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTAcc                                                                                                                        <  8:46940/148‑1 (MQ=255)                          tGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGt                                                                                                                      >  5:123119/1‑149 (MQ=255)                            gAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGt                                                                                                                      <  6:263596/147‑1 (MQ=255)                                gCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCgg                                                                                                                >  5:60673/1‑149 (MQ=255)                                 ccGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCgg                                                                                                                >  6:176117/1‑148 (MQ=255)                                  cGACAGCACGCGCCTTCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGgca                                                                                                              <  7:153128/149‑1 (MQ=255)                                   gACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAat                                                                                                                                       >  6:283962/1‑123 (MQ=255)                                   gACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAat                                                                                                                                       <  5:283962/123‑1 (MQ=255)                                     cAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGgcagca                                                                                                           <  2:67641/149‑1 (MQ=255)                                       gCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAAc                                                                                                         <  7:88903/149‑1 (MQ=255)                                       gCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTAGCCGGCAGCAAc                                                                                                         >  1:59450/1‑149 (MQ=255)                                        cACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAAc                                                                                                         >  2:89612/1‑148 (MQ=255)                                         aCGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCTGCAGGTCGat                                                                                                                                       <  7:25652/117‑1 (MQ=255)                                         aCGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCTGCAGGTCGat                                                                                                                                       >  8:25652/1‑117 (MQ=255)                                         aCGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGa                                                                                                       >  4:31299/1‑149 (MQ=255)                                             gcCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGTAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACcac                                                                                                   <  5:258375/149‑1 (MQ=255)                                              ccTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACcgca                                                                                                  <  3:31299/149‑1 (MQ=255)                                                tGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACAt                                                                                                 >  3:280413/1‑148 (MQ=255)                                                  cccGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATcac                                                                                              <  6:123119/149‑1 (MQ=255)                                                        gcgcTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATg                                                                                        <  6:60673/149‑1 (MQ=255)                                                          gcTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATgg                                                                                       >  5:339948/1‑148 (MQ=255)                                                              cATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGc                                                                                 >  6:281239/1‑149 (MQ=255)                                                                 ggCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGgt                                                                              >  5:127754/1‑149 (MQ=255)                                                                    caCCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTTCGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCgg                                                                                                                <  8:311317/112‑1 (MQ=255)                                                                                  tttCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTtg                                                              <  6:339948/148‑1 (MQ=255)                                                                                            aaCCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTgg                                                   >  8:109140/1‑149 (MQ=255)                                                                                             aCCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGc                                                  <  5:176117/149‑1 (MQ=255)                                                                                               cccTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGc                                                  >  6:79041/1‑147 (MQ=255)                                                                                                ccTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAg                                               <  4:280413/149‑1 (MQ=255)                                                                                                    gACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAgcg                                             <  3:281733/147‑1 (MQ=255)                                                                                                       cTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTgg                                                   <  4:98552/138‑1 (MQ=255)                                                                                                       cTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTgg                                                   >  3:98552/1‑138 (MQ=255)                                                                                                          cgcgTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGa                                                >  8:290509/1‑138 (MQ=255)                                                                                                          cgcgTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGa                                                <  7:290509/138‑1 (MQ=255)                                                                                                          cgcgTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTAGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTc                                      <  2:59450/148‑1 (MQ=255)                                                                                                             gaccGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGt                                    <  1:123144/145‑1 (MQ=255)                                                                                                                cGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGGCGATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTg                                <  1:89612/149‑1 (MQ=255)                                                                                                                cGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTgg                               >  8:200239/1‑149 (MQ=255)                                                                                                                     tGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCg                                                                                >  3:284193/1‑95 (MQ=255)                                                                                                                     tGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCg                                                                                <  4:284193/95‑1 (MQ=255)                                                                                                                          gtgtGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCgg                     <  5:79041/149‑1 (MQ=255)                                                                                                                             tgCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGc                   <  6:127754/148‑1 (MQ=255)                                                                                                                              gCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATg                                                                                        >  3:259268/1‑79 (MQ=255)                                                                                                                              gCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATg                                                                                        <  4:259268/79‑1 (MQ=255)                                                                                                                                      ccGAGCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCTCATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAAGAGCGGGCATATCCagag         <  2:158060/149‑1 (MQ=255)                                                                                                                                        tattgACCTGCTGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGt                                    <  3:251124/116‑1 (MQ=255)                                                                                                                                        tattgACCTGCTGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGt                                    >  4:251124/5‑120 (MQ=255)                                                                                                                                          tCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGcatc     <  2:424/149‑1 (MQ=255)                                                                                                                                             aCCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTgg                                                       <  3:176662/96‑1 (MQ=255)                                                                                                                                             aCCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTgg                                                       >  4:176662/1‑96 (MQ=255)                                                                                                                                             aCCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCg                      <  4:71011/129‑1 (MQ=255)                                                                                                                                             aCCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCg                      >  3:71011/1‑129 (MQ=255)                                                                                                                                             aCCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCatcatc  <  5:281239/149‑1 (MQ=255)                                                                                                                                              ccTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACGATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTCGGTGGCGTGCGGCGCTTCGTGCTGGCAAACATCGGGCATCTCCAGAGCatcatc  >  6:220017/1‑148 (MQ=255)                                                                                                                                                  cagcagGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGc        <  2:87351/138‑1 (MQ=255)                                                                                                                                                  cagcagGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGc        >  1:87351/1‑138 (MQ=255)                                                                                                                                                   |                                                                                                                                                ACCGCCGGCGTTCGACATCCTGTACTGGAATGCCGACAGCACGCGCCTGCCCGCCGCGCTGCATGGCGACCTGCTGGAGTTTTTCAAGCTCAACCCCTTGACCTACGCGTCCGGGCTGGAGGTGTGCGGTACGCCGATCGACCTGCAGCAGGTCAATATCGACAGCTTTACCGTGGCCGGCAGCAACGACCACATCACACCATGG‑GATGCGGTGTACCGCTCGGCCTTGCTGCTGGGTGGCGAGCGGCGCTTCGTGCTGGCCAACAGCGGGCATATCCAGAGCATCATC  >  NC_002947/5703332‑5703620

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