Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,383,801	T→C	10.3%	S47G (AGC→GGC)	PP_2091 ←	serine/threonine‑protein kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,383,801	0	T	C	10.3%	180.3 / 7.2	66	S47G (AGC→GGC)	PP_2091	serine/threonine‑protein kinase
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (3/4);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCCTGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCG  >  NC_002947/2383655‑2383944                                                                                                                                                   |                                                                                                                                                gCCCTGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGAc                                                                                                                                              <  2:689967/150‑1 (MQ=255) gCCCTGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAAGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGTCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGAc                                                                                                                                              <  1:520928/150‑1 (MQ=255)  cccTGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGAGCGCCGCCAGGCCGCCGTCGGTACATGTGCTGACg                                                                                                                                             >  2:912922/1‑150 (MQ=255)     tGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCg                                                                                                                                          >  1:335314/1‑150 (MQ=255)          gcgcAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGc                                                                                                                                     >  1:1432560/1‑150 (MQ=255)           cgcAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGcc                                                                                                                                    >  2:1187365/1‑150 (MQ=255)            gcAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCa                                                                                                                                   >  2:197433/1‑150 (MQ=255)             cAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAg                                                                                                                                  <  2:1095837/150‑1 (MQ=255)             cAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGGCGGTACATTGGCCGACGGCGTCAGCCaa                                                                                                                                  <  2:598837/150‑2 (MQ=255)              aaCCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCaat                                                                                                                                 <  1:912922/150‑3 (MQ=255)              aaCCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAgg                                                                                                                                 >  1:527970/1‑150 (MQ=255)              aaCCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAgg                                                                                                                                 <  1:1574703/150‑1 (MQ=255)               aCCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAggg                                                                                                                                <  1:1187365/150‑1 (MQ=255)                 cAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCg                                                                                                                              >  2:573419/1‑150 (MQ=255)                  aGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGa                                                                                                                             >  2:96601/1‑150 (MQ=255)                      gTTGTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCAAGGTGGCCGCCGTCGCGTCGTTCTCCCGGTCCGGTGGGTGGAGGCTGGCGTTCGCCCTCAGGGCGTCGTCGGTACATTGGTTGACGCCGGCGGCCAGGGCGAACAg                                                                                                                         <  2:596793/150‑1 (MQ=255)                       ttttGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAgg                                                                                                                        >  2:843000/1‑150 (MQ=255)                                   agcagGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCCGACGGCGTCAGc                                                                                                                                     <  1:51059/125‑1 (MQ=255)                                   agcagGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCCGACGGCGTCAGc                                                                                                                                     >  2:51059/1‑125 (MQ=255)                                     cagGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTggc                                                                                                          <  2:729193/150‑1 (MQ=255)                                     cagGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTggc                                                                                                          >  2:511989/1‑150 (MQ=255)                                        gCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGCCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTggcggc                                                                                                       >  1:1572577/1‑150 (MQ=255)                                              gAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTc                                                                                                 >  2:1288371/1‑150 (MQ=255)                                                  gCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGt                                                                                             <  1:1451636/150‑1 (MQ=255)                                                   ccTGGGCCACGCGCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTg                                                                                            >  2:1432678/1‑150 (MQ=255)                                                    cTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGc                                                                                           >  1:201161/1‑150 (MQ=255)                                                     ggggCCCCGCTTCAGGTGGCCGGCGTTGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCGGCCAGGGCGCCGTCGGTACGTTGGCTGACGCCGTCGGCGAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGcc                                                                                          <  2:30543/149‑1 (MQ=255)                                                          cACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGt                                                                                     >  2:206817/1‑150 (MQ=255)                                                               tCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCa                                                                                >  1:247410/1‑150 (MQ=255)                                                                       gCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAgg                                                                        <  1:1309114/150‑1 (MQ=255)                                                                       gCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAgg                                                                        <  1:1125110/150‑1 (MQ=255)                                                                         cGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCg                                                                                                                                       <  2:707839/85‑1 (MQ=255)                                                                         cGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCg                                                                                                                                       >  1:707839/1‑85 (MQ=255)                                                                         cGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCCGACGGCGTCAGCCaa                                                                                                                                  >  1:1315723/1‑89 (MQ=255)                                                                         cGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCCGACGGCGTCAGCCaa                                                                                                                                  <  2:1315723/90‑2 (MQ=255)                                                                             gtggcgtAGTAATCCAGGGCCCGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTNTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTcc                                                                  <  1:286943/150‑1 (MQ=255)                                                                              tggcgtAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCtt                                                                                                                 <  1:833292/102‑1 (MQ=255)                                                                              tggcgtAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCtt                                                                                                                 >  2:833292/1‑102 (MQ=255)                                                                                    agtaATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGc                                                                                                                    <  1:1540792/93‑1 (MQ=255)                                                                                    agtaATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGc                                                                                                                    >  2:1540792/1‑93 (MQ=255)                                                                                    agtaATCCAGGGCCAGCGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGtttt                                                           <  2:906716/150‑1 (MQ=255)                                                                                      taATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTcc                                                         <  2:335314/150‑1 (MQ=255)                                                                                      taATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTcc                                                         <  1:908839/150‑1 (MQ=255)                                                                                        aTCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAgg                                                                                                                        <  1:854775/85‑1 (MQ=255)                                                                                        aTCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAgg                                                                                                                        >  2:854775/1‑85 (MQ=255)                                                                                         tCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCcgcg                                                      >  2:1404794/1‑150 (MQ=255)                                                                                            aGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGaa                                                                                                                            >  1:229452/1‑77 (MQ=255)                                                                                            aGGGCCAGTGCCTGCAGGCTGGCGAGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGaa                                                                                                                            <  2:229452/77‑1 (MQ=255)                                                                                             gggCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCcgcgcgcg                                                  <  1:96601/150‑1 (MQ=255)                                                                                                ccAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTg                                                                                               <  1:1272774/102‑1 (MQ=255)                                                                                                ccAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTg                                                                                               >  2:1272774/1‑102 (MQ=255)                                                                                                    tGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGt                                           <  2:1432560/150‑1 (MQ=255)                                                                                                    tGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGt                                           <  1:1404794/150‑1 (MQ=255)                                                                                                       cTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGgc                                        <  2:464313/150‑1 (MQ=255)                                                                                                       cTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGgc                                        >  1:858114/1‑150 (MQ=255)                                                                                                            ggCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCGAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGACGCCGGGTgcgcgcga                                   <  2:565233/150‑8 (MQ=255)                                                                                                                ggCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTg                               <  1:1288371/150‑1 (MQ=255)                                                                                                                 gcgcgcgcCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTgg                              <  1:511989/150‑1 (MQ=255)                                                                                                                       gccgccAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAg                        <  1:843000/150‑1 (MQ=255)                                                                                                                         cgccAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCt                      <  2:247410/150‑1 (MQ=255)                                                                                                                          gccAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCCAAGCTTCGTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTc                     >  1:100355/1‑150 (MQ=255)                                                                                                                                   ccgTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAgg            <  2:1572577/150‑1 (MQ=255)                                                                                                                                      tCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTc         >  2:1371117/1‑150 (MQ=255)                                                                                                                                             aTTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCg  >  1:733649/1‑150 (MQ=255)                                                                                                                                              ctGGCCGACGGCGTCAGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGtt                                                             <  2:1495809/89‑1 (MQ=255)                                                                                                                                              ctGGCCGACGGCGTCAGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGtt                                                             >  1:1495809/2‑90 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCCCTGGGCGCGCAACCAGCGGTTTTGCGCCAGCAGCAGGCGGTCGAGGGCCTGGGCCACGCTCCAGGTGGCCGGCGTGGCGTAGTAATCCAGGGCCAGTGCCTGCAGGCTGGCGCGCGCCGCCAGGCCGCCGTCGGTACATTGGCTGACGCCGTCGGCCAGGGCGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCG  >  NC_002947/2383655‑2383944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑