Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,587,762	C→G	6.4%	A185A (GCC→GCG)	catA‑II →	catechol 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,587,762	0	C	G	6.4%	187.4 / 4.9	63	A185A (GCC→GCG)	catA‑II	catechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base G (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

CTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCA  >  NC_002947/3587616‑3587911                                                                                                                                                    |                                                                                                                                                      cTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAg                                                                                                                                                    <  1:643613/150‑1 (MQ=255)  tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGGCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                                   >  2:699030/1‑150 (MQ=255)  tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                                   >  1:383449/1‑150 (MQ=255)  tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAgg                                                                                                                                                   >  1:624750/1‑150 (MQ=255)   ggAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATTCCCAggg                                                                                                                                                  <  2:514227/150‑1 (MQ=255)                cTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgcg                                                                                                                                     <  2:339262/150‑1 (MQ=255)                  ggACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgcgcg                                                                                                                                   >  1:1186637/1‑150 (MQ=255)                      ccGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTTTTGGTCCAAAGCCAGTCGGCGTGCCAAGCGCAAGCTGC‑CGTTAATTCCGATGGCCAGGGGCGCTACCGCGCGCGCTcc                                                                                                                               >  1:870904/1‑150 (MQ=255)                           cGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGt                                                                                                                          <  2:624750/150‑1 (MQ=255)                                   cGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCgg                                                                                                                  >  1:972061/1‑150 (MQ=255)                                   cGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCgg                                                                                                                  <  2:1186637/150‑1 (MQ=255)                                        ccGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTAt                                                                                                             >  2:394086/1‑150 (MQ=255)                                             tGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTg                                                                                                        <  2:972061/150‑1 (MQ=255)                                             tGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTg                                                                                                        <  2:887469/150‑1 (MQ=255)                                              gCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGc                                                                                                       <  2:5156/150‑1 (MQ=255)                                                cACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACAAATACCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATAGCTGCGa                                                                                                     >  2:1068848/1‑150 (MQ=255)                                                 aCGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGAc                                                                                                    >  1:1047892/1‑150 (MQ=255)                                                       gACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATACACAGGGGCGCTAATGCGCGCTCGCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAg                                                                                              >  1:1173647/1‑150 (MQ=38)                                                          cTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAgggg                                                                                           <  1:699030/150‑1 (MQ=255)                                                              ggCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCaa                                                                                       <  2:880621/150‑1 (MQ=255)                                                               gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGAc                                                                                                    >  1:466764/1‑136 (MQ=255)                                                               gCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGAc                                                                                                    <  2:466764/136‑1 (MQ=255)                                                                 aTGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAAccc                                                                                    >  1:79237/1‑150 (MQ=255)                                                                   gCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAg                                                                                  <  2:79237/150‑1 (MQ=255)                                                                       ctACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTCCCGCGCGCGCTCCATCGTGCCAACGGGCTATGTCTGCGACCCGTATGGGCCAACCCAGGAAt                                                                              <  2:1103857/149‑1 (MQ=255)                                                                         aCCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                            >  1:802015/1‑150 (MQ=255)                                                                              tGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGa                                                                       <  2:416105/150‑1 (MQ=255)                                                                                 tACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGTTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACct                                                                    >  1:404166/1‑150 (MQ=255)                                                                                    cTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACctgct                                                                 >  2:336479/1‑150 (MQ=255)                                                                                      aCTCGGTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTgg                                                               >  2:998491/1‑150 (MQ=255)                                                                                       cTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTggg                                                              <  1:336479/150‑1 (MQ=255)                                                                                        tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTg                                                                         <  1:988741/138‑1 (MQ=255)                                                                                        tCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTg                                                                         >  2:988741/1‑138 (MQ=255)                                                                                                   ccAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCa                                                  <  2:846197/150‑1 (MQ=255)                                                                                                     aaaGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAgc                                                <  1:998491/150‑1 (MQ=255)                                                                                                       aGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAgcgc                                              >  2:527501/1‑150 (MQ=255)                                                                                                        gCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                            <  1:990730/119‑1 (MQ=255)                                                                                                        gCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGc                                                                            >  2:990730/1‑119 (MQ=255)                                                                                                        gCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGcc                                             >  1:410302/1‑150 (MQ=255)                                                                                                            gTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAggg                                                                                            >  1:502254/1‑99 (MQ=255)                                                                                                            gTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAggg                                                                                            <  2:502254/99‑1 (MQ=255)                                                                                                             tCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCCGCTGGGCCGCCAGTGCCAGCGcccttcg                                        <  2:1168686/150‑5 (MQ=21)                                                                                                                gCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCAc                                     <  1:527501/150‑1 (MQ=255)                                                                                                                  gTACAACCTCC‑GTGGGCGCATTGTTACCGATGCCCAGGGGCGCTACCGCGCGCACTCCATCGTGACATCGGGCTTTGGCTGCGACCCGCAGGGCCCAACCCAGGACTGCCAGGACCTATTGGGCCGTCATGGCCAGCGCCCGGCGCACGt                                   >  2:424585/1‑150 (MQ=18)                                                                                                                     cAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCa                                <  2:1173647/150‑1 (MQ=255)                                                                                                                        ccTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCActt                             >  1:330532/1‑150 (MQ=255)                                                                                                                          tGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCActtct                           <  2:872620/150‑1 (MQ=255)                                                                                                                              ccctgggcATCGGTAACGATGCGCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTa                        >  2:394760/7‑150 (MQ=18)                                                                                                                                ctgggcATCGGTAACGATGCGCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                              <  2:432040/90‑1 (MQ=21)                                                                                                                                ctgggcATCGGTAACGATGCGCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAAt                                                                              >  1:432040/5‑94 (MQ=21)                                                                                                                                   cgcATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCg                   <  2:870904/150‑1 (MQ=255)                                                                                                                                    gcATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCgg                  >  1:103122/1‑150 (MQ=255)                                                                                                                                         gTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGTTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCa             <  2:404166/150‑1 (MQ=255)                                                                                                                                          ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTTCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTCGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg            >  2:293387/1‑150 (MQ=255)                                                                                                                                          ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg            <  2:1047892/150‑1 (MQ=255)                                                                                                                                          ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg            <  2:410302/150‑1 (MQ=255)                                                                                                                                          ttACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTCCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg            <  1:293387/150‑1 (MQ=255)                                                                                                                                          gtaacGATGCGCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAg            <  1:394760/146‑1 (MQ=255)                                                                                                                                              cGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAggg          >  2:1123459/1‑148 (MQ=255)                                                                                                                                              cGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAggg          <  1:1123459/148‑1 (MQ=255)                                                                                                                                              cGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTa        >  1:1006520/1‑150 (MQ=255)                                                                                                                                              cGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTa        >  2:228778/1‑150 (MQ=255)                                                                                                                                              cGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTa        >  1:165980/1‑150 (MQ=255)                                                                                                                                                    ccAGGTGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCa  >  2:689424/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACGGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGC‑GTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGCGACCCGCAGGGGCCAACCCAGGAATGCCTGGACCTGCTGGGCCGTCATGGCCAGCGCCCGGCGCACGTGCACTTCTTTATCTCGGCCCCAGGGTACCGGCA  >  NC_002947/3587616‑3587911

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