Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,268,583	G→C	6.3%	A291A (GCC→GCG)	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,583	0	G	C	6.3%	199.9 / 4.6	64	A291A (GCC→GCG)	PP_2867	pyridine nucleotide‑disulphide oxidoreductase family protein
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base C (2/2);  total (28/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 = 5.48e-01

AGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCGTC  >  NC_002947/3268434‑3268715                                                                                                                                                       |                                                                                                                                     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                                                                     <  2:506833/143‑1 (MQ=255)                                                                         cccGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCaa                                                                     >  1:506833/1‑143 (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                                                      >  2:217364/1‑146 (MQ=255)                                                                                     cccAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      <  1:217364/146‑1 (MQ=255)                                                                                                  cATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCt                                     <  2:932239/150‑1 (MQ=255)                                                                                                  cATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCt                                     >  2:391958/1‑150 (MQ=255)                                                                                                               cATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCgg                        <  2:1286738/150‑1 (MQ=255)                                                                                                               cATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCgg                        <  1:1025753/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)                                                                                                                                                       ggCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      <  1:1329182/81‑1 (MQ=255)                                                                                                                                                       ggCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTc                                                      >  2:1329182/1‑81 (MQ=255)                                                                                                                                                       |                                                                                                                                     AGGCACGTCACGTACTTGGGCTGGCTGTAAGCCGTTGGTGCTACGCCAAGCAACTCGGCAGCGGCATTGTTGCCCGCATGGCGCCCCAGCGCAATCGCATGTTGGCACGACATCACGGCATGATTGCCAAGCTCGTCACAGG‑TGGCGTAGGCCACGTCTCCAGCTGCATACACGGCATCATTGCCCTTGACCTTCAGGTTGCCATCGACATGCAAGCGGCCGCGGTTGTCGCGCTCGCCGCTGATCTGTTCGGTGAGCGGGTTGGCCTTGAAGCCAACCGTC  >  NC_002947/3268434‑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: ‑