Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,005,021	G→C	6.1%	intergenic (+59/‑117)	PP_5259 → / → ydcJ	LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,005,021	0	G	C	6.1%	199.6 / 2.9	66	intergenic (+59/‑117)	PP_5259/ydcJ	LysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

GCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGA  >  NC_002947/6004873‑6005166                                                                                                                                                     |                                                                                                                                                    gcgcATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAAGGGAGTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                                                     <  5:271675/149‑1 (MQ=255)     aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTTAAGCCGCTCCTAAGGGGACCGCGTAAGCGCGGCTTTAGGCGGGa                                                                                                                                                   >  6:223654/1‑147 (MQ=255)     aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGa                                                                                                                                                   >  6:214792/1‑147 (MQ=255)      ttACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGTCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                                                <  7:141399/149‑1 (MQ=255)      ttACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaa                                                                                                                                                  >  6:311033/1‑147 (MQ=25)        aCCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAgg                                                                                                                                               <  7:231052/148‑1 (MQ=255)        aCCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaagag                                                                                                                                               >  8:167655/1‑145 (MQ=25)        aCCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaagag                                                                                                                                               <  7:167655/148‑4 (MQ=25)           tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                                           >  7:255796/1‑149 (MQ=255)             gCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTc                                                                                                                                         >  1:256948/1‑149 (MQ=255)             gCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGCCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCg                                                                                                                                           <  5:311033/147‑1 (MQ=18)                     cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                                              <  8:334860/136‑1 (MQ=255)                     cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                                              >  7:334860/1‑136 (MQ=255)                     cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGa                                                                                                                                 <  2:256948/149‑1 (MQ=255)                             cGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                                                                         <  5:214792/149‑1 (MQ=255)                              ggCAGAGGTGCCAAAGGTGCGGGCGCTGTTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                                                                        <  4:48593/149‑1 (MQ=255)                                     gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGCCCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGa                                                                                                                 >  2:84464/1‑148 (MQ=255)                                     gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGa                                                                                                                 <  6:161264/148‑1 (MQ=255)                                      tGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGaca                                                                                                               >  1:87973/1‑149 (MQ=255)                                        ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                                                    <  5:211880/127‑1 (MQ=255)                                        ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCTCCTACGGGGACCGCGAAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                                                    >  6:211880/1‑127 (MQ=255)                                                 cGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACtttt                                                                                                    <  8:255796/149‑1 (MQ=255)                                                   ggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACtttttt                                                                                                  >  8:214751/1‑149 (MQ=255)                                                     cgcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTc                                                                                                <  8:233888/149‑1 (MQ=255)                                                      gcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCg                                                                                               >  8:253419/1‑149 (MQ=255)                                                          gtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCtt                                                                                            <  5:223654/148‑1 (MQ=255)                                                                  gATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCa                                                                                     >  8:81969/1‑147 (MQ=255)                                                                  gATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCa                                                                                     <  7:81969/147‑1 (MQ=255)                                                                       gCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTAAGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAATGGCCCTCGACGTCGATCAAAAAAGGAA‑TGACACACTACCTTTTTTTCGCTTGTGAGCaaatagc                                                                               >  6:192579/1‑144 (MQ=255)                                                                          aggaGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGGAAAGTGCATTc                                                                           <  2:214938/146‑1 (MQ=255)                                                                          aGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                                                >  5:167930/1‑81 (MQ=255)                                                                          aGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                                                <  6:167930/81‑1 (MQ=255)                                                                          aGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTc                                                                           <  7:214751/149‑1 (MQ=255)                                                                                  cAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑t                                                                                                                   >  7:109436/1‑101 (MQ=255)                                                                                  cAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑t                                                                                                                   <  8:109436/101‑1 (MQ=255)                                                                                  cAGGTTTGACGGCCTCTTCGCGGCTAAAACCGCTCCTACGGGGACCGCGGAGGCCCGGCTTTAGCCGGGAAAGGGCCGTAGAGGTCCATCAAAAAAGGAA‑TGACACACTCACTTATTTTAGCATGTGAGCAAAGTGCATTCCCAGCttt                                                                   >  2:198659/1‑149 (MQ=255)                                                                                    ggTTTGACGGCCTGTTCGCGGCTAAAGCCTGTCCTACGGGGACCGGGTAGGCGGGGCTTTATCCGGGAAATGGGCGTCGAGGTCGATCAAAAAAGGAATTCACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTc                                                                  <  1:84464/149‑1 (MQ=255)                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGGGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                                                                <  5:192579/149‑1 (MQ=255)                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                                                                >  1:207402/1‑149 (MQ=255)                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                                                                <  7:253419/149‑1 (MQ=255)                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCAt                                                                <  6:83181/149‑1 (MQ=255)                                                                                     gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTAGCTTGTGAGAAAAGTGCATTACCAGCTTTCAt                                                                >  8:100985/1‑149 (MQ=255)                                                                                          aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTc                                                                  <  7:73632/142‑1 (MQ=255)                                                                                          aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                                            <  1:198659/148‑1 (MQ=255)                                                                                          aCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTAACTTTTTTTCGGTTTTGAGCAAAGTGCATTCCCAGCTTTc                                                                  >  8:73632/1‑142 (MQ=255)                                                                                                tctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTcc                                                     >  7:73327/1‑148 (MQ=255)                                                                                                  ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGacaca                                                                                                             <  4:301955/91‑1 (MQ=255)                                                                                                  ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGacaca                                                                                                             >  3:301955/1‑91 (MQ=255)                                                                                                      cgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCa                                              >  4:193638/1‑149 (MQ=255)                                                                                                        gCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAAc                                            <  2:207402/149‑1 (MQ=255)                                                                                                              gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGaa                                      <  8:73327/149‑1 (MQ=255)                                                                                                               ccGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAg                                     >  1:336256/1‑149 (MQ=255)                                                                                                               ccGCTCCCAAGGGGAACGGGTAGGGGGGGGTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAg                                     <  3:5141/149‑1 (MQ=255)                                                                                                                cgcgCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCACTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGc                                    >  2:261838/1‑149 (MQ=38)                                                                                                                 gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACttttttt                                                                                                 >  5:138428/4‑88 (MQ=25)                                                                                                                 gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACttttttt                                                                                                 <  6:138428/85‑1 (MQ=25)                                                                                                                     cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGACCCCAACCAGGAAGCTAGc                                >  2:69006/1‑148 (MQ=255)                                                                                                                      tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTGACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGa                              >  6:263364/1‑149 (MQ=255)                                                                                                                           ggACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAAcc                                                >  3:248295/1‑126 (MQ=255)                                                                                                                           ggACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAAcc                                                <  4:248295/126‑1 (MQ=255)                                                                                                                              ccGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGc                       >  7:158114/1‑148 (MQ=255)                                                                                                                                gcgTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTAACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTATCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCaa                    <  5:263364/149‑1 (MQ=255)                                                                                                                                   tAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGa                 <  7:100985/149‑1 (MQ=255)                                                                                                                                     ggCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAACGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGAt                <  1:69006/149‑1 (MQ=255)                                                                                                                                       agcgGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCACTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTc             <  1:261838/148‑1 (MQ=255)                                                                                                                                                  ccGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGa  <  3:193638/149‑1 (MQ=255)                                                                                                                                                  ccGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGa  >  1:208055/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAA‑TGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAA‑GCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCCAACGATTTCGTCAGCCCAGA  >  NC_002947/6004873‑6005166

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