Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,309 A→G L49L (CTA→CTG dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,3090AG89.7% 99.9 / 2.4 41L49L (CTA→CTGdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (19/16);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAC  >  NC_000913/12067‑12552
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
gCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTc                                                                                                                                                                                                                                              <  1:373575/250‑1 (MQ=255)
  aaaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCCCCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGc                                                                                                                                                                                                                                           >  1:82046/1‑251 (MQ=255)
    aaaaaTGAAATTGGGCAGTTGAAATCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCg                                                                                                                                                                                                                                         >  1:22319/1‑251 (MQ=255)
     aaaaTGAACTTGGGCAGTTGAAACCAGCCGTTTCGCCCCTATTCAAGACTCACAACCACATGATGACCGAATATATAGTGGAGCCGTTTAGATGGTTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTCGCGCTTATGGATGGCAACACTCCTCGTGTGCTGGAGACCGCCGCATGCGCTCGCACCCCGCCTTCTATCATTGACTATACCCAGGATGGTGAAACTCTGGTTGGGCAACCgg                                                                                                                                                                                                                                        >  2:38890/1‑251 (MQ=255)
                  gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGc                                                                                                                                                                                                                           <  2:22319/251‑1 (MQ=255)
                   cagatGAAACCAGACGTTTCGACCATATTCCAGACTCACAAACACATGATGCCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCAGAGTACTACCAACTCTTGTGTAGCGATTATGGGTGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCGTCTATACTTGCCTATACCCCGGATGGAGAAACTCTGGTTGGTCAGCCGGCTAACCGtacggc                                                                                                                                                                                                                           >  2:502715/5‑245 (MQ=255)
                     gTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGt                                                                                                                                                                                                                        <  1:323199/251‑1 (MQ=255)
                          aaCCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGCAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTCGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACg                                                                                                                                                                                                                    >  2:234150/1‑250 (MQ=255)
                              aGACGTCTCGCCCTTATATCAGACTCACAACCACATGATGACCGAATAGATAGTGGAGACGTTTAGACGGGAAAAATAATGGGTATCGACCTGGGTACTGCCAACTCGTGTTTAGCGATTATGGATGGCACCACTCCTCGTGTGCTTGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTTGTTGGTCAGCCGGCTAAACGCCAGGCAGTGACGTAccc                                                                                                                                                                                                               <  2:128602/251‑1 (MQ=255)
                                     tCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCaaaa                                                                                                                                                                                                         <  1:318996/250‑1 (MQ=255)
                                             aTTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGtt                                                                                                                                                                                                >  1:491558/1‑251 (MQ=255)
                                                        acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGCTTAAAc                                                                                                                                                                                     >  1:278604/1‑251 (MQ=255)
                                                        acaACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGGAGCGATTAGGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAAc                                                                                                                                                                                     <  2:198251/251‑1 (MQ=255)
                                                              acaTGCTGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACACCGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGACTAACCGTCAGGCAGGGGCGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTg                                                                                                                                                                                >  2:151821/1‑250 (MQ=255)
                                                                                   atgagaCGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCATCTCTTGTGTAGCGATTATGGATGGCACCACTCCGCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTGGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGAGTGGCCGCCGCGTCCAGGACg                                                                                                                                                          <  2:282808/249‑1 (MQ=255)
                                                                                             tAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTAc                                                                                                                                                 >  1:92986/1‑250 (MQ=255)
                                                                                                  gggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                                                                                           >  1:106687/1‑251 (MQ=255)
                                                                                                   ggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTgg                                                                                                                                          >  1:37624/1‑250 (MQ=255)
                                                                                                   ggTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGa                                                                                                                                          <  1:25267/251‑1 (MQ=255)
                                                                                                    gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                                                                                                         <  2:92986/251‑1 (MQ=255)
                                                                                                        aataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTGCCAGGACGAAGAAGTACAGCGTGATGttt                                                                                                                                     <  1:156198/251‑1 (MQ=255)
                                                                                                               ggTATCGACCTGGGTACTACCAACTCTTGTTTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACGCCCTGTTGGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcatcat                                                                                                                              <  2:278604/251‑1 (MQ=255)
                                                                                                                    cGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGt                                                                                                                         >  1:50247/1‑251 (MQ=255)
                                                                                                                     gACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGtt                                                                                                                        <  1:376013/251‑1 (MQ=255)
                                                                                                                                        ctTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCCATCATTGCCTATACCCAGGATGGTTAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGctgctg                                                                                                       <  1:296631/249‑1 (MQ=255)
                                                                                                                                          tgtgTAGCGAGTATGGAGGGCACCACTCCTCCTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCGATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGGACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGa                                                                                                      <  2:204541/248‑1 (MQ=255)
                                                                                                                                            tgtAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGAACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACg                                                                                                 >  1:37806/1‑251 (MQ=255)
                                                                                                                                               aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                                                                                                               >  1:37835/1‑122 (MQ=255)
                                                                                                                                               aGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                                                                                                               <  2:37835/122‑1 (MQ=255)
                                                                                                                                                        ggatggCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGCTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATggg                                                                                     >  2:153494/1‑251 (MQ=255)
                                                                                                                                                                  acTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTa                                                                            >  1:506924/1‑250 (MQ=255)
                                                                                                                                                                        cgtgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGcc                                                                      <  1:217678/247‑1 (MQ=255)
                                                                                                                                                                                          gAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGCGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAACCACCCTGTTTGCGATTAAACGCCTGATTGGACGCCGCTTCCAGGACGAAGAAGAACAGCGTGATGTTTCCATCATGCCGTTCAACATTATTGCTGCTGATAACGGCGACGCATAGGTCGAAGTTAAAGGCCAGGAAATggacacgccgc                                                    >  2:48052/1‑241 (MQ=255)
                                                                                                                                                                                                 cacgacccacCCCTTCTAACATTGCCTATACCCGGGATGGTGAACCCCTGGTGGGTCAGCCGGCTAAACGTCAGGCGGTGACGAACCCGCAAACCACCCTGTTTGCGATAAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCAAGAAAATGGCACCGCCGCAGATctct                                            <  2:106687/245‑1 (MQ=255)
                                                                                                                                                                                                      accacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTctgctg                                        >  1:251452/1‑250 (MQ=255)
                                                                                                                                                                                                               tCTATCATTCCCTATCCCCAGGATGGTGAACCTCTAGTGGGTCAGCCGGCAAAACGTCAGGCAGTGACGACCCCGCAAAACCCTCTGTTTGCGATTAACCGCCTGTTTGGCCGCCGTTTCCAGGACGAAGAAGTACAGCGTGATTTTTCCATCATGCCGTTAAAAATTATTTCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGa                               <  1:234150/250‑1 (MQ=255)
                                                                                                                                                                                                                     agtGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGaagaag                                                                                                                                                    <  2:242261/125‑1 (MQ=255)
                                                                                                                                                                                                                     aTTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGaagaag                                                                                                                                                    >  1:242261/1‑127 (MQ=255)
                                                                                                                                                                                                                          cTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATga                      >  1:187465/1‑248 (MQ=255)
                                                                                                                                                                                                                                       ggTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGCTCTCTGCTGGAGTGCTGAAAAAAATGAAGAAAACCGCTGa         >  2:117177/1‑248 (MQ=255)
                                                                                                                                                                                                                                           aaaCTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTTATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGACGTTAAAGGCCAGAAACTGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAc  >  1:185339/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
GCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAC  >  NC_000913/12067‑12552

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