Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	268,744	T→A	7.2%	intergenic (‑75/+25)	PP_0216 ← / ← desA	sensory box/GGDEF family protein/delta‑9 fatty acid desaturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	268,744	0	T	A	7.2%	201.0 / 6.8	70	intergenic (‑75/+25)	PP_0216/desA	sensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base A (2/3);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

CGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGG  >  NC_002947/268597‑268889                                                                                                                                                    |                                                                                                                                                   cGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgg                                                                                                                                                  <  1:78017/150‑1 (MQ=255)   gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                <  1:335042/150‑1 (MQ=255)    gtaaaaTCGTTGGAAGGTTCTCGGCCAGGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggggg                                                                                                                                               <  2:627079/148‑1 (MQ=255)        aaTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGt                                                                                                                                           >  2:2204597/1‑150 (MQ=255)        aaTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGt                                                                                                                                           <  2:2003438/150‑1 (MQ=255)        aaTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGt                                                                                                                                           <  1:2013277/150‑1 (MQ=255)              tGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGATAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGAt                                                                                                                                     >  2:1744444/1‑150 (MQ=255)                gAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCa                                                                                                                                   >  1:1471491/1‑150 (MQ=255)                     ttCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggc                                                                                                                              >  2:1847021/1‑150 (MQ=255)                     ttCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggc                                                                                                                              >  2:1846866/1‑150 (MQ=255)                      tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                             <  1:1106107/150‑1 (MQ=255)                      tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                             <  1:1643751/150‑1 (MQ=255)                                gtgCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGcc                                                                                                                   <  1:13479/150‑1 (MQ=255)                                 tgCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCg                                                                                                                  <  1:1923321/150‑1 (MQ=255)                                   cACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGa                                                                                                                                      >  1:204543/1‑128 (MQ=255)                                   cACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGa                                                                                                                                      <  2:204543/128‑1 (MQ=255)                                   cACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGt                                                                                                                >  2:648696/1‑150 (MQ=255)                                    aCGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                                 <  1:537501/116‑1 (MQ=255)                                    aCGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                                 >  2:537501/1‑116 (MQ=255)                                             tcagcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCAAGGAGGTGATATTCCGGATTCGAGTGGAGAGGGCGCGTCGAAGTTCCGAATGCCACGCCCCGCTGGGGGGCGTGCCGATCAGCGGCGGTCAGGCGCCGGTCGGGCGCAGc                                                                                                      <  1:772591/146‑1 (MQ=255)                                             gtttcAACGGGAAAGCGTGCGGCCAGCAACCGGTCCGTCACGCAGGCGATATTCGGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAAGGGCAGGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGGGCCGGTCGGGCGCAGc                                                                                                      <  2:1046129/147‑1 (MQ=255)                                             gTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGCGGTGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGc                                                                                                      >  1:162728/1‑150 (MQ=255)                                              tctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGGAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCg                                                                                                     >  2:548906/1‑150 (MQ=255)                                              tctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCggg                                                                                                            >  1:1803406/1‑143 (MQ=255)                                              tctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCggg                                                                                                            <  2:1803406/143‑1 (MQ=255)                                              tctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCg                                                                                                     >  1:1637425/1‑150 (MQ=255)                                               ctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGa                                                                                                    <  1:2204597/150‑1 (MQ=255)                                                   aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                <  1:1847021/150‑1 (MQ=255)                                                   aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                >  1:756715/1‑150 (MQ=255)                                                   aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                <  1:1846866/150‑1 (MQ=255)                                                     gTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGt                                                                                              <  1:1837642/150‑1 (MQ=255)                                                          aGCGTGCGGCCCGCATCCGGGCTGCCACTAATGCGATATTCCGGTTTCGAGTGGAGAGGTGCGGGCGAAGTTCCGATTTGGACGCCCCGCTTGGGGGCGTGGCGAGCAACGGCGTTCAGGCGTCGGTCTGGCGCAGCGAGTAAGTCTTGa                                                                                         <  2:1136923/150‑1 (MQ=255)                                                          aGCGTGCGGCCAGCATCCGGCCCGCCACACAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                         >  1:768557/1‑150 (MQ=255)                                                          aGCGTGCGGCCAGCATCCGGCCCGCCACACAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                         <  2:768557/150‑1 (MQ=255)                                                               gCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGGAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTtgt                                                                                    <  1:548906/150‑1 (MQ=255)                                                                   ccAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATGCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGACGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctg                                                                                >  1:1952852/1‑150 (MQ=255)                                                                           cGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTcg                                                                       <  1:1707175/150‑1 (MQ=255)                                                                             gCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTcgcg                                                                     <  2:604146/150‑1 (MQ=255)                                                                               ccGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCa                                                                   >  1:1996264/1‑150 (MQ=255)                                                                                 gCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGt                                                                 <  1:2408793/150‑1 (MQ=255)                                                                                  ccACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGGGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTg                                                                >  1:2298429/1‑150 (MQ=255)                                                                                     cGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCt                                                             <  1:2319085/150‑1 (MQ=255)                                                                                         ggCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGAt                                                         <  1:1716531/150‑1 (MQ=255)                                                                                          gCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCGCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATg                                                        >  2:944889/1‑150 (MQ=255)                                                                                                ttCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTg                                                  <  1:648696/150‑1 (MQ=255)                                                                                                ttCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTg                                                  >  2:1030606/1‑150 (MQ=255)                                                                                                  ccGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGt                                                <  2:2097381/150‑1 (MQ=255)                                                                                                   cGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGtt                                               >  1:1560465/1‑150 (MQ=255)                                                                                                     gATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTc                                             <  2:1780270/150‑1 (MQ=255)                                                                                                        tCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGc                                          <  1:1744444/150‑1 (MQ=255)                                                                                                        tCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGc                                          <  2:2298429/150‑1 (MQ=255)                                                                                                             tGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCAt                                     >  1:1822485/1‑150 (MQ=255)                                                                                                               gagaGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGt                                   >  1:1985707/1‑150 (MQ=255)                                                                                                                  agGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTAcc                                >  1:1412384/1‑150 (MQ=255)                                                                                                                    gTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCa                              <  2:1985707/150‑1 (MQ=255)                                                                                                                     tGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAg                             <  1:192915/150‑1 (MQ=255)                                                                                                                      gggggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGt                            >  2:1321765/1‑150 (MQ=255)                                                                                                                        gggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCt                          <  1:1790076/150‑1 (MQ=255)                                                                                                                            gAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCt                                        <  1:381248/132‑1 (MQ=255)                                                                                                                            gAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCt                                        >  2:381248/1‑132 (MQ=255)                                                                                                                                      tcggCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAg            >  1:1205354/3‑150 (MQ=255)                                                                                                                                       tGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGc           <  1:1050052/150‑1 (MQ=255)                                                                                                                                         gCACGCCCCGCTGGGGGGCGTGCCGATCCAGGGCCGGCCGGCGCGGTGTGTGGGTAGCGAGTCAGTCTTGAGCTGAAGAGCACACGTC‑TGAACTCCCGTGATGCCGCTGGTTTCGGGCTCATGTACCGAGTCGTTCATGTCTGCGAgggg         >  1:1151080/1‑147 (MQ=255)                                                                                                                                         gCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCa                      <  2:327024/137‑1 (MQ=255)                                                                                                                                         gCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCa                      >  1:327024/1‑137 (MQ=255)                                                                                                                                           aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGt       <  2:1471491/150‑1 (MQ=255)                                                                                                                                           aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGt       >  2:1429555/1‑150 (MQ=255)                                                                                                                                               cccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTTGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg   >  2:225129/1‑150 (MQ=255)                                                                                                                                               cccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg   >  2:554539/1‑150 (MQ=255)                                                                                                                                                ccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCAGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTgg  >  2:299057/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGC‑GAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGG  >  NC_002947/268597‑268889

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