Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,470,708	A→G	5.6%	V161A (GTG→GCG)	PP_3961 ←	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,470,708	0	A	G	5.6%	200.6 / 2.3	71	V161A (GTG→GCG)	PP_3961	transporter
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (2/2);  total (34/37)
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.76e-01

CACCGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAGC  >  NC_002947/4470560‑4470854                                                                                                                                                       |                                                                                                                                                   cacCGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCa                                                                                                                                                    <  2:163909/149‑1 (MQ=255)              tCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCTCCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCGCCGATTTcggtg                                                                                                                                        <  5:232259/148‑3 (MQ=255)              tCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggc                                                                                                                                         <  8:298090/147‑1 (MQ=255)                 tGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACTCGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGc                                                                                                                                     >  1:277717/1‑148 (MQ=255)                 tGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGc                                                                                                                                     >  1:272428/1‑148 (MQ=255)                   cccGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCa                                                                                                                                   <  6:69229/148‑1 (MQ=255)                             gcATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGt                                                                                                                        <  1:146615/149‑1 (MQ=255)                                   gCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTc                                                                                                                  <  3:118488/149‑1 (MQ=255)                                          cTGGCTGCTCGGTGCCGACGTGCACGGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACAGGGCAGCAGCCAATGTACCGGCCTTCGGAAGTACCGGCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCa                                                                                                              <  6:40078/146‑1 (MQ=255)                                          cTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACGCACCTTTGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCaa                                                                                                             <  1:131295/149‑1 (MQ=255)                                          cTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCaa                                                                                                             <  1:360805/147‑1 (MQ=255)                                          cTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCaa                                                                                                             <  4:260914/147‑1 (MQ=255)                                           tGGCGGCTCGGTGCCGACCTGCACCGCCAAAGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGaa                                                                                                          >  1:337036/1‑149 (MQ=255)                                             gcggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggc                                                                                                                                         <  4:334573/116‑1 (MQ=255)                                             gcggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATcggc                                                                                                                                         >  3:334573/1‑116 (MQ=255)                                             gcggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGaa                                                                                                          <  4:55473/147‑1 (MQ=255)                                               ggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTc                                                                                                               >  5:283006/1‑140 (MQ=255)                                               ggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTc                                                                                                               <  6:283006/140‑1 (MQ=255)                                                   cGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTc                                                                                                   >  8:158443/1‑148 (MQ=255)                                                     gTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAg                                                                                                >  4:379047/1‑149 (MQ=255)                                                      tGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                               >  5:114268/1‑149 (MQ=255)                                                        ccGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGc                                                                                               >  6:279214/1‑147 (MQ=255)                                                          gACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAg                                                                                            >  2:210706/1‑148 (MQ=255)                                                           aCCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGt                                                                                                                        >  2:102340/1‑119 (MQ=255)                                                           aCCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGt                                                                                                                        <  1:102340/119‑1 (MQ=255)                                                               gCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCt                                                                                      >  7:196949/1‑149 (MQ=255)                                                                cACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCt                                                                                      <  6:114268/148‑1 (MQ=255)                                                                cACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTg                                                                                     >  5:373021/1‑149 (MQ=255)                                                                 aCCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGt                                                                                    >  1:207309/1‑149 (MQ=255)                                                                  ccGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTc                                                                                   <  2:272428/149‑1 (MQ=255)                                                                  ccGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTc                                                                                   <  6:25001/149‑1 (MQ=255)                                                                  ccGCCAACGCTAACACC‑‑CCGATGTCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGCTTTCCAGCGAGTAGGCTGTc                                                                                   <  2:277717/149‑1 (MQ=255)                                                                   cGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCt                                                                                  >  4:32788/1‑149 (MQ=255)                                                                   cGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCt                                                                                  >  4:240226/1‑149 (MQ=255)                                                                   cGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCt                                                                                  >  8:184094/1‑149 (MQ=255)                                                                       aaGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAg                                                                              <  6:373021/149‑1 (MQ=255)                                                                         ggCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGtt                                                                            >  6:361871/1‑149 (MQ=255)                                                                           cGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTgc                                                                          <  8:196949/149‑1 (MQ=255)                                                                            gAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTgcg                                                                         >  5:22808/1‑149 (MQ=255)                                                                            gAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTgcg                                                                         <  5:279214/149‑1 (MQ=255)                                                                                acG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCtt                                                                     >  5:232329/1‑149 (MQ=255)                                                                                  g‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGc                                                                   >  1:190900/1‑149 (MQ=255)                                                                                        aTGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCcgcg                                                               <  7:158443/149‑1 (MQ=255)                                                                                           gCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGt                                                             <  7:323831/148‑1 (MQ=255)                                                                                           gCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTa                                                            <  5:184948/149‑1 (MQ=255)                                                                                             cAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTaaa                                                          >  1:224948/1‑149 (MQ=255)                                                                                              aGGAACACGGCAGCAGCCAATGTACCGGCCTGCAGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAAt                                                         >  8:55029/1‑149 (MQ=255)                                                                                                  acacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATcgcg                                                     <  5:137529/149‑1 (MQ=255)                                                                                                  acacGGCAGCAGCCAATGTACCGGCCTGCAGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATcgcg                                                     <  7:55029/149‑1 (MQ=255)                                                                                                   cacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCgg                                                    <  1:282490/149‑1 (MQ=255)                                                                                                   cacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCgg                                                    <  7:184094/149‑1 (MQ=255)                                                                                                         agcagcCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgc                                              >  7:68338/1‑149 (MQ=255)                                                                                                               cAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgcgcgcgc                                        <  5:231668/149‑1 (MQ=255)                                                                                                                     aCCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCa                                  >  3:265402/1‑149 (MQ=255)                                                                                                                      ccGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCAt                                 <  4:265402/149‑1 (MQ=255)                                                                                                                       cGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATc                                <  1:210706/149‑1 (MQ=255)                                                                                                                                gAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATg                       >  8:84902/1‑149 (MQ=255)                                                                                                                                gAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATg                       >  6:51466/1‑149 (MQ=255)                                                                                                                                    tACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTgg                   >  8:132457/1‑149 (MQ=255)                                                                                                                                      ccGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTggg                  <  3:379047/148‑1 (MQ=255)                                                                                                                                       cGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATcgcg                                                     >  6:281222/1‑112 (MQ=255)                                                                                                                                       cGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATcgcg                                                     <  5:281222/112‑1 (MQ=255)                                                                                                                                         ccACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCAt               >  6:312253/1‑148 (MQ=255)                                                                                                                                         ccACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCAt               >  7:202086/1‑148 (MQ=255)                                                                                                                                               tCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCg        >  1:330326/1‑149 (MQ=255)                                                                                                                                                 tggcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTggg                  >  6:21152/4‑137 (MQ=255)                                                                                                                                                 tggcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTggg                  <  5:21152/134‑1 (MQ=255)                                                                                                                                                 tggcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGc      >  5:389466/4‑149 (MQ=255)                                                                                                                                                 aGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGc      <  1:226732/149‑1 (MQ=255)                                                                                                                                                   ccccACCGAATTCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGa    <  1:133368/149‑1 (MQ=255)                                                                                                                                                    cccACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAg   <  2:337036/149‑1 (MQ=255)                                                                                                                                                     ccACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAgc  <  2:190900/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CACCGACATGCTGTCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAGC  >  NC_002947/4470560‑4470854

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