Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,936,329	C→G	7.1%	A24G (GCC→GGC)	PP_3469 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,936,329	0	C	G	7.1%	219.6 / 8.0	71	A24G (GCC→GGC)	PP_3469	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base G (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAG  >  NC_002947/3936185‑3936471                                                                                                                                                 |                                                                                                                                               tCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgcg                                                                                                                                           >  2:609966/1‑150 (MQ=255) tCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgcg                                                                                                                                           >  2:697520/1‑150 (MQ=255) tCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgcg                                                                                                                                           >  1:475100/1‑150 (MQ=255)   aaCTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCt                                                                                                                                         >  2:907197/1‑150 (MQ=255)     cTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTAc                                                                                                                                       <  1:1013233/150‑1 (MQ=255)      tGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACg                                                                                                                                      >  1:902852/1‑150 (MQ=255)      tGTCAAATCTATCAGTCGCCCAAGAGATATTTCTAGAACATCATCAACAACTAAATGCGGCAGGTTGCAATGACCAGCCCGCAGCATCTACGACACTTTTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACg                                                                                                                                      >  1:988695/1‑150 (MQ=255)       gTCAAAGCTATCAGTCGTCCAAGGGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGGCAGCGCttgtc                                                                                                                                     <  1:955223/150‑5 (MQ=255)        tCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGc                                                                                                                                                >  1:237778/1‑138 (MQ=255)        tCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGc                                                                                                                                                <  2:237778/138‑1 (MQ=255)         cAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGGCAGCGCttgtcag                                                                                                                                   <  2:845973/150‑7 (MQ=255)          aaaTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttct                                                                                                                                  >  2:840000/1‑150 (MQ=255)          aaaTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttct                                                                                                                                  >  1:310570/1‑150 (MQ=255)          aaaTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttct                                                                                                                                  <  1:327754/150‑1 (MQ=255)           aaTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttctt                                                                                                                                 <  2:108869/150‑1 (MQ=255)            aTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttcttc                                                                                                                                <  1:907197/150‑1 (MQ=255)            aTCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGttcttc                                                                                                                                >  1:381223/1‑150 (MQ=255)                agcggTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTg                                                                                                                            >  2:536744/5‑150 (MQ=255)                   aGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGAc                                                                                                                          >  1:717591/1‑149 (MQ=255)                   aGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGAc                                                                                                                          <  2:717591/149‑1 (MQ=255)                   aGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACa                                                                                                                         >  1:321365/1‑150 (MQ=255)                    gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGGCAGCGCttgtca                                                                                                                                    >  2:1128795/1‑133 (MQ=255)                    gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGGCAGCGCttgtca                                                                                                                                    <  1:1128795/138‑6 (MQ=255)                    gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAg                                                                                                                        <  2:988695/150‑1 (MQ=255)                    gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAg                                                                                                                        >  1:883229/1‑150 (MQ=255)                         ccAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCg                                                                                                                   >  1:834208/1‑150 (MQ=255)                           aaGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                 <  2:381223/150‑1 (MQ=255)                                  atTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTgg                                                                                                          <  1:677642/150‑1 (MQ=255)                                         gACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAatt                                                                                                   <  2:475100/150‑1 (MQ=255)                                            cATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAattatt                                                                                                <  2:902852/150‑1 (MQ=255)                                                 tCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAattatt                                                                                                >  2:436381/1‑145 (MQ=255)                                                 tCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAattatt                                                                                                <  1:436381/145‑1 (MQ=255)                                                 tCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATtgtgtg                                                                                           >  1:1862/1‑150 (MQ=255)                                                      ccAAAACCTCGGTAGATCTCAATGACCGGCCCGCAGCATCAACGACACTTCTACATGCATAAGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGcc                                                                                      >  2:171249/1‑150 (MQ=255)                                                         aaaCCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAAt                                                                                   <  1:697520/150‑1 (MQ=255)                                                          aaCCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATc                                                                                  <  1:840000/150‑1 (MQ=255)                                                              tCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATCCGCACTTACACTCGGTGCTTACAAGCGCTGCCAGCCCTACGTTCTTCCGTGACAGCTCCGTACTACTGGCAGAATTATTGTGTGAAGCCAAACAGAc                                                                              >  2:448416/1‑150 (MQ=255)                                                               cGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACa                                                                             <  1:536744/150‑1 (MQ=255)                                                                 gTAGATCTCAATGACCAGACCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGc                                                                           >  1:270588/1‑150 (MQ=255)                                                                 gTAGATCTCAATGAACAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCGCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGc                                                                           <  1:196172/150‑1 (MQ=255)                                                                    gATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGAc                                                                        <  1:609966/150‑1 (MQ=255)                                                                        tcAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGa                                                                               <  2:464464/139‑1 (MQ=255)                                                                        tcAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGa                                                                               >  1:464464/1‑139 (MQ=255)                                                                          aaTGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGAAACATTTCTGCGACTGGCCGATTTATTGTGTGAAGCCAATCAGACAGCGACCACGCt                                                                  >  2:33764/1‑150 (MQ=255)                                                                               ccAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                 >  1:353756/1‑98 (MQ=255)                                                                               ccAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                 <  2:353756/98‑1 (MQ=255)                                                                                   cccGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGctc                                                         <  2:1157460/150‑1 (MQ=255)                                                                                          taTCTACGACACTTCTACATGTATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  <  1:448416/149‑1 (MQ=255)                                                                                          caTCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  <  2:310570/150‑1 (MQ=255)                                                                                               aCGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAg                                             >  1:356317/1‑150 (MQ=255)                                                                                                    acTTCTACATGCATACGCACTTACGCTCGGTGCTGACAAGCGCTGACAGCGCTACGTTCTTCCGTGACAGCTCAGTAGTACTGGCCGAATCAGTGTGTGGAGACCATCAGACAGAGACCACGCTTGTCGGCTCTGACATGGCCAGAGCAc                                        >  2:537853/1‑150 (MQ=255)                                                                                                     cTTCGACACGCCTACGCACTTACACTTGCTGCTGACAAGCGCTGCCAGCGCCACGTTTTTCCGTGACCGCCCCGTACTACTGCCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCAcc                                       <  1:19364/150‑1 (MQ=255)                                                                                                            aTGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTg                                                                >  2:729500/1‑118 (MQ=255)                                                                                                            aTGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTg                                                                <  1:729500/118‑1 (MQ=255)                                                                                                            aTGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAAcc                                >  1:462486/1‑150 (MQ=255)                                                                                                              gCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCg                              <  2:356317/150‑1 (MQ=255)                                                                                                                aTACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                 >  2:1075145/1‑65 (MQ=255)                                                                                                                aTACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                 <  1:1075145/65‑1 (MQ=255)                                                                                                                aTACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCGGATCACCCATAACCCGcc                            <  2:321365/150‑2 (MQ=255)                                                                                                                 tACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAg                           >  1:1038330/1‑150 (MQ=255)                                                                                                                  aCGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  <  2:374569/126‑1 (MQ=255)                                                                                                                  aCGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  >  1:374569/1‑126 (MQ=255)                                                                                                                  aCGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  >  2:1088610/1‑126 (MQ=255)                                                                                                                  aCGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATg                                                  <  1:1088610/126‑1 (MQ=255)                                                                                                                    gCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTAc                        >  1:806822/1‑150 (MQ=255)                                                                                                                      aCTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACgg                      >  1:593564/1‑150 (MQ=255)                                                                                                                          acacTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGa                                                      >  2:581996/1‑114 (MQ=255)                                                                                                                          acacTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGa                                                      <  1:581996/114‑1 (MQ=255)                                                                                                                           ctcgcGGTGCTGACAATCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATGATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAac                 <  1:537853/146‑1 (MQ=255)                                                                                                                           cacTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAac                 <  2:1038330/150‑1 (MQ=255)                                                                                                                                       cgtaGCGCTGGCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGGGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCg     >  2:1139648/4‑150 (MQ=255)                                                                                                                                       aCAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCg     >  1:63230/1‑150 (MQ=255)                                                                                                                                          aGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACAATTTTTACGGAAAACACCTCCTTACCGCAg  <  1:33764/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               TCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAG  >  NC_002947/3936185‑3936471

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