Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,227,465	T→A	6.1%	intergenic (‑21/‑426)	PP_1964 ← / → holB	dGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,227,465	0	T	A	6.1%	205.4 / 4.0	66	intergenic (‑21/‑426)	PP_1964/holB	dGMP/dTMP/dhydroxyCMP kinase/DNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base A (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 = 9.32e-01

CTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCG  >  NC_002947/2227318‑2227609                                                                                                                                                    |                                                                                                                                                 cTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAttt                                                                                                                                                >  1:1000992/1‑150 (MQ=255)  ttGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTc                                                                                                                                               >  1:160002/1‑150 (MQ=255)  ttGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTc                                                                                                                                               >  2:1752426/1‑150 (MQ=255)  ttGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGAGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTc                                                                                                                                               <  1:1188181/150‑1 (MQ=255)     gCAACCTGTTTTGAGCGGATCAGCTAGGGCATATTTCGCGACCTCAGGGTGCTTCATGAGCATGTAAGCGACGGTGTCTTTACCGTATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCTCGAATTTCggg                                                                                                                                            <  2:252972/150‑1 (MQ=255)      cAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCgggg                                                                                                                                           <  1:1657093/150‑1 (MQ=255)       aaCCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggggg                                                                                                                                          >  1:1346684/1‑150 (MQ=255)         ccTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                        <  1:1044020/150‑4 (MQ=255)          cTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggggcctg                                                                                                                                       <  1:1542353/150‑5 (MQ=255)               ttGAGCGGATCAGCTAGGGCATATGTCACAGCCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAggg                                                                                                                                  >  1:1809217/1‑150 (MQ=255)                 gcgCGGGTCAGGTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGAt                                                                                                                                <  1:1485765/148‑1 (MQ=255)                        tCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAAc                                                                                                                         >  2:1215540/1‑150 (MQ=255)                        tCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAAc                                                                                                                         >  2:1215554/1‑150 (MQ=255)                               gggCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGtt                                                                                                                  >  1:1367701/1‑150 (MQ=255)                                 gCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTa                                                                                                                <  2:1247023/150‑1 (MQ=255)                                 gCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTa                                                                                                                >  1:462062/1‑150 (MQ=255)                                     atGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACAtt                                                                                                            >  2:1908931/1‑150 (MQ=255)                                       gTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTCCAGTCt                                                                                                          >  1:329287/1‑150 (MQ=255)                                         cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCttc                                                                                                        <  2:1329774/150‑2 (MQ=255)                                         cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                        <  1:1215554/150‑1 (MQ=255)                                         cgcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTa                                                                                                        <  1:1215540/150‑1 (MQ=255)                                          gcgACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAg                                                                                                       >  2:1185739/1‑150 (MQ=255)                                                    ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTaa                                                                                                    <  2:764500/143‑1 (MQ=255)                                                    ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTaa                                                                                                    >  1:764500/1‑143 (MQ=255)                                                    ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                        <  2:905468/107‑4 (MQ=255)                                                    ggTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAAATTCggggcct                                                                                                                                        >  1:905468/1‑104 (MQ=255)                                                     gTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATg                                                                                                                               <  2:384289/115‑1 (MQ=255)                                                     gTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATg                                                                                                                               >  1:384289/1‑115 (MQ=255)                                                     gTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGta                                                                                            <  2:1346684/150‑1 (MQ=255)                                                      tGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGtat                                                                                           <  1:1822121/150‑1 (MQ=255)                                                      tGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGtat                                                                                           <  1:1821871/150‑1 (MQ=255)                                                             tGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTa                                                                                                                             <  1:496633/109‑1 (MQ=255)                                                             tGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTa                                                                                                                             >  2:496633/1‑109 (MQ=255)                                                                gCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTa                                                                                 >  2:6077/1‑150 (MQ=255)                                                                             ggTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTAc                                                                    <  2:1000992/150‑1 (MQ=255)                                                                                gtCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCa                                                                 <  1:1481687/150‑1 (MQ=255)                                                                                        cGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTTCATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGaa                                                         <  2:329287/150‑1 (MQ=255)                                                                                        cGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACa                                                                         <  2:1710408/134‑1 (MQ=255)                                                                                        cGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACa                                                                         >  1:1710408/1‑134 (MQ=255)                                                                                        cGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGaa                                                         >  1:1670705/1‑150 (MQ=255)                                                                                         ggATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAAt                                                        <  1:1110751/150‑1 (MQ=255)                                                                                           aTCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATAc                                                      <  2:160002/150‑1 (MQ=255)                                                                                             ggAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACaa                                                    <  2:682567/149‑1 (MQ=255)                                                                                             cGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGGAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACaa                                                    <  2:462062/150‑1 (MQ=255)                                                                                               aGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGa                                                  >  2:1306244/1‑150 (MQ=255)                                                                                                 cTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTa                                                                           >  1:170246/1‑123 (MQ=255)                                                                                                 cTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTa                                                                           <  2:170246/123‑1 (MQ=255)                                                                                                  tACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGtat                                                                                           <  2:787261/106‑1 (MQ=255)                                                                                                  tACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGtat                                                                                           >  1:787261/1‑106 (MQ=255)                                                                                                    cTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCa                                             <  1:1752426/150‑1 (MQ=255)                                                                                                     tGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTAAATTTACTGAAGAACGAATACAAGTCGCCAc                                            >  2:111799/1‑150 (MQ=255)                                                                                                           gAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGaaa                                                                                               >  1:1886057/1‑92 (MQ=255)                                                                                                           gAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAAc                                                                                               <  2:1886057/93‑1 (MQ=255)                                                                                                           gAGGCCGATGATGGTTGGCATTTTTCGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACAATCTTAGTAAGGAACGTATAATTGGCGTAGGGTTACAGTTACTCAAGAACGAATACAAGACGCCACCTATGc                                      <  1:1306244/150‑1 (MQ=255)                                                                                                              gCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAg                                                                                      >  2:1172562/1‑99 (MQ=255)                                                                                                              gCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAg                                                                                      <  1:1172562/99‑1 (MQ=255)                                                                                                                cGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                            <  2:648474/43‑1 (MQ=255)                                                                                                                cGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCggg                                                                                                                                            >  1:648474/1‑43 (MQ=255)                                                                                                                   tgatgGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAAGTATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAgg                              <  1:425874/150‑1 (MQ=255)                                                                                                                           ggCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGaa                                                                                                >  1:1928907/1‑76 (MQ=255)                                                                                                                           ggCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGaa                                                                                                <  2:1928907/76‑1 (MQ=255)                                                                                                                                       tctaacccccGAAATTCGGGGGTTAGAGGGATGTAGAGCTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcg          >  1:1572805/7‑150 (MQ=255)                                                                                                                                       tCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcg          >  2:845663/1‑150 (MQ=255)                                                                                                                                         catgCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTgtcggt        <  1:768265/147‑1 (MQ=255)                                                                                                                                               ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCg  >  2:245048/1‑150 (MQ=255)                                                                                                                                               ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCg  >  2:244984/1‑150 (MQ=255)                                                                                                                                               ccGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCg  >  2:1174460/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTTGGCAACCTGCTTTGAGCGGATCAGCTAGGGCATATGTCGCGACCTCAGGGTGCTCCATGAGCATGGAAGCGACGGTGTCTTTACCGGATCGAGCTACTGCGGCGAGGCCGATGATGGTTGGCATTTTTGGCTCCAGGCCCCGAATTTCGGGGGTTAGAGGGATGTAGAACTCATGTTTACATTCTTAGTAAGGAACGTATAATAGGCCTATGTTTACATTTACTCAAGAACGAATACAAGACGCCACCTATGCCTATCAGGGCGGCGGTCGTGGTCTGTCGGTCGAGCG  >  NC_002947/2227318‑2227609

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