Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	827,136	G→T	5.8%	L109L (CTG→CTT)	ycaC‑I →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	827,136	0	G	T	5.8%	224.4 / 3.3	69	L109L (CTG→CTT)	ycaC‑I	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGG  >  NC_002947/826990‑827278                                                                                                                                                   |                                                                                                                                               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                                                                                    <  6:216307/149‑1 (MQ=255)                                                           ccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGgaggag                                                                                    >  4:140612/1‑149 (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                                                                               <  1:312613/144‑1 (MQ=255)                                                                     aCGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCtt                                                                               >  2:312613/1‑144 (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                                              >  3:14801/1‑149 (MQ=255)                                                                                                 aaCGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCAc                                              >  7:231050/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                                                                                                               <  4:211459/62‑1 (MQ=255)                                                                                                                       tAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGcc                                                                                                               >  3:211459/1‑62 (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                                       >  6:163881/1‑127 (MQ=255)                                                                                                                              cACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAAc                                       <  5:163881/127‑1 (MQ=255)                                                                                                                                 cGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc              <  6:332604/149‑1 (MQ=255)                                                                                                                                 cGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc              <  4:14801/149‑1 (MQ=255)                                                                                                                                   gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                           >  6:111643/1‑134 (MQ=255)                                                                                                                                   gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGt                           <  5:111643/134‑1 (MQ=255)                                                                                                                                      agaagCAGCAGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACg          >  5:74624/1‑148 (MQ=255)                                                                                                                                       gaagCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCa        <  6:175430/149‑1 (MQ=255)                                                                                                                                         agcagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATg      <  3:140612/149‑1 (MQ=255)                                                                                                                                            agcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAg   <  4:20590/149‑1 (MQ=255)                                                                                                                                             gctgctGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAgg  <  3:121905/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGG  >  NC_002947/826990‑827278

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