Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,575 T→C V181V (GTA→GTG)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,5750TC95.5% 195.3 / ‑4.0 66V181V (GTA→GTG) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (41/22);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 5.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGA  >  NC_000913/4623338‑4623819
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
aaaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AAt                                                                                                                                                                                                                                            <  1:246815/249‑1 (MQ=255)
aaaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGc                                                                                                                                                                                                                                          <  1:15805/251‑1 (MQ=255)
     gtacctCCAGCAGATGCTAGAAGGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTc                                                                                                                                                                                                                                     <  2:565100/246‑1 (MQ=255)
        cATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTcgcc                                                                                                                                                                                                                                  <  2:240777/251‑1 (MQ=255)
           ccAGCAGATGCGAGAATGCCGGAGCCTGACCGCAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTATCCGGGGAGATGATTTCCGATTATACACGCTATACATAATGGGCGCAAAAGGCCTTGGTTATGGCCTCGAATACCAGCTCATGGGTGATCCCCGGCAGTAGCTCGGTGAGACTGGTAACGCGGGAACGCACTG‑CATGAATGCCTTTCgccgcc                                                                                                                                                                                                                               >  2:28774/1‑251 (MQ=255)
             agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTTGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCGCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAg                                                                                                                                                                                                                             <  2:371130/251‑1 (MQ=255)
                   ggCGAGAATGCCGGGGCCTGACCGAAGTTCCTTTCCCAGCTTCTCTGACGGGCAAAGGTTTCGGCGAAGTTGGGCAAGGCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACCCGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGAGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCtt                                                                                                                                                                                                                       <  2:66009/250‑1 (MQ=255)
                       agaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGCTTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCtttttt                                                                                                                                                                                                                   >  2:486906/1‑250 (MQ=255)
                            gCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCg                                                                                                                                                                                                               <  1:615894/250‑1 (MQ=255)
                              cGGAGCCTGACCGAAGGTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAt                                                                                                                                                                                                            <  1:87375/251‑1 (MQ=255)
                               ggAGCCTGACCGACGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGACGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCACGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGA‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                                                                                                                                                           >  1:622251/1‑251 (MQ=255)
                               ggAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                                                                                                                                                           >  1:625020/1‑251 (MQ=255)
                                   ccTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGCTGATTTACGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGATTATGGCCTAGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTAGCCGCCAGTTTCTTTTTATcagggatgagg                                                                                                                                                                                                       >  2:640427/1‑241 (MQ=255)
                                      gACCGAAGTACCATTCCCAGCTACTGTGAGGGGCAAAGGTTTCGGGTAAGTTGGGAAAGTCTGGCGTTTTGT‑CCGGGAGATGATTTCCGCTTCTACACTCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTAGCCGGATTGAGATAg                                                                                                                                                                                                    <  1:496789/250‑1 (MQ=255)
                                           aaGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTGCGGCGAAGTTGGGCAAGTCTGGCGTTTTGTCCGGGTAGATGATTTCCGCTGCGACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGGGAGGTTTGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTTAGATAGTTTGc                                                                                                                                                                                               <  2:622251/251‑1 (MQ=255)
                                                                 gacgtgCAAAGGTTTCGGCGAAGTTTGGCAAGTTGGGTGTGTTGGCCGGGGAGATGATTGCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCCCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCTGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCAtt                                                                                                                                                                          <  2:625020/245‑1 (MQ=255)
                                                                            gTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATACTGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAAGGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGAGAGTTTGCCAAGCGGCTAAGGGCGGCCTTGACC‑AGCAAggg                                                                                                                                                             >  2:322347/1‑250 (MQ=255)
                                                                                 ggCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgt                                                                                                                                                        >  1:584754/1‑251 (MQ=255)
                                                                                 ggCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCCGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgt                                                                                                                                                        >  1:83547/1‑251 (MQ=255)
                                                                                 ggCGAAGTTTGGCAAGACTGGCGTTTTGTCCGGGGAGATGCTTTCAGCTTCTACACGCTCGCCATACTGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTAGGTGAGGTTGGTACCGCGGGAAAGCACTGAGGT‑AATGCCTTTCGCCGCCAGGTTCTTTTTATCCGGATTGAGATAGTTGGCCCGGCGGCTAAGGTCGGCATT‑TGTCAGCAAGGTGCCgg                                                                                                                                                        >  2:538409/1‑250 (MQ=255)
                                                                                    gAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAAATGATTTCCGCTTCTACACTCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCGGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgtggt                                                                                                                                                     <  2:641912/251‑1 (MQ=255)
                                                                                           ggCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTg                                                                                                                                                                         >  1:369649/1‑225 (MQ=255)
                                                                                            gCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGcca                                                                                                                                             <  1:619667/251‑2 (MQ=255)
                                                                                             cAAGTCTGTGGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCCCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGccac                                                                                                                                            <  2:545451/251‑3 (MQ=255)
                                                                                             cAAGTCTGGCGTTTTGCCCGGGGCGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCTGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGctcc                                                                                                                                            >  2:382994/1‑248 (MQ=255)
                                                                                                   tGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCtt                                                                                                                                      >  1:363176/1‑251 (MQ=255)
                                                                                                   gggCGTTTGGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTGGGTCACGCGGGAACGCACTGATGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGTCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCtt                                                                                                                                      <  1:157567/250‑1 (MQ=255)
                                                                                                     gCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTTGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTg                                                                                                                                    >  2:620229/1‑251 (MQ=255)
                                                                                                      cGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTgg                                                                                                                                   >  1:390881/1‑251 (MQ=255)
                                                                                                           ggtcgggggggATGATTTCCGCTTCTACACTGTCGACATAATGGGCGAAAAAGGCCTTGGTTAGGGCCTCGCAATCCTGCTCATGGGTGATCCCCGGCAGCAGCGCGGTGAGGTTGGTCACGCGGGAACGCTCGGAGGT‑AATTCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTGGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTc                                                                                                                               <  2:584754/245‑1 (MQ=255)
                                                                                                            gtccGGGGAGATGATTTCAGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGCTTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGCGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTcgc                                                                                                                             >  2:109943/4‑251 (MQ=255)
                                                                                                                 gggggttgatTTCCGCGTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTGATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGGGGTGGAAGCCACGATCTTTGGGTTCTCGATAg                                                                                                                        <  2:286843/245‑1 (MQ=255)
                                                                                                                     gatgatTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCg                                                                                                                    >  1:666201/1‑251 (MQ=255)
                                                                                                                         atTTCAGCTTCTACACGCGCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACAGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGGTTGAGATAGTTTGCCAGGCGGCTACTGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGACGCCACGATCTTTGCGTTCGCGATAGGCCGaacc                                                                                                                >  1:73672/1‑248 (MQ=255)
                                                                                                                                cttctacacgcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGACGAGGTGGAAGCCACGATTTTTGGTTTCGCGATAGGCCGAACCTGAGACt                                                                                                         >  2:117934/9‑251 (MQ=255)
                                                                                                                                  tctacacgcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGGGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGGCttt                                                                                                       >  2:339185/7‑251 (MQ=255)
                                                                                                                                             gCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AGTGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCtt                                                                                            >  1:295172/1‑251 (MQ=255)
                                                                                                                                             gCCATAACGGGCGAAAAAGGCCTTGGTTATGGAATAGCAAACCTGCTCATGGGGGATCACCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGACCGCACTGAGGT‑AATGCCTTTCGCGGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCCAGGGGCCGAGGTGGATGCTACGTTCTTTGGCTTCGCGCTAGGCCGAGCCTGAGACTTTGCGGtggactt                                                                                            >  2:565540/1‑245 (MQ=255)
                                                                                                                                              ccATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTc                                                                                           <  1:91461/251‑1 (MQ=255)
                                                                                                                                              ccATAATAGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCCGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCACGGTGCCGTGGTGGAAGCCCCGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTTCGGTCGCCTTc                                                                                           >  2:71047/1‑251 (MQ=255)
                                                                                                                                                     gggCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGtt                                                                                    >  1:135074/1‑251 (MQ=255)
                                                                                                                                                     gggCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGtt                                                                                    >  1:403983/1‑251 (MQ=255)
                                                                                                                                                      tgCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGCTTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCACGGCGCCGTGGTGGACGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGAATTGAGGGTCGCCTTCggagggtt                                                                                   >  2:547248/2‑244 (MQ=255)
                                                                                                                                                           aaaaGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTcacc                                                                              <  1:21174/251‑1 (MQ=255)
                                                                                                                                                                 ccttggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGGGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCCGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCCCGATCTTTGGTTTCGCGATAGGCCGCTCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTcaccccc                                                                           >  1:268086/5‑245 (MQ=255)
                                                                                                                                                                    tggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGAt                                                                       <  1:224560/248‑1 (MQ=255)
                                                                                                                                                                     ggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGgtttacacaccaaaccct                                                                    >  2:661649/1‑237 (MQ=255)
                                                                                                                                                                        tATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGTACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTa                                                                  <  1:356535/250‑1 (MQ=255)
                                                                                                                                                                        tATGGCATCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCAACAGATCGTTAc                                                                 >  1:113962/1‑251 (MQ=255)
                                                                                                                                                                         aTGGCCTCGCAAACCTGCTCATGGGTGACCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTAc                                                                 <  1:491125/250‑1 (MQ=255)
                                                                                                                                                                           ggCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGGTTTCACCACCAGATCGTTACgcc                                                              >  1:652859/1‑249 (MQ=255)
                                                                                                                                                                             ccTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTaagcccg                                                            >  1:345520/1‑245 (MQ=255)
                                                                                                                                                                              cTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGCGGTTGGTCACGCGGGACCGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCCTTGAGC‑AGCAAGGTGCCGTGGGGGAAGCAACGATCTTTGGTTTCGCGCTCGGGCGCTCCTGAGCCTTTGCGGTCGCCTTCGGCGTGTTTCACCACCGGATCGtatccccccg                                                           >  2:657649/1‑242 (MQ=255)
                                                                                                                                                                                 gCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGCCTTTGCGGTCGCCTTCGGCGGTTTTTACCCCCCGAGCGGTACGCCCGGACg                                                        >  2:170548/1‑251 (MQ=255)
                                                                                                                                                                                         gcccatggggtaTCCCCGGCGTCAGCTCGGTGAGGTTGGTCACGCGGGAACGTACGGAGGT‑AATGCCTTTCGCCGCCAGTTTCGTTTTATCCCGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCAGTGAGC‑AGCAAGGGGCCGTTGTGGAAGCCACGATCTTTGGTTTCGCGATATGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCgg                                                  <  2:577174/238‑1 (MQ=255)
                                                                                                                                                                                          cTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCAACAGATCGGTACGCCCGGACGCTTCGgcgc                                               >  1:506825/1‑251 (MQ=255)
                                                                                                                                                                                            cATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTg                                             >  1:498080/1‑251 (MQ=255)
                                                                                                                                                                                              tGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGatcc                                                                                                                >  1:13621/1‑179 (MQ=255)
                                                                                                                                                                                                  tGATCCCCGGCAGCAGCTCGGTGAGATTGGTCACGAGGGACCGCACTGAGGT‑AATGCCTTGCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGTCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGGAGCCACGATCTTTGGTTTCTCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTACCCCCCAACTCGTTCCGCCCGGACGCTTCTGCTCTCACGccc                                       >  2:336000/1‑250 (MQ=255)
                                                                                                                                                                                                       cccGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGCTCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGGTCGTTACGCCCCGGCGCTTCGgcggttgcgccaagccc                                  >  1:248254/1‑237 (MQ=255)
                                                                                                                                                                                                       cccGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAgcgc                                  >  1:469651/1‑251 (MQ=255)
                                                                                                                                                                                                           gcagcAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTc                              <  1:486906/247‑1 (MQ=255)
                                                                                                                                                                                                                aGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCCCCAGATCGTTACGCCCGGACGCTTCCGCGCCGACGCCAAGCGCGTTCAgccc                         >  1:644171/1‑249 (MQ=255)
                                                                                                                                                                                                                                 cACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGAc         <  1:135981/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcgGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGAAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGATTTCACCACCAAATCGGTACGCACGGACGATTAGGCCCTGACGCCAAGCGCGCTACGCGCATTGAGCACAAtagccgtg      >  1:124203/1‑244 (MQ=255)
                                                                                                                                                                                                                                       ggAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGGACGTTACACCCGGCCGCTTCGGCGCTGCCGCCACGCGCGTTCAGCGCagtgcgctcaacctacgttgagg  >  2:224981/1‑229 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
AAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGA  >  NC_000913/4623338‑4623819

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑