Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	827,137	A→C	5.7%	I110L (ATT→CTT)	ycaC‑I →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	827,137	0	A	C	5.7%	216.7 / 2.2	69	I110L (ATT→CTT)	ycaC‑I	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base C (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 = 9.72e-01

GAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCC  >  NC_002947/826990‑827283                                                                                                                                                    |                                                                                                                                                   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                                                                                                                                             >  8:173081/1‑112 (MQ=38)                                            aGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCAGCAATCAGCAGCTGCTTCTTGCCgg                                                                                                                                             <  7:173081/112‑1 (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                                                                                            >  3:169443/1‑147 (MQ=255)                                                          cccTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTggag                                                                                            >  7:23669/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                                                                                    >  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                                                   >  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)                                                                                                                                                    aTTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGcc  <  4:147313/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAGACGGGGCCCAACGGCCCACTGGTGCCGGAACTGAAAGAGCAGTTCCCTGATGCCCCCTACATTGCACGCCCGGGCAACATCAATGCCTGGGACAACGAAGATTTCGTCAAAGCGGTAAAAGCCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCC  >  NC_002947/826990‑827283

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