Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,285,238 G→A intergenic (‑25/+175) yjcH ← / ← acs DUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,285,2380GA95.5% 201.4 / ‑0.6 66intergenic (‑25/+175)yjcH/acsDUF485 domain‑containing inner membrane protein YjcH/acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (31/32);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCG  >  NC_000913/4284992‑4285483
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gaTCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAt                                                                                                                                                                                                                                                   <  1:387902/251‑1 (MQ=255)
       aTCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGtttt                                                                                                                                                                                                                                             >  1:426252/1‑250 (MQ=255)
                    cAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCt                                                                                                                                                                                                                               <  1:394803/251‑1 (MQ=255)
                            aTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAAATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATggg                                                                                                                                                                                                                         >  1:299420/1‑249 (MQ=255)
                                         cGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTc                                                                                                                                                                                                          <  1:81609/251‑1 (MQ=255)
                                              gtgtTCGGATTCAGCCGCGTGCCCAGCCAGCTGGGCGGGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAAACGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCACTGAACTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTa                                                                                                                                                                                                     <  1:491318/251‑1 (MQ=255)
                                              gtgtGCGGTTTCAGCGGCGTGCCCAGCCAGCCGGGCGCTAAGGCGATCAGTAAAATAAAGCCGATATAACCTGCCAGCATAATAATCGACAGGATGGTGGCAATCCTTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTTTCTTCTATCAGCTGATAAATAGTTCCATTCATCACAGATTCGCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTa                                                                                                                                                                                                     <  2:153317/251‑1 (MQ=255)
                                                           gcggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAAt                                                                                                                                                                                        <  1:24780/251‑1 (MQ=255)
                                                               cgggCCCAGCCAGCTGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAAGCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCAGCACAGATTCGCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCCTGCCGAATGTAg                                                                                                                                                                                    <  2:362565/248‑1 (MQ=255)
                                                                  gccctgcccgccagtcgcgAATGCGATAAGTAAAATAAAACCAAGATAAACGGCCAGCATAAGAATCTACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAAATCCTTGAAATGCGCATTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCTCAGATTCACCAGAGGTAATGTAAGGAGTGTTTTAATTCCCGCTCCCGGATGGGTGAAGGGTAACGCTCGGGTAACCCTTGCCGACTGTCGGc                                                                                                                                                                                  <  2:192769/235‑1 (MQ=255)
                                                                    ccagccagCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCgg                                                                                                                                                                               <  2:22010/251‑1 (MQ=255)
                                                                      agccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGAt                                                                                                                                                                             >  2:213624/1‑251 (MQ=255)
                                                                          agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAgg                                                                                                                                                                         <  1:213624/251‑1 (MQ=255)
                                                                           gccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGc                                                                                                                                                                        <  1:87824/251‑1 (MQ=255)
                                                                              gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                                                                                                                                                     <  1:482667/251‑1 (MQ=255)
                                                                                gcgcGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGttt                                                                                                                                                                    >  1:444445/1‑250 (MQ=255)
                                                                                              aGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                                      >  1:137935/1‑250 (MQ=255)
                                                                                                    aTAAAGCCGATATAAACTTCCAGCATAATAATCGACCGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGGTTCGCCAGAGGTAATGTAAGGATTGTTTTATTTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACACTTTCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcactcaat                                                                                                                                               >  2:308086/1‑251 (MQ=255)
                                                                                                        aGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                                           <  1:397924/251‑1 (MQ=255)
                                                                                                        aGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                                           >  2:82397/1‑251 (MQ=255)
                                                                                                          ccGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGTGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTg                                                                                                                                          >  1:118750/1‑250 (MQ=255)
                                                                                                            gaaataAACTGCCAGCATAATAATCGACAGGAAGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTGTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                                       <  2:205111/248‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAGATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 <  2:215633/240‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                                                      <  1:75744/251‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 >  1:215633/1‑240 (MQ=255)
                                                                                                                 aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACg                                                                                                                                  >  1:124335/1‑251 (MQ=255)
                                                                                                                   aCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                                                                <  1:131421/251‑1 (MQ=255)
                                                                                                                            ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                                       >  1:210498/1‑251 (MQ=255)
                                                                                                                            ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                                       >  1:286713/1‑251 (MQ=255)
                                                                                                                               ataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                                    >  1:60404/1‑251 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGTAATGCGCATTGTCTGCTATCCGCTGAGAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTACCGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 <  2:269120/221‑1 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 >  1:269120/1‑221 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTaac                                                                                                                   >  2:365292/1‑249 (MQ=255)
                                                                                                                                 aatCGACAGGATGGTGGCAAACCTTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                  <  1:95833/251‑1 (MQ=255)
                                                                                                                                  atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                                                 <  2:17102/251‑1 (MQ=255)
                                                                                                                                        agggtggttgCAAACCGTTGCCGTTTTTCTACTAACTCCCGGAAAAGCGCATGTTCTTCTATACGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCCATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                           <  2:124251/247‑1 (MQ=38)
                                                                                                                                        aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                           <  1:315998/251‑1 (MQ=255)
                                                                                                                                         ggATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCAGGcccac                                                                                                           >  2:268112/1‑247 (MQ=255)
                                                                                                                                           aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTCCaa                                                                                                         >  2:181615/1‑250 (MQ=255)
                                                                                                                                              gtggCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCGTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAGTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCg                                                                                                     <  2:411932/251‑1 (MQ=255)
                                                                                                                                                 gCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGtt                                                                                                   >  1:495601/1‑250 (MQ=255)
                                                                                                                                                  cAAACCGTTGCCGTTTTTCGACTAACTCACTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                                 >  1:441369/1‑251 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                          >  2:1397/1‑193 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                          <  1:1397/193‑1 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                                                >  2:309281/1‑251 (MQ=255)
                                                                                                                                                     aCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGCGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGAGGAGAAGCTGTCGCGTCTTATCAGGCCTACAAACCGGTACCg                                                                                               >  1:73585/1‑250 (MQ=255)
                                                                                                                                                       cGTTGCCGTTTTTCGACGAACTCCCTGAAATGCGCATTGGCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACt                                                                                            <  2:206880/251‑1 (MQ=255)
                                                                                                                                                               gttttTGACTAACTCCCTGAATTGCGCATTGGCTCCGAACCGCTGATAAATAGTGCTATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAGGGGAGAAGGTTAACGCGCGGGTAACCCTTGCCGAATGTATGCCGGATAAGGCGTTTACGCCGCAGCCGGCAATCAATTCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCg                                                                                     <  2:247946/249‑1 (MQ=255)
                                                                                                                                                                tttgcGACTAATTCCCTGAAATGCGCATTGCTCTCTACCCGCTGATAAATCGTGCCATTCAGTACAGATTCGCAAGAGTTAATGGAAGGATTGTTCTAATTCCCGCTCCCTTAGGGGAGAGGGTTAACGCTCGGGTATCCCGTGCCGAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCTTCTTATCAGGCCTAAAAACCGTTACCGACCCGCATCgg                                                                                    <  2:470177/246‑1 (MQ=25)
                                                                                                                                                                 tttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGACGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGc                                                                                  >  2:47172/1‑251 (MQ=255)
                                                                                                                                                                         aaCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCGTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTgg                                                                          <  1:226741/251‑1 (MQ=255)
                                                                                                                                                                         aaCTCCCTGAAATGCGCATTGTCTACTATCCGCTGAGAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTTCCGAATGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCGGGCAATTGTgg                                                                          <  2:272481/251‑1 (MQ=255)
                                                                                                                                                                              ccTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTAc                                                                     >  1:362039/1‑251 (MQ=255)
                                                                                                                                                                                 gAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCGGGCACTTGTGGGTTACGAt                                                                  >  2:503324/1‑251 (MQ=255)
                                                                                                                                                                                        gcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATcg                                                           <  2:495601/251‑1 (MQ=255)
                                                                                                                                                                                          aTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATcgcg                                                         >  1:308931/1‑251 (MQ=255)
                                                                                                                                                                                             ggcttctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCGATTGTGGGTTACGATGGCATCGCGAt                                                       <  2:138889/248‑1 (MQ=255)
                                                                                                                                                                                                           tggtAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCCTCGGGCAATTGTGGGTTACGCTGGCATCGCGATACCCTGCTTCTCtt                                         >  2:171591/4‑251 (MQ=255)
                                                                                                                                                                                                                  tAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAagc                                    >  2:192421/1‑249 (MQ=255)
                                                                                                                                                                                                                   aGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                                                                    <  1:200043/152‑1 (MQ=255)
                                                                                                                                                                                                                   aGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                                                                    >  2:200043/1‑152 (MQ=255)
                                                                                                                                                                                                                               cacaGAGTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTAc                      <  2:153781/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCATTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGACTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGACTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGcc                   >  2:193400/1‑249 (MQ=255)
                                                                                                                                                                                                                                           agagGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTCCGCCAGGATCGGcc         >  2:404743/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 aaTGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTc   <  1:16747/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCg  >  1:153404/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCG  >  NC_000913/4284992‑4285483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: