Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	827,124	G→T	6.1%	K105N (AAG→AAT)	ycaC‑I →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	827,124	0	G	T	6.1%	212.4 / 3.6	65	K105N (AAG→AAT)	ycaC‑I	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base T (2/2);  total (32/33)
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.68e-01

CCACCAGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCA  >  NC_002947/826979‑827272                                                                                                                                                  |                                                                                                                                                     ccaccaGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCaagaag                                                                                                                                                   >  6:141947/1‑149 (MQ=255) ccaccaGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCaagaag                                                                                                                                                   >  6:107944/1‑149 (MQ=255)  caccaGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCaagaag                                                                                                                                                   >  1:262203/1‑148 (MQ=255)     caGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAagcagc                                                                                                                                               >  4:178543/1‑149 (MQ=255)            gagaCGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGAt                                                                                                                                        <  2:172800/149‑1 (MQ=255)                 ggggCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCAGCAATCAGCAGCTGCTTCTTGccg                                                                                                                                   >  5:206019/1‑147 (MQ=25)                 ggggCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCCACATTGCACGCCCTGGCAACATCAATGCCTGGGACAACGATGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTg                                                                                                                                   <  1:50337/149‑1 (MQ=255)                    gCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGgtg                                                                                                                                >  1:324542/1‑149 (MQ=255)                     cccAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCAGCAATCAGCAGCTGCTTCTTGCCGgtgg                                                                                                                               <  6:206019/149‑2 (MQ=255)                      ccAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGTCAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGtgtg                                                                                                                              <  7:231417/149‑1 (MQ=255)                        aaCGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGt                                                                                                                            >  7:294808/1‑149 (MQ=255)                          cGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCa                                                                                                                          <  3:178543/149‑1 (MQ=255)                             cccACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTg                                                                                                                       <  7:144445/149‑1 (MQ=255)                                aCTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAg                                                                                                                    >  3:20590/1‑149 (MQ=255)                                       ccGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGt                                                                                                             >  5:216307/1‑149 (MQ=255)                                        cGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTg                                                                                                            <  1:293836/149‑1 (MQ=255)                                         ggAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCAc                                                                                                                         >  8:279218/1‑135 (MQ=255)                                              tGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTc                                                                                                      >  3:147313/1‑149 (MQ=255)                                                     gCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGcc                                                                                                <  2:270354/148‑1 (MQ=255)                                                      cAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCt                                                                                              >  1:101217/1‑149 (MQ=255)                                                       aGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCAGCAATCAGCAGCTGCTTCTTGCCgg                                                                                                                                  <  7:173081/112‑1 (MQ=38)                                                       aGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCAGCAATCAGCAGCTGCTTCTTGCCgg                                                                                                                                  >  8:173081/1‑112 (MQ=38)                                                         ttCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTc                                                                                           <  5:141947/149‑1 (MQ=255)                                                                 tGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTg                                                                                    <  2:262203/148‑1 (MQ=255)                                                                   cccccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTggag                                                                                 <  5:107944/149‑1 (MQ=255)                                                                   cccccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTgga                                                                                  <  2:250531/148‑1 (MQ=255)                                                                     cccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTggag                                                                                 >  7:23669/1‑147 (MQ=255)                                                                     cccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTggag                                                                                 >  3:169443/1‑147 (MQ=255)                                                                      ccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGgaggag                                                                              >  4:140612/1‑149 (MQ=255)                                                                      ccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGgaggag                                                                              <  6:216307/149‑1 (MQ=255)                                                                       cTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGaggagg                                                                             <  8:294808/149‑1 (MQ=255)                                                                             tGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTc                                                                        >  5:44579/1‑148 (MQ=255)                                                                               cACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAg                                                                     <  8:234950/149‑1 (MQ=255)                                                                               cACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCGGGTGTGGTCACTGAAGTATGCGGGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAg                                                                     <  8:190371/149‑1 (MQ=255)                                                                                aCGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCtt                                                                         >  2:312613/1‑144 (MQ=255)                                                                                aCGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCtt                                                                         <  1:312613/144‑1 (MQ=255)                                                                                    ccGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTAt                                                                 >  5:175430/1‑148 (MQ=255)                                                                                        gCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTgga                                                                                  >  8:102181/1‑127 (MQ=255)                                                                                        gCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTgga                                                                                  <  7:102181/127‑1 (MQ=255)                                                                                          aaCATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTcgtcgt                                                          >  4:121905/1‑149 (MQ=255)                                                                                             aTCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCAc                                                       <  2:324542/149‑1 (MQ=255)                                                                                              tCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATtcgtc                                                            >  2:309838/1‑143 (MQ=255)                                                                                              tCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATtcgtc                                                            <  1:309838/143‑1 (MQ=255)                                                                                                  tGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACg                                                  >  2:300247/1‑149 (MQ=255)                                                                                                    ccTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGc                                                 <  2:101217/148‑1 (MQ=255)                                                                                                     cTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCATCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGcc                                                >  1:246895/1‑148 (MQ=255)                                                                                                        ggACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTg                                            <  1:300247/149‑1 (MQ=255)                                                                                                         gACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTgg                                           <  8:23669/149‑1 (MQ=255)                                                                                                            aaCGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCAc                                        >  7:231050/1‑149 (MQ=255)                                                                                                            aaCGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCAc                                        >  3:14801/1‑149 (MQ=255)                                                                                                               gAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGtt                                     >  4:48734/1‑149 (MQ=255)                                                                                                                   aTTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                 >  5:332604/1‑149 (MQ=255)                                                                                                                    tttCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                 >  6:132472/1‑148 (MQ=255)                                                                                                                      tCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGa                               <  8:231050/148‑1 (MQ=255)                                                                                                                           aaaGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGAc                         <  2:246895/149‑1 (MQ=255)                                                                                                                                  tAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGcc                                                                                                         >  3:211459/1‑62 (MQ=255)                                                                                                                                  tAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGcc                                                                                                         <  4:211459/62‑1 (MQ=255)                                                                                                                                  tAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTACAACGAACTGACCCGTgcc                  >  2:287187/1‑147 (MQ=255)                                                                                                                                    aaaGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAggg                                                                            <  1:172824/89‑1 (MQ=255)                                                                                                                                    aaaGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAggg                                                                            >  2:172824/1‑89 (MQ=255)                                                                                                                                    aaaGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGTCGCCTCTGGCACGTTCAACGAACTGACCCGTGACTc                <  1:287187/149‑1 (MQ=255)                                                                                                                                        ccACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGcc            >  4:152230/1‑149 (MQ=255)                                                                                                                                         cACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                 <  5:163881/127‑1 (MQ=255)                                                                                                                                         cACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                 >  6:163881/1‑127 (MQ=255)                                                                                                                                            cGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc        <  4:14801/149‑1 (MQ=255)                                                                                                                                            cGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc        <  6:332604/149‑1 (MQ=255)                                                                                                                                              gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                     <  5:111643/134‑1 (MQ=255)                                                                                                                                              gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                     >  6:111643/1‑134 (MQ=255)                                                                                                                                                 agaagCAGCAGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACg    >  5:74624/1‑148 (MQ=255)                                                                                                                                                  gaagCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCa  <  6:175430/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCACCAGTTTCGAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCA  >  NC_002947/826979‑827272

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