Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	883,866	A→C	15.5%	intergenic (+101/‑227)	PP_0766 → / → PP_0767	hypothetical protein/LuxR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	883,866	0	A	C	15.5%	134.3 / 2.8	59	intergenic (+101/‑227)	PP_0766/PP_0767	hypothetical protein/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/16);  new base C (3/7);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 7.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCCATCAT  >  NC_002947/883717‑884015                                                                                                                                                      |                                                                                                                                                       cAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACa                                                                                                                                                        <  2:1226129/150‑1 (MQ=255)      ccGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTc                                                                                                                                                   <  1:665258/150‑1 (MQ=255)          tttCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGACCGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCcgcg                                                                                                                                               <  2:1324271/150‑1 (MQ=255)          tttCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCcgcg                                                                                                                                               <  2:3218189/150‑1 (MQ=255)                  ccGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGc                                                                                                                                       <  2:3097305/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  1:1137734/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  2:3268446/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  2:267838/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  2:1159410/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  1:4068520/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  1:2654314/150‑1 (MQ=255)                    gCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCgg                                                                                                                                     <  1:2239033/150‑1 (MQ=255)                       gcgcTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGtt                                                                                                                                  <  2:3408974/150‑1 (MQ=255)                          cTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGTACCCGCGCCGACAGGTTCCGCGGGGGGGGGGGGTTCAc                                                                                                                               >  2:524543/1‑150 (MQ=255)                          cTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACCGCGGAACCTGTAGGAGCGGGTTCAc                                                                                                                               <  2:3336852/150‑1 (MQ=25)                          cTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTAGCGGGGGGAACCCCACCGACACGTTACGCGGTGGGGGCggggtgac                                                                                                                               >  2:524542/1‑145 (MQ=255)                              gTCCGTCGGGTATGCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGGGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACCGGTTCCGCGGTAGGAGCGGGTTCACCcgc                                                                                                                           <  1:856329/150‑1 (MQ=255)                                    cGGGTATTCGCTAAGCTGAAATCACTGCACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTATATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCGCCTCCCGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAagag                                                                                                                     <  2:310639/150‑1 (MQ=39)                                       gTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAAACCGCACCTACAGGTTGCGCGGTAGGCGCGGGTTCACCCCCGAAGAgggc                                                                                                                  >  2:3176484/1‑148 (MQ=255)                                         ttgcGCTACGCTGAACGCACTGAACGCCGTGTGATCCCCAGACCGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTGCGGGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGTTTCCCCCGCGAAGAGGCCgg                                                                                                                <  2:1013201/147‑1 (MQ=255)                                             gCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGggcaa                                                                                                            <  2:3730104/150‑1 (MQ=255)                                             gCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGGGAACCCGCGCCTACCGGGTCCGCGGTAGGAGCGGGGTCACCCGCGAAGAGGCCGggcaa                                                                                                            >  2:3645280/1‑150 (MQ=38)                                                      aaaTCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACCGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCaa                                                                                                   <  2:2232868/150‑1 (MQ=255)                                                          cACTGACCCCCGTGTGATCCCCAGACCGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTCACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACCGGTTCCGGGGTAGGAGGGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAAcctc                                                                                               <  2:3209496/150‑1 (MQ=38)                                                           aCTGAACCCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAAcctcc                                                                                              >  2:4415538/1‑150 (MQ=255)                                                            ctgacggCCGTGTTATCCCCCCACCGCCCGGGCGCGGGGGGGCTTTTATTTTCCCTGCCCGGGCCCCTGGGGGGGTTACCCGGCGCCTCCTGGTGCGGGGTTTGTCGCGGTTACACCCGCGAAGAGGCCGGGCAAGGCAATCAAcctcct                                                                                             <  2:393806/144‑1 (MQ=255)                                                             tGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGCAGAGGCCGGGCAAGGCAAACAAcctccag                                                                                            >  2:4108123/1‑148 (MQ=255)                                                             tGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTg                                                                                            >  2:4119138/1‑150 (MQ=255)                                                                  gcggtgtGAGCCCCACAACGCCCGGGCGCGCGGGGGGGTTTCTATTTCCTGACCGGGCCTCTTCGGGGGTGAACCCGCTCCTACAGTTTCCGCGGTAGGAGGGGTTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATAc                                                                                       <  2:1253985/147‑1 (MQ=37)                                                                                        ccGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGGGGAACCCCCCCCACCCGTTGCGCGGGGGGGGCCGGGTCACCCCCCCAGGGGCCCGGCCAAGCCAGCAAACACCCGGAAACCCCGCCCAAAcccccccccccc                                                                 >  2:855630/1‑148 (MQ=255)                                                                                                       gTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGGAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCACTCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCCTCAGcc                                                 >  2:4072503/1‑150 (MQ=255)                                                                                                             aTTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCCGGTTCACCCGCGAAGAGGCCGGGCAAGGCCCTCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCAACCGCCaaggtg                                           >  2:2612010/1‑145 (MQ=38)                                                                                                              tttCCTGACCCGGCCTCTTCGCGGGTGAACCCCCTCCGACAGGGGGCGCGGGAGGAGCGGGGTAACCCCCGACCCCGCCGGGCCAGGCCAACAACCACCCGGATAACGA‑GCAAAACCCCCCCCCCCCCCCCCCTCCAACAGcccgggtgg                                          >  2:4377811/1‑143 (MQ=255)                                                                                                                tCCTGACCCGGCCTCTTCGCGGGTGAACCCCCTACTAAAGGTTCCCCGGTAGGAGCGGGGTCACCCGCGAAGAGGCCCGGCCAGGCCATCAAACACCTGCAGAACGC‑TCCAAACCCCCGCCCCAGCCCCAGTCCAGCCGCCAAGTTgggt                                        >  2:333430/1‑147 (MQ=21)                                                                                                                   tGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACCGGTTCCGCGGTAGGAGGGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCttttt                                     <  2:2146146/150‑1 (MQ=255)                                                                                                                         ggCATCGTCGGGGGTTACCCCGCTCCTACGGGTGCGGGGTTAGGCGGGGTTCCCCCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCa                               <  1:1661113/150‑1 (MQ=255)                                                                                                                                 cgcgGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTc                       <  2:1320537/150‑1 (MQ=255)                                                                                                                                    gggTGAACCCGCTCCTACCGCGGAACCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCcgc                    >  2:930927/1‑150 (MQ=25)                                                                                                                                     ggTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCcgcg                   >  1:1187184/1‑150 (MQ=255)                                                                                                                                          aCCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAt              >  1:690153/1‑150 (MQ=255)                                                                                                                                          aCCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAt              >  1:690264/1‑150 (MQ=255)                                                                                                                                           cccGCTCCTACCGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAtt             <  2:2936465/150‑1 (MQ=255)                                                                                                                                           cccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAtt             >  1:832913/1‑150 (MQ=255)                                                                                                                                           cccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAtt             >  1:2392836/1‑150 (MQ=255)                                                                                                                                           cccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTAtt             >  1:2812693/1‑150 (MQ=255)                                                                                                                                            ccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTg            >  2:3475368/1‑150 (MQ=255)                                                                                                                                            ccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTg            >  2:3836865/1‑150 (MQ=255)                                                                                                                                            ccGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTg            >  1:3581274/1‑150 (MQ=255)                                                                                                                                             cGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTgc           >  2:808855/1‑150 (MQ=255)                                                                                                                                             cGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTgc           >  2:808747/1‑150 (MQ=255)                                                                                                                                             cGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGGCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCGATTgc           >  1:2717984/1‑150 (MQ=255)                                                                                                                                             cGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACACCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTATTTGTCTCAGTCTTTTCCGCGCCTATTgc           >  2:1522203/1‑150 (MQ=255)                                                                                                                                              gCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATGCCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTGTTTGTCGCAGTCTTTTCCGCGCCTATTgcg          >  2:3384440/1‑150 (MQ=255)                                                                                                                                              gCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTgcg          >  2:3384487/1‑150 (MQ=255)                                                                                                                                                tCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGcc        >  1:2501108/1‑150 (MQ=255)                                                                                                                                                tCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGcc        >  1:2501080/1‑150 (MQ=255)                                                                                                                                                tCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGcc        >  2:1779553/1‑150 (MQ=255)                                                                                                                                                   tACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCACTCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCcat     >  2:3909313/1‑150 (MQ=255)                                                                                                                                                   tACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCcat     >  1:172549/1‑150 (MQ=255)                                                                                                                                                    aCAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCcatc    >  1:3441875/1‑150 (MQ=255)                                                                                                                                                      aGGTTCCGTGGTAGGAGTGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCcctcat  >  2:1941248/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CAAGGCCGATTTCCAGCCCGCAGCGCTGCGTCAGTCGGGTATTCGCTAAGCTGAAATCACTGAACGCCGTGTGATCCCCAGAACGCCCCGGCTCGCCGGGGCGTTTCTATTTCCTGACCCGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGTTCCGCGGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCAAGGCAATCAACCTCCTGAATACCTC‑TCCAAACCCCCGCCCCAGCCCCTGTCCATCAGCCATGTTGCTTTTTGTCGCAGTCTTTTCCGCGCCTATTGCGCCCATCAT  >  NC_002947/883717‑884015

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