Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,711,029	T→G	7.8%	intergenic (+289/+25)	PP_2373 → / ← PP_2374	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,711,029	0	T	G	7.8%	205.2 / 10.4	77	intergenic (+289/+25)	PP_2373/PP_2374	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/33);  new base G (2/4);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

ACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTCGCG  >  NC_002947/2710881‑2711170                                                                                                                                                      |                                                                                                                                              aCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTc                                                                                                                                              >  2:32394/1‑150 (MQ=255)  cccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCa                                                                                                                                             >  1:751193/1‑150 (MQ=255)   ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAc                                                                                                                                            >  1:995702/1‑150 (MQ=255)     gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAccc                                                                                                                                          >  2:997609/1‑150 (MQ=255)     gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAccc                                                                                                                                          >  1:279882/1‑150 (MQ=255)     gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAccc                                                                                                                                          >  2:1133838/1‑150 (MQ=255)         aaGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTGAACCcgcgc                                                                                                                                      <  2:545187/150‑2 (MQ=255)          aGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGaa                                                                                                                                     <  1:668812/150‑1 (MQ=255)          aGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTGAACCcgcgcc                                                                                                                                     <  2:342435/150‑3 (MQ=255)                gggccTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCCCGGGTTCACCCGCGCAGAGGcc                                                                                                                               >  2:174653/1‑150 (MQ=255)                   ccTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGc                                                                                                                            >  1:915378/1‑150 (MQ=255)                     tgctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGcac                                                                                                                          <  1:1006548/150‑1 (MQ=255)                      gctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCGGCTCCCG‑TATTCGCGGGTAAACCCGCTTCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGcaca                                                                                                                         <  2:785561/150‑1 (MQ=255)                      gctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAc                                                                                                                                            <  2:487225/131‑1 (MQ=255)                      gctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCAc                                                                                                                                            >  1:487225/1‑131 (MQ=255)                       ctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCcg                                                                                                                                         >  2:1019124/1‑133 (MQ=255)                       ctgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCcg                                                                                                                                         <  1:1019124/133‑1 (MQ=255)                        tgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAg                                                                                                                        >  2:810925/1‑149 (MQ=255)                        tgGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAg                                                                                                                        <  1:810925/149‑1 (MQ=255)                            cAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAg                                                                                                                   >  2:1348235/1‑150 (MQ=255)                             agtttCCACCTGCCGTTCACCGCTGCCGCAACTCCCGCTTTTCG‑GTTGAAGTGCGCTCCTACGCCTGAACCTGGGGGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGc                                                                                                                  <  1:672353/148‑1 (MQ=21)                              aTTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCg                                                                                                                 <  1:702983/150‑1 (MQ=255)                               tttCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGTGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCgg                                                                                                                <  1:1217335/150‑1 (MQ=255)                                  ccACCGGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTc                                                                                                             <  2:500882/150‑1 (MQ=255)                                   cACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCa                                                                                                            <  2:915378/150‑1 (MQ=255)                                        gCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGc                                                                                                       <  1:1348235/150‑1 (MQ=255)                                        gCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGc                                                                                                       >  1:595139/1‑150 (MQ=255)                                           gTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCAAAGGATAGCAGGTCAGGAGCatg                                                                                                    >  1:21618/1‑148 (MQ=255)                                               aCCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCa                                                                                                >  1:465943/1‑150 (MQ=255)                                                       cgccgcTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCACCTGCa                                                                                          >  1:1158464/1‑148 (MQ=255)                                                       cgccgcTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCa                                                                                          <  2:1158464/148‑1 (MQ=255)                                                             tACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGCCACCGGACGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCg                                                                                  >  2:1267103/1‑150 (MQ=255)                                                             tACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCg                                                                                  >  1:496511/1‑150 (MQ=255)                                                                   tATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCaa                                                                             >  1:1122226/1‑150 (MQ=255)                                                                               aaCCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGaa                                                                 <  2:445158/150‑1 (MQ=255)                                                                                 cccGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAaca                                                               >  2:104075/1‑150 (MQ=255)                                                                                  ccGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAacac                                                              >  1:902327/1‑150 (MQ=255)                                                                                         tACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTg                                                       >  1:305906/1‑150 (MQ=255)                                                                                            gCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGaa                                                                 >  2:565813/1‑137 (MQ=255)                                                                                            gCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGaa                                                                 <  1:565813/137‑1 (MQ=255)                                                                                            gCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAAc                                                    >  2:400240/1‑150 (MQ=255)                                                                                            gCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAAc                                                    >  1:756459/1‑150 (MQ=255)                                                                                             ccTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACa                                                   <  1:997609/150‑1 (MQ=255)                                                                                              cTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAAt                                                                    >  2:711363/1‑132 (MQ=255)                                                                                              cTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAAt                                                                    <  1:711363/132‑1 (MQ=255)                                                                                               tGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAgc                                                 <  2:1122226/150‑1 (MQ=255)                                                                                                        gTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAgcgc                                               <  1:947325/143‑1 (MQ=255)                                                                                                        gTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAgcgc                                               >  2:947325/1‑143 (MQ=255)                                                                                                            acacAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCaa                                    <  2:279882/150‑1 (MQ=255)                                                                                                             cacaGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAAt                                   <  2:595139/150‑1 (MQ=255)                                                                                                                aGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                          <  2:1350457/124‑1 (MQ=255)                                                                                                                aGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAgg                                                          >  1:1350457/1‑124 (MQ=255)                                                                                                                aGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGaa                                <  1:400240/150‑1 (MQ=255)                                                                                                                 ggTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAg                               >  1:1280078/1‑150 (MQ=255)                                                                                                                  gTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGACGTCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGt                              >  2:555430/1‑150 (MQ=255)                                                                                                                      ggggTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCAc                                                                                         >  2:341682/1‑87 (MQ=255)                                                                                                                      ggggTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCAc                                                                                         <  1:341682/87‑1 (MQ=255)                                                                                                                       gggTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGtggctg                         <  2:1325381/150‑1 (MQ=255)                                                                                                                          ttGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGtggctggc                       <  1:722789/149‑1 (MQ=255)                                                                                                                          ttGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGtggctggc                       >  2:722789/1‑149 (MQ=255)                                                                                                                             gACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTg                   >  2:255200/1‑150 (MQ=255)                                                                                                                              acacAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCAc                                                                                         >  2:864632/1‑79 (MQ=255)                                                                                                                              acacAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCAc                                                                                         <  1:864632/79‑1 (MQ=255)                                                                                                                              acacAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGc                  >  1:138561/1‑150 (MQ=255)                                                                                                                                  aGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGcc                      <  2:863906/142‑1 (MQ=255)                                                                                                                                  aGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGcc                      >  1:863906/1‑142 (MQ=255)                                                                                                                                   ggCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTgg             <  2:756459/150‑1 (MQ=255)                                                                                                                                   ggCCGCGGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTgg             <  1:1267103/150‑1 (MQ=255)                                                                                                                                    gtgtggctAGGGGCGGGTTCCCCCGCGAAAAGGACGGCACAGGTTAGCGGGTCAAGCGCAGGGCCAAGTACAAGAGCCCGATGAATAGCGCAATTAACAGCAGGGTGAACAGCGTACACAGCGCAATGAAGTGGCTGGCCTTGCAGTggg            >  2:599502/8‑150 (MQ=255)                                                                                                                                       gCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGt                              <  2:1071406/129‑1 (MQ=255)                                                                                                                                       gCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGt                              >  1:1071406/1‑129 (MQ=255)                                                                                                                                        cTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGaa        <  1:1073068/150‑1 (MQ=255)                                                                                                                                          ctcttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGt                                                        >  2:505371/6‑100 (MQ=38)                                                                                                                                          ctcttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGt                                                        <  1:505371/95‑1 (MQ=38)                                                                                                                                             ttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTgg             <  2:1090841/138‑1 (MQ=255)                                                                                                                                             ttcgcgGGTGAACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTgg             >  1:1090841/3‑140 (MQ=255)                                                                                                                                              gcgcGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTcgcg  <  2:902327/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              ACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACGCCTGAACGAGGGTGGACACAGGTCGGGGGTTGGACACAGGCCGCTGTAGGCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGTTAGCGGGTCAGGTGCAGGGCCAACTGCACCAGCCCGATCAATACCGCAATGAACACCAGGGTGAACAGCGCCCCCAGCGCAATGAAGTGGCTGGCCTTGCCGTGGGTGAAGTCGCG  >  NC_002947/2710881‑2711170

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