Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,470,714	A→T	6.4%	I159N (ATT→AAT) ‡	PP_3961 ←	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,470,714	0	A	T	6.4%	209.3 / 3.0	72	I159N (ATT→AAT) ‡	PP_3961	transporter
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (2/3);  total (35/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 = 8.58e-01

TCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAGCGCGGCC  >  NC_002947/4470573‑4470860                                                                                                                                                |                                                                                                                                                   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                                                                                                                   <  4:260914/147‑1 (MQ=255)                             cTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCaa                                                                                                                   <  1:360805/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                                                                                                                     <  6:283006/140‑1 (MQ=255)                                  ggcTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTc                                                                                                                     >  5:283006/1‑140 (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                                                                                                                              <  1:102340/119‑1 (MQ=255)                                              aCCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGt                                                                                                                              >  2:102340/1‑119 (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                                                                                        >  8:184094/1‑149 (MQ=255)                                                      cGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCt                                                                                        >  4:240226/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                                                          <  7:184094/149‑1 (MQ=255)                                                                                      cacGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCgg                                                          <  1:282490/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                        <  5:21152/134‑1 (MQ=255)                                                                                                                                    tggcccGCCGATTTCGGTGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTggg                        >  6:21152/4‑137 (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)                                                                                                                                            cGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAGCGCgg    >  3:201617/1‑149 (MQ=255)                                                                                                                                              aaaTCGGCGGGCCAGCAATGAAGGTGACTTCGTCGAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgc                                                    <  1:378701/99‑1 (MQ=255)                                                                                                                                              aaaTCGGCGGGCCAGCAATGAAGGTGACTTCGTCGAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAgc                                                    >  2:378701/1‑99 (MQ=255)                                                                                                                                              aaaTCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGATCGCGGcc  >  1:277235/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   TCTTGCCCGGTGGCGGCATCGGCCTCTACTGGCGGCTCGGTGCCGACCTGCACCGCCAAGGCGAACACG‑‑CCGATGGCCAGGAACACGGCAGCAGCCAATGTACCGGCCTGCGGAAGTACCGCCACGCTCAGCCCCACCGAAATCGGCGGGCCAGCAATGAAGGTGACTTCGTCAAGAACGGTTTCCAGCGAGTAGGCTGTCTTCAGTTGCGCCTTGCCGCGGTAAATCGCGGTCCAGCGCGCGCGCACCATCGCCGACATGCTGGGCATGAAGCCGGCGAGCGCGGCC  >  NC_002947/4470573‑4470860

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