Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 606,299 A→T 50.7% V362E (GTG→GAG)  ybdK ← weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,2990AT50.7% 2.5 / 291.5 193V362E (GTG→GAG) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base A (44/50);  new base T (57/40);  total (101/92)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 193 total aligned reads displayed.

ACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCG  >  NC_000913/606201‑606397
                                                                                                  |                                                                                                  
aCATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAc                                                                                                   >  1:1584837/1‑100 (MQ=255)
 cATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAcc                                                                                                  <  1:3710514/100‑1 (MQ=255)
 cATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCAc                                                                                                   >  1:1899882/1‑99 (MQ=255)
   ttGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAg                                                                                                <  2:3877386/100‑1 (MQ=255)
   ttGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAg                                                                                                <  1:1350911/100‑1 (MQ=255)
    tGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGc                                                                                               <  2:1536739/100‑1 (MQ=255)
     gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGcc                                                                                              >  2:753107/1‑100 (MQ=255)
     gTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGcc                                                                                              >  2:1507764/1‑100 (MQ=255)
      ttAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGccc                                                                                             >  2:1824829/1‑100 (MQ=255)
        aaaTGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCaa                                                                                           >  2:1960671/1‑100 (MQ=255)
        aaaTGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCaa                                                                                           <  1:1659636/100‑1 (MQ=255)
           tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCa                                                                                        >  1:3752183/1‑100 (MQ=255)
           tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCa                                                                                        >  1:2864108/1‑100 (MQ=255)
           tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCa                                                                                        >  2:1265767/1‑100 (MQ=255)
           tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCa                                                                                        <  2:1614683/100‑1 (MQ=255)
           tGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTGTTCTCCAGCCCAATCa                                                                                        >  1:2493285/1‑100 (MQ=255)
            gTAAGTTTAGACCGCATTGCCAGGTATATTCGTCCGTAACGCAATTCCCGGCTTCTCCACCGGCCCAGATCGCACCATGCTTTTTCACCAGCCCACTc                                                                                         <  1:4242823/98‑1 (MQ=255)
            gTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAg                                                                                       <  1:3701915/100‑1 (MQ=255)
            gTAAGTTAAAACACAATTGTCAGGGATATTCTTCTGTAGGGCAATTCACGGCTTTGTCACCGGCCCAGGTCTCACAATGCTTTTTCTCTAGCCCAATCa                                                                                        >  1:1031013/1‑99 (MQ=255)
               aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGa                                                                                    >  1:3906845/1‑100 (MQ=255)
               aGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGa                                                                                    >  2:3494040/1‑100 (MQ=255)
                gTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAg                                                                                   >  2:1634466/1‑100 (MQ=255)
                gTTAAAACACCATTGTCAGAAATATTCTTCTGTAAGGAAATTACCGGCTTAGTCACCGACCCAAATCTCAAAATGATTTTTATCAAACCAAATCAGCg                                                                                     >  1:472443/1‑98 (MQ=255)
                  tAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAgcc                                                                                 <  1:884984/100‑1 (MQ=255)
                  tAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCACTCAGCGAgcc                                                                                 >  2:3140030/1‑100 (MQ=255)
                   aaaaCACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCATTTTCACCAGCCCAATCAGCGAgccg                                                                                <  2:2597711/100‑1 (MQ=255)
                      acacCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCa                                                                             >  2:3229913/1‑100 (MQ=255)
                           aTTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGc                                                                        <  2:548813/100‑1 (MQ=255)
                             tGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcga                                                                      <  1:2979940/100‑1 (MQ=255)
                             tGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcga                                                                      >  1:3082214/1‑100 (MQ=255)
                                cAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGcgacga                                                                   >  2:3623088/1‑100 (MQ=255)
                                 aGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGaa                                                                  >  2:4064441/1‑100 (MQ=255)
                                  gggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGaaa                                                                 >  2:1706550/1‑100 (MQ=255)
                                  gggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGaaa                                                                 >  1:3855153/1‑100 (MQ=255)
                                  gggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGaaa                                                                 <  1:206574/100‑1 (MQ=255)
                                   ggATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAAt                                                                >  2:440659/1‑100 (MQ=255)
                                     atatTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcg                                                              <  1:2786087/100‑1 (MQ=255)
                                      tatTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATcgc                                                             >  1:3018564/1‑100 (MQ=255)
                                         tcttctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCACTCAGCGAGCCGCCATAGGCGACGAAATCGCGCa                                                          >  2:2485225/1‑100 (MQ=255)
                                            tctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTa                                                       >  2:205314/1‑100 (MQ=255)
                                            tctGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTa                                                       >  1:2135373/1‑100 (MQ=255)
                                              tGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGc                                                     >  2:1393017/1‑100 (MQ=255)
                                               gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCt                                                    >  2:2156629/1‑100 (MQ=255)
                                               gTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAACCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCt                                                    <  1:2626979/100‑1 (MQ=255)
                                                 aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                                  >  1:2609194/1‑100 (MQ=255)
                                                 aaGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgc                                                  <  2:298499/100‑1 (MQ=255)
                                                  aGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgcg                                                 >  1:3093500/1‑100 (MQ=255)
                                                   ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgcgc                                                >  1:2173672/1‑100 (MQ=255)
                                                   ggCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCtaccc                                                >  1:3082398/1‑96 (MQ=255)
                                                    tcatttCCCGGCTTAGCGACCGGCCCAGACCACACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTgcg                                                 <  1:683287/94‑1 (MQ=255)
                                                    gCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCt                                               <  2:3519430/100‑1 (MQ=255)
                                                     cAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCtt                                              <  2:2405996/100‑1 (MQ=255)
                                                       aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                            <  2:79159/100‑1 (MQ=255)
                                                       aTTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttca                                            <  2:25694/100‑1 (MQ=255)
                                                        ttCCCGGCTTAGTCACCGGCCCAGCTCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcat                                           >  1:400137/1‑100 (MQ=255)
                                                        ttCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcat                                           <  1:1314099/100‑1 (MQ=255)
                                                          cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattc                                         >  2:985467/1‑100 (MQ=255)
                                                          cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattc                                         <  2:136724/100‑1 (MQ=255)
                                                          cccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattc                                         <  2:31782/100‑1 (MQ=255)
                                                          cccGGCTTAGTCACCGGCCCAGATCCCACAATGCTCTTCCCCCAGCCCAATCAGCGAGCCGCCATCGGCGACGACATCGCGCATTAGCTGCGCttcattc                                         <  1:3148612/100‑1 (MQ=255)
                                                           ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                        <  1:293759/100‑1 (MQ=255)
                                                           ccGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                        >  2:3944255/1‑100 (MQ=255)
                                                           atctcTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattca                                        <  1:2965367/96‑1 (MQ=255)
                                                            cGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAg                                       <  1:1087812/100‑1 (MQ=255)
                                                               cTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCttcattcaagat                                    >  2:2965379/1‑96 (MQ=255)
                                                                 tAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGc                                  <  1:1724325/100‑1 (MQ=255)
                                                                   gTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTg                                >  1:2442271/1‑100 (MQ=255)
                                                                      aCCGGCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg                             <  1:1060881/100‑1 (MQ=255)
                                                                         ggCCCAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGAcc                          <  1:2135133/100‑1 (MQ=255)
                                                                         ggCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGAcc                          <  1:3231315/100‑1 (MQ=255)
                                                                         ggCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGAcc                          <  2:3543283/100‑1 (MQ=255)
                                                                          gCCCAGATCTCACAATGCTTTTTTTCTATCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCt                         >  1:1386918/1‑100 (MQ=255)
                                                                            ccAGATCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTgg                       >  2:1863233/1‑100 (MQ=255)
                                                                                aTCTCACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGAt                   >  1:4185602/1‑100 (MQ=255)
                                                                                aTCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGAt                   <  1:1769037/100‑1 (MQ=255)
                                                                                   tcACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCa                >  2:2257010/1‑100 (MQ=255)
                                                                                   tcACAATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCa                <  1:4085516/100‑1 (MQ=255)
                                                                                    cacaATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAg               >  2:3428310/1‑100 (MQ=255)
                                                                                     acaATGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAgg              <  2:370894/100‑1 (MQ=255)
                                                                                       aaTGCTTTTTCCCCAGCCAAAAAAGCGAGCCGCCATCGGCGACGAAATCGCGCATGAGGTGAGCTTCATTCAGACCGCTGACGATCTGGCGATGCAggg             <  1:1172418/99‑1 (MQ=255)
                                                                                       aaTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAggg             <  2:2470202/99‑1 (MQ=255)
                                                                                        aTGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  2:2014086/100‑1 (MQ=255)
                                                                                        aTGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  2:1800036/100‑1 (MQ=255)
                                                                                        aTGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  1:3975685/100‑1 (MQ=255)
                                                                                        aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  2:3093512/100‑1 (MQ=255)
                                                                                        aTGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  2:890279/100‑1 (MQ=255)
                                                                                         tGCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGcc           <  1:1208805/99‑1 (MQ=255)
                                                                                         tGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCt          <  2:1936884/100‑1 (MQ=255)
                                                                                          gCTTTTTCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTc         >  2:1212961/1‑100 (MQ=255)
                                                                                          gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTc         <  2:392785/100‑1 (MQ=255)
                                                                                          gCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTc         <  2:1664289/100‑1 (MQ=255)
                                                                                            tttttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGAGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCaa       >  1:3547601/1‑100 (MQ=255)
                                                                                            tttttCACCAGCCCAATCAGCGAACCGCCATCGACAACGAAATCGCACATTATATACGCTTAATTCAAACCACTAAAGACATTACGATGCAgaacctaaa       >  1:7668/1‑92 (MQ=255)
                                                                                             ttttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt      >  2:2662560/1‑100 (MQ=255)
                                                                                             ttttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt      >  1:4145110/1‑100 (MQ=255)
                                                                                             ntttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAAt      >  2:2949399/2‑100 (MQ=255)
                                                                                              tttCTCCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATc     <  2:1099159/100‑1 (MQ=255)
                                                                                              tttCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATc     <  2:3829865/100‑1 (MQ=255)
                                                                                                 cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATcgcg  >  1:4154224/1‑100 (MQ=255)
                                                                                                 cACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCACACCGCTGACGACCTGGCGATGCAGGGCCTCAATcgcn  <  1:2289296/100‑2 (MQ=255)
                                                                                                  |                                                                                                  
ACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACGACCTGGCGATGCAGGGCCTCAATCGCG  >  NC_000913/606201‑606397

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