Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,848,831	G→C	6.1%	G77A (GGC→GCC)	cysM →	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,848,831	0	G	C	6.1%	207.9 / 4.6	66	G77A (GGC→GCC)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (2/2);  total (32/34)
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.44e-01

CAGCAATACCCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGAT  >  NC_002947/1848682‑1848969                                                                                                                                                      |                                                                                                                                           cAGCAATACCCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACggg                                                                                                                                            >  1:654996/1‑150 (MQ=255)     aaTACCCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGATAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATc                                                                                                                                        <  1:513717/150‑1 (MQ=255)          cctcctGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCTAAGCCACCTCCGGCAACACGGGCATCGCCCt                                                                                                                                   >  2:1073196/1‑150 (MQ=255)           ctcctGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGTCGATACCCTGATCGAATCCACCTCCTGCAACACGGGCATCGCCCtc                                                                                                                                  >  2:654770/1‑149 (MQ=255)                gCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGAt                                                                                                                             <  1:284019/150‑1 (MQ=255)                   cAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggc                                                                                                                          <  1:893488/150‑1 (MQ=255)                   cAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATggc                                                                                                                          >  1:643097/1‑150 (MQ=255)                    aaGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATgccc                                                                                                                         <  1:863941/150‑4 (MQ=255)                          gAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGgcggcggcg                                                                                                                   >  2:581989/1‑150 (MQ=255)                                 aCAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAgg                                                                                                            >  1:762075/1‑150 (MQ=255)                                        tGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACaa                                                                                                     >  1:765236/1‑150 (MQ=255)                                        tGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACaa                                                                                                     <  2:643097/150‑1 (MQ=255)                                        tGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACaa                                                                                                     >  1:622882/1‑150 (MQ=255)                                          ccGGCTCGGTGAACGACCGCCCGGCCATGTCGATGATCGCGCGCGCCGAATTGCGGGGCCAGATTGAACAACGGGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAga                                                                                                   <  1:247289/150‑1 (MQ=255)                                          ccGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAga                                                                                                   <  1:137844/150‑1 (MQ=255)                                            ggCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAgatg                                                                                                 <  1:1297117/150‑1 (MQ=255)                                               tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATc                                                                                              <  1:581989/150‑1 (MQ=255)                                               tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgccc                                                                                                                         >  2:316112/1‑120 (MQ=255)                                               tCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgccc                                                                                                                         <  1:316112/123‑4 (MQ=255)                                                  gTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTg                                                                                           >  2:1348368/1‑150 (MQ=255)                                                    gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGCTACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGAt                                                                                         >  2:462931/1‑150 (MQ=255)                                                    gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGAt                                                                                         <  2:654996/150‑1 (MQ=255)                                                    gAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGAt                                                                                         >  1:643614/1‑150 (MQ=255)                                                     aaaGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATc                                                                                        >  1:415365/1‑150 (MQ=255)                                                     aaaGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                       >  2:634384/1‑114 (MQ=255)                                                     aaaGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGCCATCGCCAGGGCGATgcccgt                                                                                                                       <  1:634384/119‑6 (MQ=255)                                                         aCCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAggg                                                                                                           >  2:185213/1‑127 (MQ=255)                                                         aCCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAggg                                                                                                           <  1:185213/127‑1 (MQ=255)                                                                ggCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACaa                                                                             <  2:622882/150‑1 (MQ=255)                                                                   ccTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTc                                                                          >  2:771145/1‑150 (MQ=255)                                                                   ccTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTc                                                                          >  1:509961/1‑150 (MQ=255)                                                                     tGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCa                                                                        <  1:1348368/150‑1 (MQ=255)                                                                       tCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCAcc                                                                      <  2:1323469/150‑1 (MQ=255)                                                                       tCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCAcc                                                                      >  2:1337786/1‑150 (MQ=255)                                                                                    gcgcCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAgg                                                         >  2:9576/1‑150 (MQ=255)                                                                                          aaTTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCa                                                   <  2:509961/150‑1 (MQ=255)                                                                                                ggggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGTGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACgg                                             <  1:1294195/150‑1 (MQ=255)                                                                                                ggggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAgc                                                               <  2:1267069/132‑1 (MQ=255)                                                                                                ggggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAgc                                                               >  1:1267069/1‑132 (MQ=255)                                                                                                ggggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACgg                                             >  1:519386/1‑150 (MQ=255)                                                                                                 gggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGc                                            <  2:854582/150‑1 (MQ=255)                                                                                                 gggCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGc                                            <  2:128421/150‑1 (MQ=255)                                                                                                        aTCAAGCCAGGCGATACCCTGCTCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGCCCATGCCCGACAACTCCACCGCCGAGCGCAACGCCGCCATGACGGCCTATGgc                                     <  1:1135028/150‑1 (MQ=255)                                                                                                         tCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGgcg                                    >  2:532956/1‑150 (MQ=255)                                                                                                           aaGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGcc                                  <  2:670202/150‑1 (MQ=255)                                                                                                               cAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGc                              >  1:769746/1‑150 (MQ=255)                                                                                                                aGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAgatga                                                                                                <  2:1272814/83‑1 (MQ=255)                                                                                                                aGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAgatga                                                                                                >  1:1272814/1‑83 (MQ=255)                                                                                                                 ggCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTg                            >  1:1344949/1‑150 (MQ=255)                                                                                                                  gCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCAGAGCTGa                           <  2:1344949/150‑1 (MQ=255)                                                                                                                    aaTACCCTCATCGCAGCCACCTCCGGCAACCCGGGCATCGCGCTGGCGATGTCGGCGGCGATCAAGGGCTACAAGATGATCCTCATTATGCCCGAAACCTCCACCGCCGAGAGCAATGCCGCCATCACGGCCTATGGCTCCGAGCTGATc                         <  2:797587/149‑1 (MQ=255)                                                                                                                     aTACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATcc                        <  1:532956/150‑1 (MQ=255)                                                                                                                        cccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTgg                     <  2:643614/150‑1 (MQ=255)                                                                                                                         ccTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGt                    >  1:452816/1‑150 (MQ=255)                                                                                                                           tGATCGAAGCCACCTCCGGCAACACTGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAgatga                                                                                                >  2:380977/1‑72 (MQ=255)                                                                                                                           tGATCGAAGCCACCTCCGGCAACACTGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAgatga                                                                                                <  1:380977/72‑1 (MQ=255)                                                                                                                           tGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGa                  >  1:939487/1‑150 (MQ=255)                                                                                                                            gATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGa                                                                 <  2:483610/102‑1 (MQ=255)                                                                                                                            gATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGa                                                                 >  1:483610/1‑102 (MQ=255)                                                                                                                               cGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCaa              <  1:654770/150‑1 (MQ=255)                                                                                                                                 aaGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAACGGCTACAAGATGATCCTAATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAagg            <  2:1404803/150‑1 (MQ=255)                                                                                                                                     cACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAaggagga        >  2:55548/1‑150 (MQ=255)                                                                                                                                        cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCt                             >  1:422330/1‑126 (MQ=255)                                                                                                                                        cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCt                             <  2:422330/126‑1 (MQ=255)                                                                                                                                        cTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAggg     <  2:939487/150‑1 (MQ=255)                                                                                                                                          ccGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGa   >  1:320794/1‑150 (MQ=255)                                                                                                                                           cGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGAt  >  1:1191442/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                           CAGCAATACCCTCCTGCTCAAGCTCGAAGGTAACAATCCTGCCGGCTCGGTGAAAGACCGCCCGGCCCTGTCGATGATCGCTCGCGCCGAATTGCGGGGCCAGATCAAGCCAGGCGATACCCTGATCGAAGCCACCTCCGGCAACACGGGCATCGCCCTGGCGATGGCGGCGGCGATCAAGGGCTACAAGATGATCCTGATCATGCCCGACAACTCCACCGCCGAGCGCAAGGCCGCCATGACGGCCTATGGCGCCGAGCTGATCCTGGTGACCAAGGAGGAGGGGAT  >  NC_002947/1848682‑1848969

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