Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 18.0% intergenic (+19/+21) metE → / ← ysgA 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3330AT18.0% 424.7 / 83.3 184intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (71/79);  new base T (11/22);  total (82/102)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 187 total aligned reads displayed.

TGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATG  >  NC_000913/4015235‑4015431
                                                                                                  |                                                                                                  
tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                   >  2:1106307/1‑100 (MQ=255)
tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                   >  2:3672099/1‑100 (MQ=255)
 gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                  <  2:2545963/100‑1 (MQ=255)
   aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                <  2:3304341/100‑1 (MQ=255)
      cgcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                             >  1:3509693/1‑100 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                          <  1:3445698/100‑1 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                          <  2:3428738/100‑1 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                          <  1:1618391/100‑1 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                          <  2:3202123/100‑1 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                          <  2:2185889/100‑1 (MQ=255)
         gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGATGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                          <  1:1450443/100‑1 (MQ=255)
          cgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                         <  1:3170559/100‑1 (MQ=255)
           ggCTGGGCGCAAGCCCGCGCGGCACTGGGGAACATGGGGCCGGCGGCGCGGAACTTGCGTCGGGGGTAGCATCGTACCGGGGTGGTAGTACCAGCCgg                                                                                          <  2:4009852/98‑1 (MQ=255)
           ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTc                                                                                        <  1:864745/100‑1 (MQ=255)
            gCTGGCCAGAAACCCGCGCGGCACTGGCGAACTTGGTGCTGGGGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                       <  2:2941866/100‑1 (MQ=255)
            gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGttt                                                                                       <  1:3898458/100‑3 (MQ=255)
            gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGttt                                                                                       <  2:58936/100‑3 (MQ=255)
            gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                       >  2:2034221/1‑100 (MQ=255)
            gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                       <  1:702452/100‑1 (MQ=255)
             cTGGCCATAAACCCGCGCGGCACTGTTGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGCGGTAAAATCCAAACCGGGTGGTAATACCTCCCGGTCtt                                                                                      <  1:1040720/100‑1 (MQ=255)
             cTGGCCAGAACACCGCTCGGCACTGGCGCACATGGTGCCGGCCGCGCAGAACTTGCCTCGGTGGTAACATCCAAACCGGGTGGTATTACCACCCgc                                                                                          <  2:794622/96‑2 (MQ=255)
             cTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttg                                                                                      <  1:1731795/100‑4 (MQ=255)
              tGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                         <  2:1692084/96‑1 (MQ=255)
                gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                         <  2:3782276/94‑1 (MQ=255)
                gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                   <  1:3643635/100‑1 (MQ=255)
                gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                   <  1:3310362/100‑1 (MQ=255)
                 ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCATGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttggatt                                                                                  <  1:537506/100‑8 (MQ=38)
                 ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                  >  2:4187119/1‑100 (MQ=255)
                 ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                  >  2:1288519/1‑100 (MQ=255)
                  cAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                 <  1:169292/100‑1 (MQ=255)
                   aGAAAGCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                <  2:1775173/100‑1 (MQ=255)
                   aGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                <  2:2276457/100‑1 (MQ=255)
                   aGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                <  1:2059467/100‑1 (MQ=255)
                    gAAGCCGGCGCGGCTCTGGCGAAGATGGTGCAGGCGGCGCAGAACGTGCGTCGGGGGTAAAGTCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                 <  2:2447234/98‑1 (MQ=255)
                      aaCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTa                                                                             <  2:520921/100‑1 (MQ=255)
                        cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                           >  2:3236950/1‑100 (MQ=255)
                         ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                          >  2:246397/1‑100 (MQ=255)
                         ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                          >  1:2467256/1‑100 (MQ=255)
                         ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGATCTTTTCTCATTACAg                                                                          >  2:2449169/1‑100 (MQ=255)
                           gcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGcn                                                                        <  2:1278109/100‑2 (MQ=255)
                              cgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGAct                                                                     >  2:2774798/1‑100 (MQ=255)
                                 cACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACTACCCGGTCTTTTCTCATTACAGCGActtct                                                                  <  2:2756634/100‑1 (MQ=255)
                                   cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACttcttc                                                                <  1:792975/100‑1 (MQ=255)
                                    tGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTcc                                                               <  2:2349080/100‑1 (MQ=255)
                                       cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccac                                                            <  2:451538/100‑1 (MQ=255)
                                       cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccac                                                            <  2:1167538/100‑1 (MQ=255)
                                        gAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGTGTAAACTCCAAACCGGGTGGTAATACCGCCCGGTCTTTTCTCATTACAGCGACTTCTTCccacc                                                           <  2:3010067/100‑1 (MQ=255)
                                           cATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATa                                                        <  1:83757/100‑1 (MQ=255)
                                            aTGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                       <  1:1895427/100‑1 (MQ=255)
                                             tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACt                                                      <  2:673250/100‑1 (MQ=255)
                                              ggTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                     >  1:3799860/1‑100 (MQ=255)
                                               gTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCCAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGc                                                    <  1:3014491/100‑1 (MQ=255)
                                                 gCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCtt                                                  >  2:1535022/1‑100 (MQ=255)
                                                  cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTa                                                 >  1:2333099/1‑100 (MQ=255)
                                                    ggcgggGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                               <  2:2525792/100‑1 (MQ=255)
                                                    ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGCGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTATTCCCACCGTCATCCTTaaa                                               >  1:2869980/1‑100 (MQ=255)
                                                    ggcagcGCAGCACTTGCGTAGGGGGTAAAATCCAAAGCGGGAGGTAACACCATCCGGTCTTTTATCATAACAGCGACTTCTTACCCCCATCCTGATGaaa                                               >  2:1445974/1‑100 (MQ=255)
                                                     gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                              <  2:1318500/100‑1 (MQ=255)
                                                      cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                             <  1:4070011/100‑1 (MQ=255)
                                                        gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                           >  2:3104035/1‑100 (MQ=255)
                                                         cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAtt                                          <  1:1930016/100‑1 (MQ=255)
                                                         cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAtt                                          >  2:3161693/1‑100 (MQ=255)
                                                              agcgtGCGTCGGAGGTAAAATCTAAACCAGGTGGTAATGCCACCCGTTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAg                                      <  2:3966748/95‑1 (MQ=255)
                                                               aCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCa                                    >  1:2187372/1‑100 (MQ=255)
                                                                cTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCAt                                   >  1:2689628/1‑100 (MQ=255)
                                                                    cGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                               >  1:2758662/1‑100 (MQ=255)
                                                                     gTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCt                              >  1:1471756/1‑100 (MQ=255)
                                                                      tCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTg                             <  1:1351122/100‑1 (MQ=255)
                                                                       cGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgc                            <  1:1179423/100‑1 (MQ=255)
                                                                         ggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTATCATTACAGCGACTGCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcca                          >  2:2312030/1‑100 (MQ=255)
                                                                          gggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                         >  2:1806662/1‑100 (MQ=255)
                                                                           ggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagc                        <  1:964225/100‑1 (MQ=255)
                                                                            gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                       >  1:1797735/1‑100 (MQ=255)
                                                                            gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                       >  1:1249922/1‑100 (MQ=255)
                                                                             tAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACgccgcc                           >  2:2632572/1‑95 (MQ=255)
                                                                             tAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcca                      >  1:3957014/1‑100 (MQ=255)
                                                                                aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtct                   >  2:1138932/1‑100 (MQ=255)
                                                                                 gaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctt                  >  1:3086209/7‑100 (MQ=255)
                                                                                 gaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctt                  >  1:163319/7‑100 (MQ=255)
                                                                                  tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCCTACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                 >  2:333009/1‑100 (MQ=255)
                                                                                  tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                 <  2:1631448/100‑1 (MQ=255)
                                                                                    cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg               >  1:422945/1‑100 (MQ=255)
                                                                                    cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg               <  1:871618/100‑1 (MQ=255)
                                                                                    cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg               <  1:981741/100‑1 (MQ=255)
                                                                                     agaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc              >  1:3116911/3‑100 (MQ=255)
                                                                                     agaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc              >  2:2397247/3‑100 (MQ=255)
                                                                                     aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc              <  1:1535025/100‑1 (MQ=255)
                                                                                      aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa             <  2:387739/100‑1 (MQ=255)
                                                                                        ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa           >  1:3146181/1‑100 (MQ=255)
                                                                                        ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa           >  1:1485888/1‑100 (MQ=255)
                                                                                        ccGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg            >  2:3539301/1‑99 (MQ=255)
                                                                                          gggTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATg         <  1:3959084/100‑1 (MQ=255)
                                                                                             tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg      >  1:3297332/1‑100 (MQ=255)
                                                                                             tggtAATACCACCCGGTCTTTTCTCATTACAGCGACGTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg      <  1:3993756/100‑1 (MQ=255)
                                                                                              ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc     <  2:3509676/100‑1 (MQ=255)
                                                                                              ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc     <  2:1227638/100‑1 (MQ=255)
                                                                                                tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt   <  1:2405105/100‑1 (MQ=255)
                                                                                                tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt   >  1:1117124/1‑100 (MQ=255)
                                                                                                tAATACCACCCGGTCTTTTCTCATCACAGCGACTTCTTCCCACGATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt   <  2:2685025/100‑1 (MQ=255)
                                                                                                 aaTACCACCCGGGCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATg  >  1:1389540/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATG  >  NC_000913/4015235‑4015431

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: