Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 627,348 G→A V427M (GTG→ATG)  entE → 2,3‑dihydroxybenzoate‑AMP ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913627,3480GA82.1% 126.2 / 23.8 56V427M (GTG→ATG) entE2,3‑dihydroxybenzoate‑AMP ligase
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base A (23/23);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

ATGATAGCGCGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCG  >  NC_000913/627103‑627586
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
atgatAGCGCGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAg                                                                                                                                                                                                                                           <  1:405141/251‑1 (MQ=255)
  gatAGAGCGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAggg                                                                                                                                                                                                                                         >  2:468038/1‑251 (MQ=255)
    tAGCGCGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGCAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGgc                                                                                                                                                                                                                                       >  1:119867/1‑251 (MQ=255)
     aGCGCGGAGAAAATTATCCATACCCAGTGTTGCCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGgc                                                                                                                                                                                                                                       <  2:314515/250‑1 (MQ=255)
      gcgcGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGTGCCGTACTCCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCACCGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGgcgc                                                                                                                                                                                                                                     <  1:113302/251‑1 (MQ=255)
             gaAAATTATCCATACCCAGGATTACCCAATGTGTCCGGATGACGAAGGATGGGGTGCCGATGCCGAAGGAAATCCACTGCGGTAAGGCGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCAGCGGATATTACACAAGTCCACAGCACAATGCCAGCGACTTTGCTGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTCTTGATCTAGAGGTCTACGTCACCATGCAGGGGCTCgagac                                                                                                                                                                                                                                >  2:393310/1‑248 (MQ=255)
                           cccAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGt                                                                                                                                                                                                                <  2:130174/251‑1 (MQ=255)
                            ccAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTg                                                                                                                                                                                                               >  1:97705/1‑251 (MQ=255)
                                           gtgtCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCgaga                                                                                                                                                                                                      >  2:93070/1‑245 (MQ=255)
                                             gtCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCt                                                                                                                                                                                              <  2:269294/251‑1 (MQ=255)
                                                      aCGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAagag                                                                                                                                                                                     <  1:106454/251‑1 (MQ=255)
                                                               gggTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAAc                                                                                                                                                                            >  2:482737/1‑251 (MQ=255)
                                                                    gccgatgccgaAGGAAATCCACTGCCGCAAGGGGAAGGCTGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATGGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACctgc                                                                                                                                                                        <  1:465506/250‑1 (MQ=255)
                                                                             gaAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGACTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGcc                                                                                                                                                               >  2:354033/1‑250 (MQ=255)
                                                                                       cACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTg                                                                                                                                                    <  2:4719/251‑1 (MQ=255)
                                                                                           gCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATc                                                                                                                                                 <  1:12471/250‑1 (MQ=255)
                                                                                                     gAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACt                                                                                                                                      <  2:399375/251‑1 (MQ=255)
                                                                                                      aaGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTg                                                                                                                                     <  2:380878/251‑1 (MQ=255)
                                                                                                             aCGCCTGATGACGCGCGGGCCGTACACCTTCCGCGTCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCCCCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCa                                                                                                                              <  1:250133/251‑1 (MQ=255)
                                                                                                               gCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGAGCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAGCCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATg                                                                                                                            <  2:205953/251‑1 (MQ=255)
                                                                                                                   gatgaCGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGaa                                                                                                                         <  1:448077/250‑1 (MQ=255)
                                                                                                                     tgaCGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGgcg                                                                                                                                                                                                                                      >  1:466871/1‑139 (MQ=255)
                                                                                                                     tgaCGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGgcg                                                                                                                                                                                                                                      <  2:466871/139‑1 (MQ=255)
                                                                                                                             gggCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGAt                                                                                                              <  1:468038/251‑1 (MQ=255)
                                                                                                                                 cGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGGGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGAGATCTACGCCGCACTGGTGAGCATGGCAGATGAGCTGATGGGc                                                                                                          >  2:83237/1‑251 (MQ=255)
                                                                                                                                  gTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCg                                                                                                         >  1:333018/1‑251 (MQ=255)
                                                                                                                                   tACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCg                                                                                                         >  1:259020/1‑250 (MQ=255)
                                                                                                                                   tACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGACAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCCCCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGCAGATGAGCTGATGGGCg                                                                                                         >  1:291950/1‑250 (MQ=255)
                                                                                                                                                          aaGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGCTGGGCGGAAAAAGCTGCGCTTATCtggtgg                                                                                 >  2:41673/1‑251 (MQ=255)
                                                                                                                                                                   gCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGc                                                                        <  2:202347/251‑1 (MQ=255)
                                                                                                                                                                          gCCAGCGCCTTTGATTCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCTCCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTgc                                                                  <  2:113018/250‑1 (MQ=255)
                                                                                                                                                                            cAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGCTGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTgcgcg                                                               >  1:417600/1‑251 (MQ=255)
                                                                                                                                                                             aGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTgcgcgc                                                              >  2:275040/1‑251 (MQ=255)
                                                                                                                                                                                   tttGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCCGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGCAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAaga                                                                          >  1:205696/1‑233 (MQ=255)
                                                                                                                                                                                   tttGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAaga                                                                          <  2:205696/233‑1 (MQ=255)
                                                                                                                                                                                     tGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCtggt                                                                                   >  2:50966/1‑222 (MQ=255)
                                                                                                                                                                                     tGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCtggt                                                                                   <  1:50966/222‑1 (MQ=255)
                                                                                                                                                                                         gCCAACGGTTTTTACTGGTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCTCTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGc                                                   <  2:20240/250‑1 (MQ=255)
                                                                                                                                                                                             aCGGTATTTACTGTTCCGGCGATATGATCTCTATTGATCCAGAGGGCTACATCACCGTGCAGGGGCGCGAGCAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGGGCTTATCTGGTGGTAAAAGAGCCACTGCGCGCGGTGCACGAGcgtcgt                                              >  2:61180/1‑251 (MQ=255)
                                                                                                                                                                                                  ttttACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCt                                         >  2:9406/1‑251 (MQ=255)
                                                                                                                                                                                                    ttACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTg                                        <  2:374524/250‑1 (MQ=255)
                                                                                                                                                                                                     tACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGc                                       <  1:93070/250‑1 (MQ=255)
                                                                                                                                                                                                        tGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTg                                    >  2:145782/1‑250 (MQ=255)
                                                                                                                                                                                                                 gatctgagctCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGATGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGGGCAGGTGCGTCGTTTCCTGCGTGAACAgggaa                           >  2:297358/1‑248 (MQ=255)
                                                                                                                                                                                                                      gatctCTATTGATCCAGAGGGTTACATCTCCATGCAGGGGCTCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCCCTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTa                           <  1:377364/245‑1 (MQ=255)
                                                                                                                                                                                                                      gatctCTATTGATCCAGAGGGTTACATCACCATGCATGTGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTa                           <  1:11504/245‑1 (MQ=255)
                                                                                                                                                                                                                       atctCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGc                                                   >  2:411479/1‑220 (MQ=255)
                                                                                                                                                                                                                       atctCTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGc                                                   <  1:411479/220‑1 (MQ=255)
                                                                                                                                                                                                                        tctcTATTGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTTTCTGGTGGGAAAAGAGCCGCTGCGCGCGGTGCCGGTGCGTCGTTTCCTGCGGGAACAGGGAATTGCCGac                   >  2:348893/1‑250 (MQ=255)
                                                                                                                                                                                                                              ttGATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTaa              >  1:26192/1‑250 (MQ=255)
                                                                                                                                                                                                                                gATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAAt            <  2:429899/250‑1 (MQ=255)
                                                                                                                                                                                                                                gATCCAGAGGGTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGAGATCTACGCCGCACTGGAGAGAATGGAAGAACCGCTGATGGGAGAACAAAGCTTCGCTTATCTGGGGGTAAAAAAGCCGCTGCGAGCGGTGCAGGAGCCTCGTTTACTGCGTGAACAGGGTATTGCCGCATTTAAAtt           >  1:300487/1‑251 (MQ=255)
                                                                                                                                                                                                                                     agagGGTTCCATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCgg      <  2:381905/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gagGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCgg      >  1:366009/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gggTTACATCACCATGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATc   <  1:354033/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTTACATCACCATGCAGGGGCGCGAGAAAGATAAGATTAACCGTGGCGGAGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGCTCTACGCCGCACTGGTGAGCATGGACGATGAGCTGCTGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGCATTTAAATTACCGGATcg  >  1:443669/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
ATGATAGCGCGGAGAAAATTATCCATACCCAGGGTTACCCAATGTGTCCGGATGACGAAGTATGGGTTGCCGATGCCGAAGGAAATCCACTGCCGCAAGGGGAAGTCGGACGCCTGATGACGCGCGGGCCGTACACCTTCCGCGGCTATTACAAAAGTCCACAGCACAATGCCAGCGCCTTTGATGCCAACGGTTTTTACTGTTCCGGCGATCTGATCTCTATTGATCCAGAGGGTTACATCACCGTGCAGGGGCGCGAGAAAGATCAGATTAACCGTGGCGGCGAGAAGATCGCTGCCGAAGAGATCGAAAACCTGCTGCTGCGCCACCCGGCGGTGATCTACGCCGCACTGGTGAGCATGGAAGATGAGCTGATGGGCGAAAAAAGCTGCGCTTATCTGGTGGTAAAAGAGCCGCTGCGCGCGGTGCAGGTGCGTCGTTTCCTGCGTGAACAGGGTATTGCCGAATTTAAATTACCGGATCG  >  NC_000913/627103‑627586

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