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

GATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCG  >  NC_000913/4284992‑4285483
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gaTCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAt                                                                                                                                                                                                                                                   <  1:387902/251‑1 (MQ=255)
       aTCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGtttt                                                                                                                                                                                                                                             >  1:426252/1‑250 (MQ=255)
                    cAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCt                                                                                                                                                                                                                               <  1:394803/251‑1 (MQ=255)
                            aTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAAATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATggg                                                                                                                                                                                                                         >  1:299420/1‑249 (MQ=255)
                                         cGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTc                                                                                                                                                                                                          <  1:81609/251‑1 (MQ=255)
                                              gtgtTCGGATTCAGCCGCGTGCCCAGCCAGCTGGGCGGGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAAACGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCACTGAACTGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTa                                                                                                                                                                                                     <  1:491318/251‑1 (MQ=255)
                                              gtgtGCGGTTTCAGCGGCGTGCCCAGCCAGCCGGGCGCTAAGGCGATCAGTAAAATAAAGCCGATATAACCTGCCAGCATAATAATCGACAGGATGGTGGCAATCCTTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTTTCT‑TCTATCAGCTGATAAATAGTTCCATTCATCACAGATTCGCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTa                                                                                                                                                                                                     <  2:153317/251‑1 (MQ=255)
                                                           gcggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAAt                                                                                                                                                                                        <  1:24780/251‑1 (MQ=255)
                                                               cgggCCCAGCCAGCTGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAAGCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCAGCACAGATTCGCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCCTGCCGAATGTAg                                                                                                                                                                                    <  2:362565/248‑1 (MQ=255)
                                                                  gccctgcccgccagtcgcgAATGCGATAAGTAAAATAAAACCAAGATAAACGGCCAGCATAAGAATCTACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAAATCCTTGAAATGCGCATTGTCT‑TCTATCCGCTTATAAATAGTGCCATTCATCTCAGATTCACCAGAGGTAATGTAAGGAGTGTTTTAATTCCCGCTCCCGGATGGGTGAAGGGTAACGCTCGGGTAACCCTTGCCGACTGTCGGc                                                                                                                                                                                  <  2:192769/235‑1 (MQ=255)
                                                                    ccagccagCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCgg                                                                                                                                                                               <  2:22010/251‑1 (MQ=255)
                                                                      agccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGAt                                                                                                                                                                             >  2:213624/1‑251 (MQ=255)
                                                                          agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAgg                                                                                                                                                                         <  1:213624/251‑1 (MQ=255)
                                                                           gccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGc                                                                                                                                                                        <  1:87824/251‑1 (MQ=255)
                                                                              gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                                                                                                                                                     <  1:482667/251‑1 (MQ=255)
                                                                                gcgcGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGttt                                                                                                                                                                    >  1:444445/1‑250 (MQ=255)
                                                                                              aGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                                      >  1:137935/1‑250 (MQ=255)
                                                                                                    aTAAAGCCGATATAAACTTCCAGCATAATAATCGACCGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGGTTCGCCAGAGGTAATGTAAGGATTGTTTTATTTCCCGCTCCCTTATGGGAGAAGGGTAACGCTCGGGTAACACTTTCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcactcaat                                                                                                                                               >  2:308086/1‑251 (MQ=255)
                                                                                                        aGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                                           <  1:397924/251‑1 (MQ=255)
                                                                                                        aGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                                                                                                                                           >  2:82397/1‑251 (MQ=255)
                                                                                                          ccGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGTGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTg                                                                                                                                          >  1:118750/1‑250 (MQ=255)
                                                                                                            gaaataAACTGCCAGCATAATAATCGACAGGAAGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTGTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATg                                                                                                                                       <  2:205111/248‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAGATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 <  2:215633/240‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATCAATGCCTGATGc                                                                                                                                      <  1:75744/251‑1 (MQ=255)
                                                                                                             atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 >  1:215633/1‑240 (MQ=255)
                                                                                                                 aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACg                                                                                                                                  >  1:124335/1‑251 (MQ=255)
                                                                                                                   aCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCt                                                                                                                                <  1:131421/251‑1 (MQ=255)
                                                                                                                            ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                                       >  1:210498/1‑251 (MQ=255)
                                                                                                                            ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                                                                       >  1:286713/1‑251 (MQ=255)
                                                                                                                               ataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                                    >  1:60404/1‑251 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGTAATGCGCATTGTCT‑GCTATCCGCTGAGAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTACCGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 <  2:269120/221‑1 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatca                                                                                                                                                 >  1:269120/1‑221 (MQ=255)
                                                                                                                                taatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTaac                                                                                                                   >  2:365292/1‑249 (MQ=255)
                                                                                                                                 aatCGACAGGATGGTGGCAAACCTTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                  <  1:95833/251‑1 (MQ=255)
                                                                                                                                  atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAg                                                                                                                 <  2:17102/251‑1 (MQ=255)
                                                                                                                                        agggtggttgCAAACCGTTGCCGTTTTTCTACTAACTCCCGGAAAAGCGCATGTTCT‑TCTATACGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCCATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                           <  2:124251/247‑1 (MQ=38)
                                                                                                                                        aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTAc                                                                                                           <  1:315998/251‑1 (MQ=255)
                                                                                                                                         ggATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGACGCGTCTTATCAGGcccac                                                                                                           >  2:268112/1‑247 (MQ=255)
                                                                                                                                           aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTCCaa                                                                                                         >  2:181615/1‑250 (MQ=255)
                                                                                                                                              gtggCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCG‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGAGTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCg                                                                                                     <  2:411932/251‑1 (MQ=255)
                                                                                                                                                 gCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGtt                                                                                                   >  1:495601/1‑250 (MQ=255)
                                                                                                                                                  cAAACCGTTGCCGTTTTTCGACTAACTCACTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAc                                                                                                 >  1:441369/1‑251 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                          >  2:1397/1‑193 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                          <  1:1397/193‑1 (MQ=255)
                                                                                                                                                   aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTAcc                                                                                                >  2:309281/1‑251 (MQ=255)
                                                                                                                                                     aCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGCGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGAGGAGAAGCTGTCGCGTCTTATCAGGCCTACAAACCGGTACCg                                                                                               >  1:73585/1‑250 (MQ=255)
                                                                                                                                                       cGTTGCCGTTTTTCGACGAACTCCCTGAAATGCGCATTGGCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACt                                                                                            <  2:206880/251‑1 (MQ=255)
                                                                                                                                                               gttttTGACTAACTCCCTGAATTGCGCATTGGCT‑CCGAACCGCTGATAAATAGTGCTATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAGGGGAGAAGGTTAACGCGCGGGTAACCCTTGCCGAATGTATGCCGGATAAGGCGTTTACGCCGCAGCCGGCAATCAATTCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCg                                                                                     <  2:247946/249‑1 (MQ=255)
                                                                                                                                                                tttgcGACTAATTCCCTGAAATGCGCATTG‑CTCTCTACCCGCTGATAAATCGTGCCATTCAGTACAGATTCGCAAGAGTTAATGGAAGGATTGTTCTAATTCCCGCTCCCTTAGGGGAGAGGGTTAACGCTCGGGTATCCCGTGCCGAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCTTCTTATCAGGCCTAAAAACCGTTACCGACCCGCATCgg                                                                                    <  2:470177/246‑1 (MQ=25)
                                                                                                                                                                 tttCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGACGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGc                                                                                  >  2:47172/1‑251 (MQ=255)
                                                                                                                                                                         aaCTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCGTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTgg                                                                          <  1:226741/251‑1 (MQ=255)
                                                                                                                                                                         aaCTCCCTGAAATGCGCATTGTCT‑ACTATCCGCTGAGAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTTCCGAATGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCGGGCAATTGTgg                                                                          <  2:272481/251‑1 (MQ=255)
                                                                                                                                                                              ccTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTAc                                                                     >  1:362039/1‑251 (MQ=255)
                                                                                                                                                                                 gAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGCCTCGCATCGGGCACTTGTGGGTTACGAt                                                                  >  2:503324/1‑251 (MQ=255)
                                                                                                                                                                                        gcATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATcg                                                           <  2:495601/251‑1 (MQ=255)
                                                                                                                                                                                          aTTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATcgcg                                                         >  1:308931/1‑251 (MQ=255)
                                                                                                                                                                                             ggct‑tctATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAAGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCGATTGTGGGTTACGATGGCATCGCGAt                                                       <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCT‑TCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCG  >  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: ‑