Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,336 G→C T261S (ACC→AGC)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,3360GC98.6% 226.7 / ‑3.8 70T261S (ACC→AGC) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (34/35);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

CGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACC  >  NC_000913/4623094‑4623578
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
cagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGt                                                                                                                                                                                                                                            >  1:283879/9‑251 (MQ=255)
 agcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGtt                                                                                                                                                                                                                                           >  1:575533/8‑251 (MQ=255)
   cgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTCTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTc                                                                                                                                                                                                                                          >  1:565100/6‑250 (MQ=255)
     agttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATc                                                                                                                                                                                                                                       >  1:66009/4‑251 (MQ=255)
       ttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCCGCTAAAGCGTTCATccc                                                                                                                                                                                                                                     >  2:286819/2‑250 (MQ=255)
           cTTACAGCCCACGCTGTCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCCACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGAGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGGGGTATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCCCACGCGTTCATCcaaca                                                                                                                                                                                                                                  >  2:546301/1‑247 (MQ=255)
             taccgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCcagcag                                                                                                                                                                                                                                 <  1:405603/245‑1 (MQ=255)
              accgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGAt                                                                                                                                                                                                                               <  1:199255/247‑1 (MQ=255)
              accgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGAGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATg                                                                                                                                                                                                                              >  1:656407/4‑251 (MQ=255)
               tcgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGATCCGAGAAAGCAACCAACAGCGCTTCGCACTCTTGTTGCAGAATATCCGCGCGGTCCAGGCAGCCTTGCAGTCGCCCGGCGCGGTCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAACCACCTGGGCGCGGGCGATCTGGCCTTTTTCAACGTCGCACTACCGTTCCACGCCGCCCCAGCTAAAGCTATCATCcagccgttgg                                                                                                                                                                                                                             >  2:496789/3‑245 (MQ=255)
                              acgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGcc                                                                                                                                                                                                              >  1:293981/9‑251 (MQ=255)
                               cgtcgacaccttccGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCCGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCt                                                                                                                                                                                                             <  1:473921/239‑1 (MQ=255)
                               cgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACCCCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCt                                                                                                                                                                                                             <  1:206825/244‑1 (MQ=255)
                                          cccGTAGTTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCGTGCAGTCGCCCGGCGAGGGCTTCCAGCGGAGCGGGGTTGAGGCTGTCGGTAAACACCTGTGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTc                                                                                                                                                                                                  <  2:596114/251‑1 (MQ=255)
                                            cGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTcc                                                                                                                                                                                                 >  1:545451/1‑250 (MQ=255)
                                               aGCTCTTTTTCCTTGTCCGGGAAGTCATCCCACCGCGCCTCGCCCTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTc                                                                                                                                                                                             <  2:451381/251‑1 (MQ=255)
                                                   ctTTGTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTccc                                                                                                                                                                                           <  2:334180/249‑1 (MQ=255)
                                                      tttCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCt                                                                                                                                                                                       >  1:439276/1‑250 (MQ=255)
                                                               ccGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACgg                                                                                                                                                                             <  1:345929/251‑1 (MQ=255)
                                                                       tCAACCAACAGCGCTTCGCACTCTGGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGTCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCAACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGCGCCTGACCGAATTTCGATTCCCAGCTACTCTGACGGGCAAAGGt                                                                                                                                                                     >  2:221577/1‑251 (MQ=255)
                                                                         aaccaacAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAgggtt                                                                                                                                                                   >  1:381111/1‑248 (MQ=255)
                                                                           ccaccAGCCCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGTGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCTGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCAGCCAGCAGATGCGAGAAGGCCGGAGCCTGACCGAAGTTCCATTCCCAGTTATTCTGACGGGCAAATGTTTCg                                                                                                                                                                 <  2:269690/247‑1 (MQ=255)
                                                                              acAGCGCTGCGCACTCGGGTTGCAGCATACCCGCGCGGGACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCATCCAGCGGCGCGGGGTTGAGGCTGTCGGTAACCACCTGGGCGCTGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAACTAAAGCGTTCATCCAGCATATGCGAGAATGCCGTGGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCg                                                                                                                                                              <  2:52185/251‑1 (MQ=255)
                                                                                    cgtCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGttt                                                                                                                                                        <  1:361422/249‑1 (MQ=255)
                                                                                      tCGCACTCGTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGttt                                                                                                                                                        <  2:283879/249‑1 (MQ=255)
                                                                                       cGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTgg                                                                                                                                                      <  1:267694/250‑1 (MQ=255)
                                                                                         taCTCCTGTTGCAGCATATCTGCGCGGTGCAAGCTGCCTTGCAGTCGTCCGCCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTGAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGCAAAGCAGTTCCACGTCACACCAGCTAAAGCGTTCATCCAGCAGATGAGAGACTGCCGGAGCCTGCCCTAATATCCATTCCCGGCTCCTCTGGCGGGCAAAGGTGTCGGCgtcgtttttgt                                                                                                                                                    >  2:21174/2‑240 (MQ=255)
                                                                                               tGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTg                                                                                                                                              <  2:155507/250‑1 (MQ=255)
                                                                                                      gcggatccgcgcGGTTCAGGCAGCATTGCAGCAGGCCGGCGAGGGCTGCCAGCGGCGCGTGGTTGAGGCTGTCGGTGAACACCTGGGCGCGGGTGATATGGCCTGTTTCAACGTCGAAAGGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCAGCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTg                                                                                                                                      <  2:191076/243‑1 (MQ=255)
                                                                                                         tatCCGCGCGGTCCAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAAATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAATCTGGCGTTTTGtcc                                                                                                                                   >  2:376925/1‑249 (MQ=255)
                                                                                                                     acgggCAGCCTTGCAGTCGCCCGGCGAGGGCTGCCAGCGGCGCGTGGTGGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCGTTTTCAACGTCGAAAGTCAGGTCCAGGCCGCCCCATCGATAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCGAACCGAAGTTCCAGTCCCAGCTACCCTGATGCGCAAAGGTGGCGGCGAATTTTGGCAAGTCTGGCGTTTTGTCCGGGGAgatgat                                                                                                                        <  1:32171/247‑1 (MQ=255)
                                                                                                                       aGGCAGCCTTGCAGTCTCCCGGCGCGGGCTTCGAGCGGCGCGGGGTTGAGGCTGTCGGGAAACACCCGGGCGCGGGTGATATGGCCGTTTTCAACGTCGAAATGCAGGTACACGCCGCCCCAGCTAAATCGTGCAACCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCATCGACTCTGAAGGGCAAAGGTTTCGGCGAAGTGTGGCAAGTCTGGCGGTTTTTCCGGGGAGTTGAttta                                                                                                                     <  2:421071/251‑2 (MQ=255)
                                                                                                                       aGGCAGCCTTGCAGGCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTc                                                                                                                     <  2:278554/251‑1 (MQ=255)
                                                                                                                           aGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTAATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCt                                                                                                                 >  1:94912/1‑251 (MQ=255)
                                                                                                                                 gcaggcgcccGGCGTGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACGCCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTctaca                                                                                                           <  2:225313/243‑5 (MQ=255)
                                                                                                                                  cAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTctacac                                                                                                          >  1:338344/1‑246 (MQ=255)
                                                                                                                                  cAGTCGCCCGGCGAGGGCCTCCAGCGGCGGTGGGTTGAGGCTGTCGGTATACTCCCTGGCGCGGGTGATATGTCTTCTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAGGCGTTCTTCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCTGCGAAGTTTGGCAAGTCGGGCGTTTTGTCCGGTGAGATGATTTCCTCTTctacac                                                                                                          <  2:295484/251‑6 (MQ=255)
                                                                                                                                          cgccGAGGGCTTCCAGAGGCGCGTGGATGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGACATGCAGTTCCACGCCGCACCAGCTAAAGCGTTCATCCAGCAGATGCGCGAATGCCGGAGCCTGACCGCAGTTCAATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGATTGGCAAGTATAGCGTATTGTCCGGGGAGATGATTTCCGcatctacacacccgcc                                                                                                   >  2:91461/4‑235 (MQ=255)
                                                                                                                                           ggCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGttt                                                                                                                                                                   >  2:513445/1‑185 (MQ=255)
                                                                                                                                           ggCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGttt                                                                                                                                                                   <  1:513445/185‑1 (MQ=255)
                                                                                                                                                          cggcgCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAgg                                                                                  >  1:514301/1‑251 (MQ=255)
                                                                                                                                                           ggcgCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGc                                                                                 >  1:286843/1‑251 (MQ=255)
                                                                                                                                                                   gTTGATGCTGTAGGTAAACACCTGGGCGCGGGAGATATGGCCTTTTTCAACGTCGACATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCCTTCCCAGATACTCTGACGGGCAAGGGTTTCGGca                                                                                                                                                              >  2:378471/1‑165 (MQ=255)
                                                                                                                                                                   gTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCg                                                                                                                                                              <  1:378471/166‑1 (MQ=255)
                                                                                                                                                                   gTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGtt                                                                          >  1:240777/1‑250 (MQ=255)
                                                                                                                                                                      gAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATgg                                                                      <  1:96869/251‑1 (MQ=255)
                                                                                                                                                                        ggCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGcc                                                                    >  1:320403/1‑251 (MQ=255)
                                                                                                                                                                            gTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTc                                                                  <  1:286819/249‑1 (MQ=255)
                                                                                                                                                                             ccGGTAAATACCCGGGCGCGGGTGATATGTCCTTTTTCAACGACGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCCTTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGCCCGTGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGTCTCGCa                                                               <  2:673131/250‑1 (MQ=255)
                                                                                                                                                                                  aaaCACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTCCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAACCCt                                                          >  1:254737/1‑251 (MQ=255)
                                                                                                                                                                                   aaCACCTGGGCGCGGGTGAGATGCCCTTTGTCAACTGCGAAATGCAGTCCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGGTTCGGCGAAGTTTGGCTAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCGACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTAGGGCCTCGCTAACCTg                                                         <  2:48782/251‑1 (MQ=255)
                                                                                                                                                                                            gagagggTGATATGGACTTTTGAAACGTCTAGAAGCAGTTCCACGCCGACCCAGCAAAAGCGTGCAGCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGGTCCATTCCCAGCTCCTCTGACTGGAGAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTGCCGCTTCTACACGCTCGCCATAAGGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGt                                                <  2:212243/247‑1 (MQ=255)
                                                                                                                                                                                               cgGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGAt                                             >  1:171663/1‑251 (MQ=255)
                                                                                                                                                                                                  gTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGTTTATGGCCTCGCAAACCTGCTCATGGGTGATccc                                          <  2:205669/251‑1 (MQ=255)
                                                                                                                                                                                                   tGATATGGCCTTTTTCAACGTAGAAATGCAGTTTCACGCCGCCCCAGCTAAAGCGTTCTTCCAGCAGATGCGAGACTGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCCAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGTTGATTTCCGCTTCTACACGCTCGCCCTAATGGACGCAAAAGGCCTTGTTTTTGGCCTCGCAACCCTGCTCATGTGTGCTccac                                         >  2:96547/1‑249 (MQ=255)
                                                                                                                                                                                                                 taaacgtctgaaTGCAGTTCCACGTCGCTTCAGCTAATGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATCCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGAGTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGTTTTGGTTATGGTCTCGCAAACCTGCTCATGGGTGATCCCTGGCAGCAGCTCgt                            <  1:546319/240‑2 (MQ=255)
                                                                                                                                                                                                                 tCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACAAGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCgg                            <  1:546301/250‑1 (MQ=255)
                                                                                                                                                                                                                  cAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTg                          <  1:221577/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGt                      >  1:641912/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaaTGCCGTCCCACGCCGCCCCAGCTAAAGCTTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTGTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTGCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCGCGGTGAGGTTGGt                  <  2:439276/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAACGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCGCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCACACCTGCTCATGGGTGATCCCCGGCAGCAGCACggggcgggtggtc                 >  2:224560/1‑240 (MQ=255)
                                                                                                                                                                                                                           aaTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTc                 >  1:564699/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTc                 >  1:46376/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTc                 >  1:371130/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCAcgc             >  1:405487/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGCAAAAGGCCTTGGTTATGGCTTCGCAAACCTGCTCCTGGGTGATCCCCGGCAGCAGCCCGGTGGGGTTGGTCCCGCgg           >  2:246815/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAAc       <  1:253548/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gccgccCCAGCTAAAGTGTTCAACCAGCAGATGCGAGAATGCCGGTGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAATGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgca    <  2:514301/251‑3 (MQ=255)
                                                                                                                                                                                                                                          cgccCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgcact  <  1:468321/251‑5 (MQ=255)
                                                                                                                                                                                                                                          cgccCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACgcact  <  1:376925/251‑5 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
CGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACC  >  NC_000913/4623094‑4623578

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