Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,268,580	G→C	6.1%	Y292* (TAC→TAG)	PP_2867 ←	pyridine nucleotide‑disulphide oxidoreductase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,268,580	0	G	C	6.1%	205.3 / 3.5	65	Y292* (TAC→TAG)	PP_2867	pyridine nucleotide‑disulphide oxidoreductase family protein
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base C (2/2);  total (28/37)
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.68e-01

TCCAGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCGTC  >  NC_002947/3268431‑3268715                                                                                                                                                       |                                                                                                                                        tCCAGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCg                                                                                                                                         <  2:612902/150‑1 (MQ=255)    aGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAg                                                                                                                                      <  2:85268/150‑1 (MQ=255)      gCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGc                                                                                                                                    <  1:937776/150‑1 (MQ=255)      gCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGc                                                                                                                                    >  1:1311206/1‑150 (MQ=255)                tACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCa                                                                                                                          >  2:1025753/1‑150 (MQ=255)                  cTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGc                                                                                                                        <  2:1044372/150‑1 (MQ=255)                   ttGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCt                                                                                                                       <  1:1055350/150‑1 (MQ=255)                   ttGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCt                                                                                                                       >  1:1235039/1‑150 (MQ=255)                   ttGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCt                                                                                                                       >  1:932239/1‑150 (MQ=255)                        ctggctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATa                                                                                                                  >  1:1028530/1‑150 (MQ=255)                        ctggctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATa                                                                                                                  >  2:887881/1‑150 (MQ=255)                        ctggctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATa                                                                                                                  >  1:747736/1‑150 (MQ=255)                          ggctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATaca                                                                                                                <  2:839233/150‑1 (MQ=255)                           gctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATacac                                                                                                               >  1:1306438/1‑150 (MQ=255)                           gctgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATacac                                                                                                               <  2:51636/150‑1 (MQ=255)                             tgTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACgg                                                                                                             >  2:1024101/1‑150 (MQ=255)                                 aGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGcatc                                                                                                         >  1:405125/1‑150 (MQ=255)                                    cGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCAtt                                                                                                      <  1:1207683/150‑1 (MQ=255)                                       tGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGcc                                                                                                   <  2:1235039/150‑1 (MQ=255)                                          tGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCtt                                                                                                <  2:1306438/150‑1 (MQ=255)                                            cTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGa                                                                                              <  2:1028530/150‑1 (MQ=255)                                                   aaGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAg                                                                                       >  2:1285946/1‑150 (MQ=255)                                                    aGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAgg                                                                                      <  1:1383628/150‑1 (MQ=255)                                                     gCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGt                                                                                     >  2:143733/1‑150 (MQ=255)                                                      cAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGtt                                                                                    <  2:1069613/150‑1 (MQ=255)                                                           cGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATAGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCAACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCAt                                                                               <  1:8419/150‑1 (MQ=255)                                                            ggCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGTAATCGCATGTTTGCACGACATCACGGCATGATTGCCAAGCTCGTCACACG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATc                                                                              <  2:85401/150‑1 (MQ=255)                                                            ggCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATc                                                                              <  1:408755/150‑1 (MQ=255)                                                                gCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACa                                                                          <  1:1332979/150‑1 (MQ=255)                                                                 cGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCACGCAGCAAGGAGGGCACTTCATGAGTTGCAATCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACAt                                                                         >  1:739511/1‑150 (MQ=255)                                                                   gCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGc                                                                       >  2:355582/1‑150 (MQ=255)                                                                     aTTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCaa                                                                     <  1:770610/150‑1 (MQ=255)                                                                         ttgCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTAGTCAAGGGCAGACGT‑GGCCACGTCTGCTGTGTCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCgg                                                                 <  1:1242243/150‑1 (MQ=255)                                                                            cccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTg                                                                                   >  2:450521/1‑129 (MQ=255)                                                                            cccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTg                                                                                   <  1:450521/129‑1 (MQ=255)                                                                            cccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCaa                                                                     >  1:506833/1‑143 (MQ=255)                                                                            cccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCaa                                                                     <  2:506833/143‑1 (MQ=255)                                                                             ccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCcgc                                                             >  1:731488/1‑150 (MQ=255)                                                                                 aTGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGtt                                                         <  1:887881/150‑1 (MQ=255)                                                                                    gcgcCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      <  1:611861/150‑1 (MQ=255)                                                                                        cccAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGTTTGCCATCGACATGCAAGCGGCCGCGGTTGTcgcgg                                                  >  1:159478/1‑149 (MQ=255)                                                                                        cccAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      <  1:217364/146‑1 (MQ=255)                                                                                        cccAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      >  2:217364/1‑146 (MQ=255)                                                                                                     cATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCt                                     <  2:932239/150‑1 (MQ=255)                                                                                                     cATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCt                                     >  2:391958/1‑150 (MQ=255)                                                                                                                  cATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCgg                        <  1:1025753/150‑1 (MQ=255)                                                                                                                  cATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCgg                        <  2:1286738/150‑1 (MQ=255)                                                                                                                   aTCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGc                                                                       <  2:1077479/102‑1 (MQ=255)                                                                                                                   aTCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGc                                                                       >  1:1077479/1‑102 (MQ=255)                                                                                                                   aTCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCAACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTCCGGTGAGCggg                       >  1:142149/1‑150 (MQ=255)                                                                                                                      aCGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCACCGCTGATCTGTTCGGTGAGCGGGTTg                    >  2:1053186/1‑150 (MQ=255)                                                                                                                        ggCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGAATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGc                  <  1:1379914/150‑1 (MQ=255)                                                                                                                           aTGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCtt               <  2:1311206/150‑1 (MQ=255)                                                                                                                              aTTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGaa            <  2:405125/150‑1 (MQ=255)                                                                                                                              aTTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGaa            <  1:1285946/150‑1 (MQ=255)                                                                                                                              aTTGCAAAGCACGTCAGAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGaa            <  2:747736/150‑1 (MQ=255)                                                                                                                                   cAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCaa       <  1:143733/150‑1 (MQ=255)                                                                                                                                     aGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTCCGGTGAGCGGGTTGGCCTTGAAGCCAAcc     <  2:142149/150‑1 (MQ=255)                                                                                                                                       cTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAg                                                                                       <  2:1060662/66‑1 (MQ=255)                                                                                                                                       cTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAg                                                                                       >  1:1060662/1‑66 (MQ=255)                                                                                                                                       cTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCGt   >  1:1062945/1‑150 (MQ=255)                                                                                                                                        tCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCgcc  >  2:1077684/1‑148 (MQ=255)                                                                                                                                              gacg‑TGGCCTACGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTg                    >  2:286157/4‑126 (MQ=255)                                                                                                                                              gacg‑TGGCCTACGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTg                    <  1:286157/123‑1 (MQ=255)                                                                                                                                               acg‑TGGCCTACGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGaa            <  2:254241/131‑1 (MQ=255)                                                                                                                                               acg‑TGGCCTACGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGaa            >  1:254241/3‑133 (MQ=255)                                                                                                                                                       |                                                                                                                                        TCCAGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCGTC  >  NC_002947/3268431‑3268715

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