Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	399,528	A→G	6.4%	V19A (GTG→GCG)	PP_0332 ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	399,528	0	A	G	6.4%	181.5 / 2.6	62	V19A (GTG→GCG)	PP_0332	acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

TAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACA  >  NC_002947/399382‑399674                                                                                                                                                   |                                                                                                                                                   tAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACg                                                                                                                                                  <  6:54987/149‑1 (MQ=255)       tGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTg                                                                                                                                            <  2:268049/149‑1 (MQ=255)          acaGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatc                                                                                                                                         <  1:202829/149‑3 (MQ=255)          acaGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatc                                                                                                                                         <  6:341480/149‑3 (MQ=255)                    gATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAggg                                                                                                                               <  6:257750/149‑1 (MQ=255)                       ggCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGggtgg                                                                                                                            >  7:34677/1‑149 (MQ=255)                        gCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGggtgg                                                                                                                            <  3:368207/148‑1 (MQ=255)                         cTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTGGAGCAGGGCGAAGACGCGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGggtgg                                                                                                                            >  4:254545/1‑147 (MQ=255)                              tCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGTCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTga                                                                                                                      <  1:71944/148‑1 (MQ=255)                              tCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTga                                                                                                                      <  6:208088/148‑1 (MQ=255)                              tCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTga                                                                                                                      <  5:300654/148‑1 (MQ=255)                                gATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGaa                                                                                                                   >  7:46702/1‑149 (MQ=255)                                gATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGaa                                                                                                                   >  5:116945/1‑149 (MQ=255)                                    tCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGc                                                                                                                <  1:309786/148‑1 (MQ=255)                                    tCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGc                                                                                                                <  6:301542/148‑1 (MQ=255)                                     ccGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGcc                                                                                                               <  8:34677/148‑1 (MQ=255)                                      cGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCgg                                                                                                             >  7:26818/1‑149 (MQ=255)                                            ggTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTcag                                                                                                         >  8:304746/1‑147 (MQ=255)                                                tCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGc                                                                                                    <  8:26818/148‑1 (MQ=255)                                                       cccacATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCAt                                                                                             >  5:225562/1‑148 (MQ=255)                                                        ccacATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATa                                                                                            >  4:2237/1‑148 (MQ=255)                                                         cacaTCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATaaa                                                                                          <  6:225562/149‑1 (MQ=255)                                                             tCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAgaga                                                                                      >  6:135443/1‑149 (MQ=255)                                                              cACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGAc                                                                                     >  3:290812/1‑149 (MQ=255)                                                               aCCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTTAGCAGGAGCAGGGTGGTAGTGAGAATGcc                                                                                                               <  4:78013/122‑1 (MQ=255)                                                               aCCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGcc                                                                                                               >  3:78013/1‑122 (MQ=255)                                                               aCCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGAct                                                                                    <  4:290812/149‑1 (MQ=255)                                                                ccGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGActc                                                                                   <  6:116945/149‑1 (MQ=255)                                                                   cccACGACGCGTAGTCGCGGCCACGGCCTTGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttg                                                                                >  3:74827/1‑149 (MQ=255)                                                                   cccACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttg                                                                                >  5:313740/1‑149 (MQ=255)                                                                         acgCGTCGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTgg                                                                          <  6:268294/149‑1 (MQ=255)                                                                          cgcgTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTggg                                                                         <  5:343848/149‑1 (MQ=255)                                                                           gcgTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGga                                                                        <  3:89767/149‑1 (MQ=255)                                                                              tAGTCGCGGCCACGGCCTTGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                     <  4:74827/149‑1 (MQ=255)                                                                              tAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                     <  6:343411/149‑1 (MQ=255)                                                                              tAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                     <  7:304746/149‑1 (MQ=255)                                                                               aGTCGCGGCCACGGCCTGGCACCACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAATAGACTATTGTTGTGGGAGAt                                                                     >  4:338313/1‑148 (MQ=255)                                                                                   gcgGCCACGGCCTGGCAACACCAGCTTTAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                     <  3:59170/144‑1 (MQ=255)                                                                                   gcgGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                                     >  4:59170/1‑144 (MQ=255)                                                                                                ggCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGAt                                                    >  4:98556/1‑148 (MQ=255)                                                                                                   aaCACCAGCTTGCGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGc                                                <  6:313740/149‑1 (MQ=255)                                                                                                   aaCACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAg                                                 >  8:20805/1‑148 (MQ=255)                                                                                                            ttGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTg                                       >  7:137164/1‑149 (MQ=255)                                                                                                            ttGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACAACGGTATTGAGCAGCAGCACGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAATCGGTTGGAAATGATAACAAGGGCGGTg                                       >  5:166840/1‑149 (MQ=255)                                                                                                              gtgCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTTTTGAGCAGCAGCAGGGGGGTAGGGAGTATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTg                                       <  3:2237/145‑1 (MQ=255)                                                                                                                gCAGGGCGAAGACCAGCAGCTGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCt                                   <  8:4076/149‑1 (MQ=255)                                                                                                                gCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGcc                                    >  3:333991/1‑148 (MQ=255)                                                                                                                      cGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGcc                              >  8:312099/1‑148 (MQ=255)                                                                                                                      cGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGcc                              >  7:101974/1‑148 (MQ=255)                                                                                                                       gAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGCTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCa                             >  4:80305/1‑148 (MQ=255)                                                                                                                             cagcagCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTc                       >  3:55170/1‑148 (MQ=255)                                                                                                                             cagcagCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTc                       <  4:55170/148‑1 (MQ=255)                                                                                                                               gcagcGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCaa                    >  4:108523/1‑149 (MQ=255)                                                                                                                                cagcGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttgtt                                                                              <  3:349999/90‑1 (MQ=255)                                                                                                                                cagcGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttgtt                                                                              >  4:349999/1‑90 (MQ=255)                                                                                                                                        tgctCAATACCGTGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAcg           >  3:194713/4‑149 (MQ=255)                                                                                                                                        tgctCAATACCGTGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTAAGCCAAATGTCCAATGCCCGAcg           >  8:101418/4‑149 (MQ=255)                                                                                                                                         gATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAcgg          >  8:332166/1‑149 (MQ=255)                                                                                                                                           tCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAcggcg        <  5:135443/149‑1 (MQ=255)                                                                                                                                            cAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGGGGGAGATAAAAGGTTGGCAAGGATAAGCAGGGCGGTGGCCTACGCCAATTGTCCAAGGCCCGAcggcgg       >  4:47616/1‑149 (MQ=255)                                                                                                                                                  cACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACa  >  4:152115/1‑148 (MQ=255)                                                                                                                                                  cACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCCGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACa  >  5:318761/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TAGGGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGGCGGGGACA  >  NC_002947/399382‑399674

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