Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,673,189	G→T	9.9%	S168I (AGC→ATC)	PP_2340 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,673,189	0	G	T	9.9%	273.8 / 15.4	90	S168I (AGC→ATC)	PP_2340	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/38);  new base T (4/5);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

CTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTGGTGGG  >  NC_002947/2673046‑2673335                                                                                                                                                 |                                                                                                                                                   cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                               <  7:342957/149‑1 (MQ=255) cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCAtcgcgc                                                                                                                                               <  7:146867/149‑1 (MQ=255) cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCAtcgcg                                                                                                                                                >  4:327986/1‑148 (MQ=255)    tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                              <  3:386822/147‑1 (MQ=255)    tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                             <  6:50457/148‑1 (MQ=255)    tGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                             <  1:187193/148‑1 (MQ=255)      gCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaa                                                                                                                                          <  3:327986/149‑3 (MQ=255)              ttGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑AACGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGc                                                                                                                                  <  1:335939/149‑1 (MQ=255)                 ctgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCTGCGTTATATCGCCCCCAGCTTCCTGCGCTATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCg                                                                                                                                <  7:259509/148‑1 (MQ=255)                  tgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCggg                                                                                                                              >  7:387423/1‑149 (MQ=255)                   gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                                   <  6:332014/127‑1 (MQ=255)                   gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgc                                                                                                                                                   >  5:332014/1‑127 (MQ=255)                   gcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                             >  1:184536/1‑149 (MQ=255)                    ccGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCg                                                                                                                            <  6:375034/149‑1 (MQ=255)                            ccGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTc                                                                                                                    >  6:42842/1‑149 (MQ=255)                                  cTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                                      <  6:996/125‑1 (MQ=255)                                  cTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGt                                                                                                                                      >  5:996/1‑125 (MQ=255)                                     ggAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAg                                                                                                                                             >  8:331667/1‑115 (MQ=255)                                      gAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                          >  7:241678/1‑149 (MQ=255)                                      gAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                          >  3:131471/1‑149 (MQ=255)                                       aaCTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCATGAGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCCGTCGCCGGTGTAtt                                                                                                         >  8:167035/1‑149 (MQ=255)                                            ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                                    >  2:19023/1‑149 (MQ=255)                                            ccccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCaa                                                                                                    <  8:387423/149‑1 (MQ=255)                                                aaTTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGc                                                                                                >  6:409110/1‑149 (MQ=255)                                                 aactcgcgAT‑CGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCt                                                                                               >  1:154276/5‑149 (MQ=255)                                                     cgcgATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTg                                                                                           >  5:171040/1‑149 (MQ=255)                                                      gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                             <  5:2771/114‑1 (MQ=255)                                                      gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGc                                                                                                                             >  6:2771/1‑114 (MQ=255)                                                        gatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCCGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAg                                                                                        <  3:242327/149‑1 (MQ=255)                                                          tgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                                         <  3:62726/114‑1 (MQ=255)                                                          tgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAAt                                                                                                                         >  4:62726/1‑114 (MQ=255)                                                             gatcg‑atcGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTc                                                                                   >  2:279175/1‑149 (MQ=255)                                                              atcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCg                                                                                  >  6:340809/1‑149 (MQ=255)                                                               tcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGCCAGCGCGCAGCGCCAGTGAGCCGGGCGAATCAGTCGCCCGTGTAt                                                                                                          >  6:210777/1‑124 (MQ=255)                                                               tcg‑atcgCGAGTGCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                          <  5:210777/124‑1 (MQ=255)                                                                    tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCTGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                             <  8:241678/149‑1 (MQ=255)                                                                    tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                             >  6:278358/1‑149 (MQ=255)                                                                    tcgCGAGTTCATCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                             <  4:29490/149‑1 (MQ=255)                                                                        gAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCg                                                                          >  6:112593/1‑148 (MQ=255)                                                                            ccttcCTGCTTTCAGTGTCGCACTACCAGCGCTATATCGCCCCCAGGTTCCTGCGCGATGCCATGGGCAGCGCGGAGGGCCAGTGAGGCGGGCGAATCATTCGCCGGTGTATTCGCAACCGTTGGTGCAGGTTTCGTGGATTCGAACCt                                                                     <  1:219745/148‑1 (MQ=255)                                                                             cttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCt                                                                     <  8:325667/148‑1 (MQ=255)                                                                              ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                                   >  3:292433/1‑149 (MQ=255)                                                                              ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                                   >  1:350120/1‑149 (MQ=255)                                                                              ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                                   >  6:367394/1‑149 (MQ=255)                                                                              ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                                         >  6:422798/1‑76 (MQ=21)                                                                              ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCAggaag                                                                                                                                         <  5:422798/79‑4 (MQ=21)                                                                                cttGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCCGCTTCCTGCGCTATGCCTTGGGCAGGGCGCAGGGCCAGTGAGCCTGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                                    <  2:184536/144‑1 (MQ=255)                                                                                ccTGGTTTCAGAGTCGCACTACCAGCGGTATATCGCCCCGAGCTTCCTGCGCTTAGCCATGGGCAGCGCGAAGGGCCAGTGAGCGGTGCGAATCAGTCGCCGGTGTATTCGCAACCGGTGGTGCAGGTTTCGt                                                                                 <  7:232913/133‑1 (MQ=255)                                                                                ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCtt                                                                    <  4:131471/146‑1 (MQ=255)                                                                                ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                                   <  1:19023/147‑1 (MQ=255)                                                                                ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                                   <  4:292433/147‑1 (MQ=255)                                                                                ccTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTACTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                                 >  8:232913/1‑133 (MQ=255)                                                                                 cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAc                                                                <  4:380811/149‑1 (MQ=255)                                                                                  tGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                                 >  3:349276/1‑147 (MQ=255)                                                                                    cTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGc                                                             >  4:127750/1‑149 (MQ=255)                                                                                          gagGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                                       <  6:171040/149‑1 (MQ=255)                                                                                          gagGCGCACTACCAGCGCTATATCG‑CCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                                      >  7:387922/1‑149 (MQ=255)                                                                                           agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                                       <  2:350120/148‑1 (MQ=255)                                                                                           agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                                      >  6:2643/1‑149 (MQ=255)                                                                                            ggCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCAt                                                     >  7:195906/1‑149 (MQ=255)                                                                                             gcgcACTACCAGCGATATATCCCCCCCCGCGTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATc                                                    <  2:349671/149‑1 (MQ=255)                                                                                              cgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                                   <  8:195906/149‑1 (MQ=255)                                                                                               gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAg                                                  >  1:222671/1‑149 (MQ=255)                                                                                                     ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGc                                              <  1:279175/147‑1 (MQ=255)                                                                                                     ccAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                            <  3:127750/149‑1 (MQ=255)                                                                                                       aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                                          >  4:17232/1‑149 (MQ=255)                                                                                                        gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                               >  5:240185/1‑127 (MQ=255)                                                                                                        gcgcTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGAc                                         <  5:340809/149‑1 (MQ=255)                                                                                                            tataTCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTgg                                     >  2:175996/1‑149 (MQ=255)                                                                                                              tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                            >  4:32912/1‑140 (MQ=255)                                                                                                              tatCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCtt                                            <  3:32912/140‑1 (MQ=255)                                                                                                                tCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                                  >  3:212559/1‑148 (MQ=255)                                                                                                                  gCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTcc                                  >  2:141569/1‑146 (MQ=255)                                                                                                                   cccccAGCTTCCTGCGCGTTGCCATGGGAAGCGCGCAGGGCCAGTGAGCCCGGCGAATCAGTCCCCGGTGTATTAGCCACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTGCGGCATCAGCGGCTTGACCTGGTcc                                  >  5:84316/1‑145 (MQ=255)                                                                                                                   cccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                              <  5:278358/149‑1 (MQ=255)                                                                                                                   cccccAGCTTCCCGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGa                              <  5:367394/149‑1 (MQ=255)                                                                                                                    ccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGAt                             <  1:175996/149‑1 (MQ=255)                                                                                                                      ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                              >  1:142972/1‑82 (MQ=255)                                                                                                                      ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                              <  2:142972/82‑1 (MQ=255)                                                                                                                        aGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCAc                         <  5:409110/149‑1 (MQ=255)                                                                                                                            cccTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTgg                     >  2:61501/2‑149 (MQ=255)                                                                                                                              cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                        <  1:232139/64‑1 (MQ=21)                                                                                                                              cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                        >  2:232139/1‑64 (MQ=21)                                                                                                                              cTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGc                    <  4:212559/148‑1 (MQ=255)                                                                                                                               tGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGa                  >  8:46422/1‑149 (MQ=255)                                                                                                                                 cgcgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt                 <  5:42842/148‑1 (MQ=255)                                                                                                                                   cgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCAc              <  7:46422/149‑1 (MQ=255)                                                                                                                                    gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc                <  6:251632/146‑1 (MQ=255)                                                                                                                                    gATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATc                >  5:251632/1‑146 (MQ=255)                                                                                                                                        ccATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGc         >  5:285927/1‑149 (MQ=255)                                                                                                                                        ccATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGc         >  3:367669/1‑149 (MQ=255)                                                                                                                                         cATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCt        <  5:270363/149‑1 (MQ=255)                                                                                                                                         cATGGGCAGCGAGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGAAAGCTCCGGCATCATCGGCTTTACCTGGTCCCAGCTACACTTGGCGATCAATTctct        >  6:326193/1‑149 (MQ=255)                                                                                                                                               cAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGTGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTGGTggg  >  7:178168/1‑149 (MQ=255)                                                                                                                                               cAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTGGTggg  >  5:391130/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTGGTGGG  >  NC_002947/2673046‑2673335

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