Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,457,146	G→A	6.0%	Y111Y (TAC→TAT)	PP_1275 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,457,146	0	G	A	6.0%	252.1 / 3.6	83	Y111Y (TAC→TAT)	PP_1275	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/35);  new base A (2/3);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ATGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTT  >  NC_002947/1457000‑1457291                                                                                                                                                   |                                                                                                                                                  aTGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTTCTCGGCTCGGAATTCAAGGGTCTGCGGCCGATATCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCACGGCTGCAATGCACCTTGCACGCGGTa                                                                                                                                                 >  7:319283/1‑149 (MQ=255) aTGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTa                                                                                                                                                 >  6:89952/1‑149 (MQ=255)   ggaggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGATacc                                                                                                                                               >  5:282386/1‑147 (MQ=255)     aggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGCTCTGCGGCCGATAGCCCGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCCGCTCCAATACACCTGGCCCGCGGAATcgc                                                                                                                                             >  5:335978/1‑149 (MQ=255)      ggaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATcgcg                                                                                                                                            <  4:179127/149‑1 (MQ=255)       gaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCgg                                                                                                                                           >  6:221427/1‑149 (MQ=255)         aaCACATCAAGTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCCGTGGCGGGCGGTCTGGCCACGCGGCCGCCGGCGGCAATTCCCCTGGCCCGCGGTAGCGCgggg                                                                                                                                         >  2:169209/1‑149 (MQ=255)           aacaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGATACCGCgggcca                                                                                                                                       <  6:282386/148‑4 (MQ=255)            acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGAAAGCCGGAACGCACCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCga                                                                                                                                      >  3:8717/1‑148 (MQ=255)             caTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCCCCGCGATACCGCgggccagg                                                                                                                                     <  4:391841/149‑6 (MQ=255)              aTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTg                                                                                                                                    >  6:28881/1‑149 (MQ=255)              aTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCCCCGCGATACCGCgggccagg                                                                                                                                     >  2:33949/1‑143 (MQ=255)              aTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCCCCGCGATACCGCgggccagg                                                                                                                                     <  1:33949/148‑6 (MQ=255)                 aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAg                                                                                                                                 <  5:142112/149‑1 (MQ=255)                 aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAg                                                                                                                                 >  6:109537/1‑149 (MQ=255)                 aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAg                                                                                                                                 >  5:50416/1‑149 (MQ=255)                 aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAg                                                                                                                                 <  5:283725/149‑1 (MQ=255)                      gCAAACGTGCCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGgtgta                                                                                                                            >  6:12608/1‑148 (MQ=255)                      gCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                            >  3:306618/1‑149 (MQ=255)                        aaaCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTtct                                                                                                                          <  7:210773/149‑1 (MQ=255)                        aaaCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTtct                                                                                                                          >  5:73418/1‑149 (MQ=255)                            gTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGc                                                                                                                      <  8:368990/149‑1 (MQ=255)                                 ttGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGGCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCg                                                                                                                                       <  8:375666/127‑1 (MQ=255)                                 ttGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCg                                                                                                                                       >  7:375666/1‑127 (MQ=255)                                    cAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAgg                                                                                                                  <  3:193683/145‑1 (MQ=255)                                    cAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAgg                                                                                                                  >  4:193683/1‑145 (MQ=255)                                     aGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCAc                                                                                                             <  5:221427/149‑1 (MQ=255)                                      ggTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGGCGGAACGCAGCGGCGGGAATGGCAGGCAGTCTTGCCAGTCGGCCCCCGGCTGCAATTCACCTGGCCCGCTGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTGACg                                                                                                            <  6:335978/149‑1 (MQ=255)                                         gctcggctcgGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTa                                                                                                         >  3:377135/1‑149 (MQ=255)                                          ctcggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTAt                                                                                                        <  5:28881/149‑1 (MQ=255)                                           tcggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATc                                                                                                       >  1:316701/1‑149 (MQ=255)                                            cggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCACCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCg                                                                                                      >  6:206948/1‑149 (MQ=255)                                                  ggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCTGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGt                                                                                                <  4:359270/149‑1 (MQ=255)                                                  ggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCgg                                                                                                 <  4:306618/148‑1 (MQ=255)                                                         cAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTAcg                                                                                         <  3:109853/149‑1 (MQ=255)                                                             gTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTAcgccgc                                                                                     <  6:355007/149‑1 (MQ=255)                                                                 gCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACgg                                                                                                           <  1:384287/123‑1 (MQ=255)                                                                 gCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACgg                                                                                                           >  2:384287/1‑123 (MQ=255)                                                                 gCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGGTCACGGTATCGCCCGGTGCGTACGCCGCCAcg                                                                                 >  7:408942/1‑149 (MQ=255)                                                                         tAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAgcg                                                                         >  2:40893/1‑149 (MQ=255)                                                                           gCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                       <  2:316701/149‑1 (MQ=255)                                                                            ccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAg                                                                      >  5:397932/1‑149 (MQ=255)                                                                            ccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAg                                                                      >  5:403317/1‑149 (MQ=255)                                                                            ccGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAg                                                                      >  7:294106/1‑149 (MQ=255)                                                                              ggAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGtt                                                                    <  6:73418/149‑1 (MQ=255)                                                                                aaCGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCg                                                                  <  3:57028/149‑1 (MQ=255)                                                                                  cGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGt                                                                <  3:127838/149‑1 (MQ=255)                                                                                     aGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAg                                                             >  3:361615/1‑149 (MQ=255)                                                                                           gggaaTGGCAGGCAGTCTTGCCGCTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTGTCGGCCAGGTACACGGTATCGCCCGGTGCGTACCCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAc                                                       >  7:274755/4‑149 (MQ=255)                                                                                              aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAcgtc                                                    <  6:50416/149‑1 (MQ=255)                                                                                              aaTGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAcgtc                                                    <  1:169209/149‑1 (MQ=255)                                                                                                 ggcaggcaGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGtt                                                 >  4:356632/1‑149 (MQ=255)                                                                                                   caggcagTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGc                                               <  8:294106/149‑1 (MQ=255)                                                                                                    aggcagTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGcc                                              >  5:162844/1‑149 (MQ=255)                                                                                                    aggcagTCTGGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGcc                                              <  6:263083/149‑1 (MQ=255)                                                                                                            ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCa                                      >  7:399739/1‑149 (MQ=255)                                                                                                              gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCa                                      <  2:373471/147‑1 (MQ=255)                                                                                                              gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATg                                    >  6:381447/1‑149 (MQ=255)                                                                                                                 aCTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGc                                                                                              <  1:286807/88‑1 (MQ=255)                                                                                                                 aCTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGc                                                                                              >  2:286807/1‑88 (MQ=255)                                                                                                                 aCTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCg                                  >  6:348064/1‑148 (MQ=255)                                                                                                                  cTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTa                                <  1:40893/149‑1 (MQ=255)                                                                                                                     ggCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCt                             <  5:348064/149‑1 (MQ=255)                                                                                                                       cccccGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGc                           <  5:206948/149‑1 (MQ=255)                                                                                                                           cGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAgg                       >  3:91471/1‑149 (MQ=255)                                                                                                                            ggCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGc                      <  5:109537/149‑1 (MQ=255)                                                                                                                             gCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCCCGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCa                     >  4:109938/1‑149 (MQ=255)                                                                                                                                gCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                       <  7:291631/96‑1 (MQ=255)                                                                                                                                gCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                       >  8:291631/1‑96 (MQ=255)                                                                                                                                 cAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTAGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGt                 >  6:414303/1‑149 (MQ=255)                                                                                                                                      cACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcg            >  1:178304/1‑149 (MQ=255)                                                                                                                                      cACCTGGCCCGCGGTATCGCAGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcg            >  8:93331/1‑149 (MQ=255)                                                                                                                                       aCCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggcg            <  5:381447/148‑1 (MQ=255)                                                                                                                                         cTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa         >  3:196540/1‑149 (MQ=255)                                                                                                                                          tGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa         >  4:365014/1‑148 (MQ=255)                                                                                                                                            gCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa         <  3:365014/146‑1 (MQ=255)                                                                                                                                            gCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGAc        <  2:145683/147‑1 (MQ=255)                                                                                                                                            gCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACa       <  6:162844/148‑1 (MQ=255)                                                                                                                                            gCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAg      >  3:245745/1‑149 (MQ=255)                                                                                                                                             cccGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgt     >  1:212009/1‑149 (MQ=255)                                                                                                                                               cgcgGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgtgt   >  3:201784/1‑149 (MQ=255)                                                                                                                                               cgcgGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgtgt   >  1:407253/1‑149 (MQ=255)                                                                                                                                                gcgGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTGCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGtt  <  5:12608/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATGGAGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTT  >  NC_002947/1457000‑1457291

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