Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,821 A→G T99T (ACT→ACC lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,8210AG98.0% 160.2 / ‑4.5 51T99T (ACT→ACClplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/20);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

TAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCA  >  NC_000913/4623583‑4624065
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
tAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGTGCCCAATCGACGTGGAGATGGTTTTATCggac                                                                                                                                                                                                                                          >  1:58619/1‑248 (MQ=255)
    gCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACTCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTcc                                                                                                                                                                                                                                       <  1:46772/250‑1 (MQ=255)
        tgcgccgcCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGc                                                                                                                                                                                                                                    <  2:138294/247‑1 (MQ=255)
              gccAGTTTCTTTTTATCCGGAGTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTGgccagc                                                                                                                                                                                                                            <  2:138339/251‑1 (MQ=255)
               ccAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGGTTTATCGTACTCCGGCtttcccgcc                                                                                                                                                                                                                           >  1:288666/1‑244 (MQ=255)
                  gTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCCCCAGATCGTTACGCCCGGACGCCTCGGCGCTGGCGCCCAGCGCGTTCAGCGCATTGAGCACAATCGCCGTGGGGATGGTTTTATCGTACTCCGGCTTGCCCGCCAt                                                                                                                                                                                                                         >  2:109956/1‑250 (MQ=255)
                                 gATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCCTAAAGGTAAAGCAggg                                                                                                                                                                                                          >  1:304941/1‑249 (MQ=255)
                                               gCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGt                                                                                                                                                                                            <  2:289518/250‑1 (MQ=255)
                                               gCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGt                                                                                                                                                                                            >  1:289518/1‑250 (MQ=255)
                                                ccgggcggcTAGGGTCGGCATTGAGCAGCAAGGGGCGGTGGGGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGTCTTTGCGGTCGCCTTCGGCGGTTGTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCGGCGCATGGAGCCCAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTgg                                                                                                                                                                                          <  2:403962/248‑1 (MQ=255)
                                                ccAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGCGATCGTgg                                                                                                                                                                                          >  1:109812/1‑251 (MQ=255)
                                                      ggcTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGcc                                                                                                                                                                                                             >  1:70630/1‑225 (MQ=255)
                                                      ggcTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCa                                                                                                                                                                                                             <  2:70630/226‑1 (MQ=255)
                                                                ggCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCAc                                                                                                                                                                          <  1:432474/251‑1 (MQ=255)
                                                                                      tggtggAAGCCACGATCTTTGGTTTCGAGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGCCGCTTCGGCGCTGACGCCAAGCGCGTTCCACGCATTGAGAACAATCGACGGGGAGATGGTTTTATCGTACTCCGGCTTGCCCGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCAACTGCGCCGGGCCAGGCg                                                                                                                                                    >  1:171556/1‑251 (MQ=255)
                                                                                               ccacgATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGGTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTTCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTAtt                                                                                                                                           >  1:205130/4‑250 (MQ=255)
                                                                                                      cTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCCCCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCCAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCCGGCGGGCGTTATTTTCTTcc                                                                                                                                    >  1:165984/1‑251 (MQ=255)
                                                                                                               cgcgATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCCAGCGCGTTCAGCGCATTTAGCACAATCGACGGGGAGATGGTTTTATCGTACTCCCGCTTGCCCGCCATACAGGGAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCCCCGCTACTGCGACGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCgg                                                                                                                           >  1:468380/1‑251 (MQ=255)
                                                                                                                 cgATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTa                                                                                                                                                                                                                  <  1:459400/162‑1 (MQ=255)
                                                                                                                         cgatccTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCCGATCGTTACGCCCCGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCAATGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCCGGCCAGGCGGACGGTATCTTCTTCCATCCGCCGGGTATTTCAct                                                                                                                 >  1:166404/5‑251 (MQ=255)
                                                                                                                         cgatccTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGGGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCCCGCCGCCCCCGCTACTGCGCCGTGCCAGGCGGACGTTATTTTCTTCCCTCCGCCGGGTATTacacc                                                                                                                 >  2:340556/5‑250 (MQ=255)
                                                                                                                             ccTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGGGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCttt                                                                                                             >  1:299678/1‑251 (MQ=255)
                                                                                                                                      ttGCGGTCGCCTTCGTCGGTTTGCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGt                                                                                                     <  1:191981/250‑1 (MQ=255)
                                                                                                                                        gCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTc                                                                                                   <  2:398472/250‑1 (MQ=255)
                                                                                                                                              gCCTTCGGCGGTTTTCACCACCAGATTGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCCGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTgcgcgc                                                                                            >  1:355354/1‑251 (MQ=255)
                                                                                                                                                    ggcGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCCCGGGTTCTGCGCGCGCCCa                                                                                       >  2:63562/3‑250 (MQ=255)
                                                                                                                                                              accaccAGATAGTTACGCCCGGCCGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGATCTCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGCCCAATCCCTACCgt                                                                            >  2:70846/1‑251 (MQ=255)
                                                                                                                                                                   caGATCTTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCgg                                                                       <  1:454594/251‑1 (MQ=255)
                                                                                                                                                                      aTCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCAt                                                                    <  1:70846/251‑1 (MQ=255)
                                                                                                                                                                      aTCGGTACGCCCGTACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGACAGCCATAAAGGTAAAGCAGGTATTGCCGAGATAGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCCCTACCGTGTCGGCAt                                                                    >  2:409850/1‑251 (MQ=255)
                                                                                                                                                                         gatACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGCGATGGGTTTATAGTACTCCGGCTTGCCAGCCATAAAGGGAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCAACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGATTCTGCGCGCGACCAATCACTACAGTGCTGGCATTgc                                                                 >  2:285257/3‑251 (MQ=255)
                                                                                                                                                                           tACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTgcg                                                                <  1:351784/250‑1 (MQ=255)
                                                                                                                                                                                      gCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCCATCACTACCGTGTCGGCATTGCGCCAGAGAAAc                                                      >  1:434238/1‑249 (MQ=255)
                                                                                                                                                                                      gCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGACGGCATTGCGCCAGAGAAAc                                                      >  2:391751/1‑249 (MQ=255)
                                                                                                                                                                                         tCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAAc                                                 >  1:238165/1‑251 (MQ=255)
                                                                                                                                                                                               cTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGCCGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATCAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTCTTACACTCtt                                                                                                              >  2:250168/1‑184 (MQ=255)
                                                                                                                                                                                               cTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCtt                                                                                                              <  1:250168/184‑1 (MQ=255)
                                                                                                                                                                                                   cGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTg                                       >  2:253141/1‑251 (MQ=255)
                                                                                                                                                                                                    gccaaGCGCGTTCAGCGCATGGAGCACAATCGACGTGGAGATGGTTTTATCGGTCGCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTgg                                      <  2:225178/247‑1 (MQ=255)
                                                                                                                                                                                                      caaGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCg                                    >  1:201684/3‑251 (MQ=255)
                                                                                                                                                                                                        aGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCggg                                  >  2:54817/1‑251 (MQ=255)
                                                                                                                                                                                                            cgTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCAtc                              >  2:425347/1‑250 (MQ=255)
                                                                                                                                                                                                             gTTCAGCGCATTGAGCACAATCGACGTGGAGATGGGTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGTCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCAtct                             <  1:14487/251‑3 (MQ=255)
                                                                                                                                                                                                             gTTCAGCGAATTGAGCACAATCGACGGGGAGATGGTTTTATCGTACCCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCAtct                             <  2:355354/251‑3 (MQ=255)
                                                                                                                                                                                                                  gcgcATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCAGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGa                        <  1:147779/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggACGGGGAGATGTTTTTACCGTACCCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCa        <  2:299678/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGCCGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATcact                                                                                 >  2:383249/1‑178 (MQ=255)
                                                                                                                                                                                                                                  cGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATcact                                                                                 <  1:383249/178‑1 (MQ=255)
                                                                                                                                                                                                                                  cGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGACCGCGCTGCGTGGCGGGCATTTGGCGAAAATTACACtcctccc        >  1:333514/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACCCCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCAcc      <  2:288666/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTCCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCCTTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAACAATACACCCCTCCAcccccc  >  2:336341/1‑247 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
TAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCA  >  NC_000913/4623583‑4624065

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