Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,111,188	G→C	11.9%	P108A (CCA→GCA)	PP_5361 ←	protein P47K

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,111,188	0	G	C	11.9%	154.4 / 13.5	59	P108A (CCA→GCA)	PP_5361	protein P47K
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base C (3/4);  total (28/31)
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.32e-01

GTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111040‑6111336                                                                                                                                                       |                                                                                                                                                     gTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTgg                                                                                                                                                      <  4:23634/149‑1 (MQ=255)      gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTg                                                                                                                                                 <  5:182745/149‑1 (MQ=255)      gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTg                                                                                                                                                 >  2:154113/1‑149 (MQ=255)           ggCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGGCCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGCCAGTGGCtga                                                                                                                                            <  4:290921/149‑3 (MQ=255)             cTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCaga                                                                                                                                          <  7:106601/149‑1 (MQ=255)               ggCCAGGCTGGCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCTGAGCCACTGCCAGTGGCTgaga                                                                                                                                          <  6:287374/149‑1 (MQ=38)                  cAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                                     <  4:153961/149‑1 (MQ=255)                  cAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                                     >  4:57431/1‑149 (MQ=255)                    ggCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCg                                                                                                                                   <  1:154113/149‑1 (MQ=255)                         tCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTg                                                                                                                                <  2:295683/147‑1 (MQ=255)                         tCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGa                                                                                                                              >  3:179490/1‑149 (MQ=255)                                ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGCCCATGGTGTCCAGTCGCGCCATATCGGACAGCCTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑CGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGAt                                                                                                                        >  1:315788/1‑148 (MQ=255)                                ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGCGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGa                                                                                                                         <  6:196401/147‑1 (MQ=255)                                 cTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGAt                                                                                                                        >  1:273418/1‑147 (MQ=255)                                   gATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCaga                                                                                                                                          >  2:231464/1‑127 (MQ=255)                                   gATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCaga                                                                                                                                          <  1:231464/127‑1 (MQ=255)                                                gAAGTTCAAGCCATCGCCCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGCCAGCGCCCCTGCTCGTCGCGGACAGTGAAGGTCTC‑‑AGCCCCTGGCAGTGGCTCAGCGAGACCGGTGGCCGCGATCAGCAAAGAATCGaacc                                                                                                       >  2:169709/1‑147 (MQ=255)                                                gAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGCCAGTGGCTGag                                                                                                                                           <  7:158247/113‑1 (MQ=255)                                                gAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGCCAGTGGCTGag                                                                                                                                           >  8:158247/1‑113 (MQ=255)                                                    ttCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGccc                                                                                                   <  3:57431/149‑1 (MQ=255)                                                     tCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGccc                                                                                                   <  8:164240/148‑1 (MQ=255)                                                       aaGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGAATGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCctc                                                                                                >  6:92293/1‑149 (MQ=255)                                                         gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGc                                                                                                     <  4:301616/142‑1 (MQ=255)                                                         gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGc                                                                                                     >  3:301616/1‑142 (MQ=255)                                                         gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCtt                                                                                              >  2:219113/1‑149 (MQ=255)                                                           cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCg                                                                                            >  3:3510/1‑149 (MQ=255)                                                            aTCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCgg                                                                                           <  2:273418/149‑1 (MQ=255)                                                             tCGAACACGGTGACCATGGTGTCCAGTCGCGCAATATCGGACAGAATGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGATAAGAGATACCGGTGGACTCGATCAGCAAATAATCg                                                                                                           >  5:29091/1‑132 (MQ=255)                                                                  cACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCTGTGGCTCAGAGATACCGGTTGACTCGATCAGCAAATAATCGAAACGc                                                                                                     <  3:16335/133‑1 (MQ=255)                                                                  cACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATAAGCAAATAATCGAAACga                                                                                                     >  4:16335/1‑132 (MQ=255)                                                                               gtgtCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACctcctc                                                                        <  8:120559/149‑1 (MQ=255)                                                                                    cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCa                                                                                        >  1:168159/1‑128 (MQ=255)                                                                                    cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCa                                                                                        <  2:168159/128‑1 (MQ=255)                                                                                    cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAg                                                                     >  4:307983/1‑147 (MQ=255)                                                                                    cAGGCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGG‑CAGTCGCGCAACCTCCTCGAGCAg                                                                  <  8:42692/149‑1 (MQ=255)                                                                                     aGTCGCGCCATATCGGCCAGACGGCGACCCTGCGCGGCGCGGAAAGTGAAGGGCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAATCGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAg                                                                  <  1:169709/149‑1 (MQ=255)                                                                                            ccATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGc                                                                                                                    >  1:7184/1‑92 (MQ=255)                                                                                            ccATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGc                                                                                                                    <  2:7184/92‑1 (MQ=255)                                                                                                   ggACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAAc                                                     <  8:197079/148‑1 (MQ=255)                                                                                                    gACAGACTGCGACCCTGATCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATAAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGAGCAACCTCCACGAGCAGGTCTTCGCGCAAcct                                                   >  2:58376/1‑147 (MQ=255)                                                                                                      cAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgc                                                 <  3:307983/149‑1 (MQ=255)                                                                                                      cAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgc                                                 <  5:92293/149‑1 (MQ=255)                                                                                                      cAGAATGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑ATCCACTGGCAGTGGCTCAGAGATACCGTTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgc                                                 <  1:58376/149‑1 (MQ=255)                                                                                                        gACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgcag                                               <  4:3510/149‑1 (MQ=255)                                                                                                        gACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgc                                                 >  6:201983/1‑147 (MQ=255)                                                                                                        gACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgc                                                 <  5:201983/147‑1 (MQ=255)                                                                                                             cGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCACCCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAAt                                          >  6:103310/1‑149 (MQ=255)                                                                                                              gACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATg                                         <  2:9221/149‑1 (MQ=255)                                                                                                                cccTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATg                                         >  1:83457/1‑147 (MQ=255)                                                                                                                        gtcgCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGc                               >  8:260419/1‑149 (MQ=255)                                                                                                                                   ggtatctctg‑‑AGCCACTGCCAGTGGCTGAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGAc                     >  3:161875/7‑148 (MQ=255)                                                                                                                                         tctg‑‑AGCCACTGCCAGTGGCTGGGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtc              >  5:115087/1‑149 (MQ=255)                                                                                                                                           tc‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGCCCAGTTtctc            >  6:201298/1‑149 (MQ=255)                                                                                                                                           tc‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtctc            >  6:79620/1‑149 (MQ=255)                                                                                                                                               aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtctc            <  7:260419/147‑1 (MQ=255)                                                                                                                                                 ccACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa        <  1:219113/149‑1 (MQ=255)                                                                                                                                                 ccACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGGCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa        <  2:315788/149‑1 (MQ=255)                                                                                                                                                 ccACTGCCAGTGGCTGGGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa        <  6:115087/149‑1 (MQ=255)                                                                                                                                                   aCTGGCAGTGGCTCAGAGATACCGGTGGACTCGCTCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAgc      >  1:32166/1‑149 (MQ=255)                                                                                                                                                     tGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAgcgc    >  2:329503/1‑149 (MQ=255)                                                                                                                                                       gCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCgg  >  2:164956/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTC‑‑AGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111040‑6111336

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