Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	415,540	T→C	5.6%	L102S (TTG→TCG)	waaC →	ADP‑heptose‑‑LPS heptosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	415,540	0	T	C	5.6%	207.0 / 2.4	71	L102S (TTG→TCG)	waaC	ADP‑heptose‑‑LPS heptosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CCAGGTCATCCCGGTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATGGCCCTGG  >  NC_002947/415391‑415689                                                                                                                                                      |                                                                                                                                                         ccAGGTCATCCCGGTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTTAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGtt                                                                                                                                                          >  1:436653/1‑150 (MQ=255)    ggTCATCCCGGTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGAc                                                                                                                                                       <  1:1856243/150‑1 (MQ=255)     gTCATCCCGGTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGAcc                                                                                                                                                      >  2:668727/1‑150 (MQ=255)             ggTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGt                                                                                                                                              <  2:273618/150‑1 (MQ=255)                ggCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCACGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGATTAAACCGTTagctcaa                                                                                                                                           <  2:1529007/150‑7 (MQ=255)                                 gcAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGcc                                                                                                                       >  1:870308/1‑150 (MQ=255)                                 gcAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCATTCAATAATGGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGACTCTAGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGcc                                                                                                                       >  1:1321015/1‑150 (MQ=255)                                  cAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCt                                                                                                                      <  1:1656469/150‑1 (MQ=255)                                     gAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACCCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGa                                                                                                                   >  2:201424/1‑150 (MQ=255)                                     gAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAACGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGa                                                                                                                   <  1:31380/150‑1 (MQ=255)                                       aCCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑c                                                                                                                                   <  1:82461/132‑1 (MQ=255)                                       aCCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑c                                                                                                                                   >  2:82461/1‑132 (MQ=255)                                              gCAAACCATCAAGAGAGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTc                                                                                                          >  2:1041551/1‑150 (MQ=255)                                                   ccATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑c                                                                                                                                   <  1:365420/120‑1 (MQ=255)                                                   ccATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑c                                                                                                                                   >  2:365420/1‑120 (MQ=255)                                                        aaGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGaa                                                                                                >  1:1038548/1‑150 (MQ=255)                                                            gTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCt                                                                                            >  1:848191/1‑150 (MQ=255)                                                             tGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTg                                                                                           >  2:1637548/1‑150 (MQ=255)                                                              ggCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGAACTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTgg                                                                                          <  1:1391327/150‑1 (MQ=255)                                                                    tGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCGAGc                                                                                    >  1:653769/1‑150 (MQ=255)                                                                      gAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGGACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCg                                                                                  >  1:535899/1‑150 (MQ=255)                                                                       aaGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGc                                                                                 <  2:878849/150‑1 (MQ=255)                                                                       aaGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTCGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGc                                                                                                               >  1:1757941/1‑120 (MQ=255)                                                                       aaGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTCGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGc                                                                                                               >  1:1757804/1‑120 (MQ=255)                                                                       aaGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTCGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGc                                                                                                               <  2:1757804/120‑1 (MQ=255)                                                                       aaGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTCGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGc                                                                                                               <  2:1757941/120‑1 (MQ=255)                                                                        aGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCt                                                                                >  1:1727959/1‑150 (MQ=255)                                                                          gCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTAAAGGC‑‑A‑CCCGTACCCGGCCGGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTc                                                                              >  1:785391/1‑150 (MQ=255)                                                                            gTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTa                                                                            >  1:659144/1‑150 (MQ=255)                                                                            gTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTa                                                                            >  2:664994/1‑150 (MQ=255)                                                                                aaGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTAcgac                                                                        <  2:848191/150‑1 (MQ=255)                                                                                  gcagcGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCNGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCAGCCCGTGTCCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTAcgaccg                                                                      >  1:119584/1‑150 (MQ=255)                                                                                  gcagcGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTAcgaccg                                                                      >  1:107485/1‑150 (MQ=255)                                                                                    agcGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTAcgaccgac                                                                    <  1:668727/150‑1 (MQ=255)                                                                                     gcgcGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACgaccgacc                                                                   >  2:1500670/1‑150 (MQ=255)                                                                                          tGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGt                                                              <  2:659144/150‑1 (MQ=255)                                                                                             gcgAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCgg                                                           <  1:363236/150‑1 (MQ=255)                                                                                              cgAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGcc                                                                                   <  1:265858/125‑1 (MQ=255)                                                                                              cgAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGcc                                                                                   >  2:265858/1‑125 (MQ=255)                                                                                                   gcaaccaCGACCCTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCCGATTG‑CCCAT‑‑CCTGATGCCTAATAATCTAGCCGTCCTCCATTGCTACTCGGCACGCGAAGGCTGGGCCAGCGGCTTCTACGACCGACCTCCGTCGGTTGCCa                                                     <  1:1407539/144‑1 (MQ=255)                                                                                                   gcAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCa                                                     <  1:1637548/150‑1 (MQ=255)                                                                                                    cAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTCGCCAc                                                    <  1:1041551/150‑1 (MQ=255)                                                                                                     aaGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTg                                                                                                                     >  1:1631314/1‑84 (MQ=255)                                                                                                     aaGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTg                                                                                                                     <  2:1631314/84‑1 (MQ=255)                                                                                                     aaGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCAcc                                                   <  1:867026/150‑1 (MQ=255)                                                                                                        tACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGc                                                >  1:216743/1‑150 (MQ=255)                                                                                                        tACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGc                                                <  2:436653/150‑1 (MQ=255)                                                                                                              cTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCAc                                          >  1:1783104/1‑150 (MQ=255)                                                                                                               ttGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGG‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCAcg                                         >  1:297256/1‑150 (MQ=255)                                                                                                                     tCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTgg                                   >  1:1203611/1‑150 (MQ=255)                                                                                                                       gACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCCTGCACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAg                                 >  1:1180758/1‑150 (MQ=255)                                                                                                                       gACGCACAGGGACTAATCCAGTCTTCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAg                                 <  2:1321015/150‑1 (MQ=255)                                                                                                                         cGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCTAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCg                               <  1:535199/150‑1 (MQ=255)                                                                                                                           cccAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCAACGGCCTCCACGGGGTGTAGCtgg                             <  2:1443045/150‑1 (MQ=255)                                                                                                                            gcAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGt                            >  1:378035/2‑150 (MQ=255)                                                                                                                            ccAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGCAg                                                                                               <  1:1734965/83‑1 (MQ=255)                                                                                                                            ccAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAg                                                                                               >  2:1734965/1‑83 (MQ=255)                                                                                                                                ggCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGc                        >  1:870512/1‑150 (MQ=255)                                                                                                                                 gCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGcc                       <  2:216743/150‑1 (MQ=255)                                                                                                                                         aggTCGGCCTGGTTGGCCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTCCGACCGGCCAGCACGCGGTGGAGCGGGTAAGCCAACTGttt               <  2:297256/148‑1 (MQ=255)                                                                                                                                         aaGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCTACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGttt               >  2:1431626/1‑150 (MQ=255)                                                                                                                                         aaGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGttt               <  1:664994/150‑1 (MQ=255)                                                                                                                                         aaGTCGCCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CGCTTAGCCGGCCTGGACCGCTACTCGGCGCGCGAATGCTGGGCCAGCCGCTTCTACTACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGttt               <  1:201424/150‑1 (MQ=255)                                                                                                                                          aGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTg              <  2:1716415/150‑1 (MQ=255)                                                                                                                                            tCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGttt               <  1:695726/147‑1 (MQ=255)                                                                                                                                            tCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGttt               >  2:695726/1‑147 (MQ=255)                                                                                                                                             cGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCa           <  2:870512/150‑1 (MQ=255)                                                                                                                                             cGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCa           <  2:107485/150‑1 (MQ=255)                                                                                                                                                ccTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATgg        <  1:1500670/150‑1 (MQ=255)                                                                                                                                                ccTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATgg        <  2:870308/150‑1 (MQ=255)                                                                                                                                                  tGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATGGcc      >  2:859136/1‑150 (MQ=255)                                                                                                                                                    gTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCGGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATGGCCCt    >  2:1756748/1‑150 (MQ=255)                                                                                                                                                      tGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCGGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATGGCCCTgg  <  1:1756748/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         CCAGGTCATCCCGGTGGCGATCCGTCGCTGGCGCAAGAACCTCTGGCAAACCATCAAGAGTGGCGAGTGGAAGGCGTTCAAGCAGCGCGTGCGCGAGCGCAAGTACGACCTTGTGATCGACGCCCAGGGCCTGGTCAAGTCGGCCTGGTTGACCCGTTACGTGAAGGC‑‑A‑CCCGTAGCCGGCCTGGACCGCTACTCGGCGCGCGAAGGCTGGGCCAGCCGCTTCTACGACCGACCTCTGTCGGTTGCCACCGGCCAGCACGCGGTGGAGCGGGTACGCCAACTGTTTGCCATGGCCCTGG  >  NC_002947/415391‑415689

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