Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 15.9% 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,3330AT15.9% 355.1 / 50.9 131intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (54/56);  new base T (14/7);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGA  >  NC_000913/4015235‑4015432
                                                                                                  |                                                                                                   
tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                    >  1:1236272/1‑100 (MQ=255)
 gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                   <  2:802787/100‑1 (MQ=255)
  gtgccgcgcgggtttctggccagccgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                  >  1:516649‑M2/31‑100 (MQ=255)
  aaaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                  >  1:1100033/1‑100 (MQ=255)
   aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                 <  2:854300/100‑1 (MQ=255)
    gccgcgcgggtttctggccagccgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacca                                                                                                >  2:365928‑M2/29‑100 (MQ=255)
    gccgcgcgggtttctggccagccgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGCCGGTAATacca                                                                                                >  2:2612826‑M2/29‑100 (MQ=255)
    aaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGATGCAGGCGGCGCAGAACTTGCGTCGGGTGTAAAATCCAAAACGGGTGGTAATac                                                                                                  <  2:511734/98‑1 (MQ=255)
     ccgcgcgggtttctggccagccgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAatacccc                                                                                               >  2:417445‑M2/28‑98 (MQ=255)
      CGCGCGGCTGGCCAGAAACCCGCGCGGCACTGTCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                              >  1:1998956‑M1/1‑100 (MQ=255)
      CGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                              <  2:1992384‑M1/100‑1 (MQ=255)
      CGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                              <  1:1454409‑M1/100‑1 (MQ=255)
       GCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATaccacc                                                                                              <  1:2112216‑M1/99‑1 (MQ=255)
         gcgGCTGGTCAGAAACCCGTGCTGCACTTGCTAATATGGTGTAGGTGGTGTAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                           <  1:863134/100‑1 (MQ=255)
         GCGGCTGGCCAGTAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           >  1:953933‑M1/1‑100 (MQ=255)
         GCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           >  1:1527429‑M1/1‑100 (MQ=255)
         GCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           >  1:231692‑M1/1‑100 (MQ=255)
           GGCTGGCCAGAAACCCGCGCGGCTCTGGCTAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                          >  1:2348893‑M1/1‑99 (MQ=255)
            GCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                        <  2:1306636‑M1/100‑1 (MQ=255)
             cTGGCCAGAAACCCGCGCGGCACTTGTGAACATGGTGTAGGTGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCGCCCGGttt                                                                                        <  1:48012/99‑3 (MQ=255)
             cTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAACATCCAAACCGGGTGGTAATACCACCCGGTCtt                                                                                       <  1:2514852/100‑1 (MQ=255)
              tGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttgg                                                                                      <  1:1428871/100‑5 (MQ=255)
                gtcAGAACGCCGCGCGGCACTGGCGAACGTTTTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttggat                                                                                    <  2:2063691/98‑7 (MQ=37)
                gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttggat                                                                                    <  2:2263214/100‑7 (MQ=255)
                gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                    <  1:2316252/100‑1 (MQ=255)
                 ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTtct                                                                                   >  1:44673/1‑100 (MQ=255)
                  cacgaccccGCGAGGCACTGGAGAACGTGGTGCAGGCGGCGCAGAACTTGAGTCGGGGGTAAAATCCCAACCGGGTCGTAATACCCCCCGGTCTTTtct                                                                                   <  2:2385812/93‑1 (MQ=255)
                  cAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                  <  2:1351066/100‑1 (MQ=255)
                  cAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                  <  2:597358/100‑1 (MQ=255)
                   aGAAACCCGCGCGGCACCGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                 <  1:2546783/100‑1 (MQ=255)
                    gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCAt                                                                                >  1:2227144/1‑100 (MQ=255)
                    gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCAt                                                                                >  1:1926203/1‑100 (MQ=255)
                        cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCATTACa                                                                            >  1:2033604/1‑100 (MQ=255)
                          cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGc                                                                          <  2:1154965‑M2/94‑1 (MQ=255)
                          cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGc                                                                          <  2:1182204‑M2/94‑1 (MQ=255)
                           gcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCg                                                                         <  2:665140‑M2/95‑1 (MQ=255)
                            cgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGCGGTAATACCACCCGGTCTTTTCTCATTACAGCGa                                                                        >  1:1721996‑M2/5‑100 (MQ=255)
                                gCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtc                                                                    >  1:2367035/1‑100 (MQ=255)
                                 cACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtct                                                                   <  2:1149478/100‑1 (MQ=255)
                                   cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACttcttc                                                                 <  2:2446556/100‑1 (MQ=255)
                                     ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTccc                                                               <  2:896761/100‑1 (MQ=255)
                                        gAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTTTTTTCTCATTACAGCGACTTCTTCccacc                                                            <  2:1482141/100‑1 (MQ=255)
                                         aaCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccacca                                                           <  2:535181/100‑1 (MQ=255)
                                             tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACt                                                       <  2:151123/100‑1 (MQ=255)
                                              ggTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCAGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                      >  1:822863/1‑100 (MQ=255)
                                                 tctggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCtt                                                   <  2:272583/97‑1 (MQ=255)
                                                 gCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCtt                                                   >  1:1914917/1‑100 (MQ=255)
                                                    ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaga                                                >  1:272585/1‑98 (MQ=255)
                                                    ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                >  1:1346528/1‑100 (MQ=255)
                                                    ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                >  1:1991578/1‑100 (MQ=255)
                                                     gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                               <  1:365930/100‑1 (MQ=255)
                                                     gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                               >  1:1051470/1‑100 (MQ=255)
                                                      cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                              <  2:1527400/100‑1 (MQ=255)
                                                        gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAACGACTTCTTCCCACCATACTGCTTAAACCAt                                            >  1:2150962/1‑100 (MQ=255)
                                                          gcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTc                                          <  2:2102442/100‑1 (MQ=255)
                                                             gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAg                                       >  2:1293468/1‑100 (MQ=255)
                                                                cTTGCGTCGGGGGTAAAATTCTAACCTTGTGGTAATATCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGc                                      <  1:609752/98‑1 (MQ=255)
                                                                    cGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                                >  1:2363411/1‑100 (MQ=255)
                                                                      tCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTg                              <  1:316502/100‑1 (MQ=255)
                                                                       cGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgc                             <  1:2475451/100‑1 (MQ=255)
                                                                        gggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcc                            >  1:1978305/1‑100 (MQ=255)
                                                                        gggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcc                            >  2:1331690/1‑100 (MQ=255)
                                                                          gggTAAAATCCAAACCGGGTGGTAATACCGCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                          <  2:1236260/100‑1 (MQ=255)
                                                                          gggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                          >  2:2182988/1‑100 (MQ=255)
                                                                           ggTAAAATCCAAATCGGGTGGTAATACCACCCGGTCTTTTCTCATTATAGCGAATTCTTCGCACCATACTGCTTAAACCATTCCAGCATACGCTgccagc                         >  1:2647886/1‑100 (MQ=255)
                                                                             tAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcca                       <  2:2014567/100‑1 (MQ=255)
                                                                             tAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcca                       >  1:201392/1‑100 (MQ=255)
                                                                             tAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcca                       <  2:516645/100‑1 (MQ=255)
                                                                              tgagaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAt                      >  2:2198313/10‑100 (MQ=38)
                                                                               aaaTCCAAACCTTTTATTAATACCACCCTTTCTTTTCTCATTACATCTACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtc                     <  1:1525851/100‑1 (MQ=255)
                                                                               aaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtc                     >  1:1231708/1‑100 (MQ=255)
                                                                               aaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtc                     >  1:403588/1‑100 (MQ=255)
                                                                                agaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtct                    >  2:2144231/8‑100 (MQ=38)
                                                                                aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTACTTAAACCATTCCAGCATACGCTGCCAGCCAtct                    >  1:2145800/1‑100 (MQ=255)
                                                                                 gaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctt                   >  2:278620/7‑100 (MQ=255)
                                                                                 gaaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctt                   >  2:21662/7‑100 (MQ=255)
                                                                                 aTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctt                   >  2:1012198/1‑100 (MQ=255)
                                                                                  tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                  <  1:1992403/100‑1 (MQ=255)
                                                                                  tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                  <  1:1923187/100‑1 (MQ=255)
                                                                                  aaaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                  >  1:1774454/6‑100 (MQ=255)
                                                                                   ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTTATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttct                 <  1:1597505/100‑1 (MQ=255)
                                                                                   ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCTCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttct                 <  1:1507089/100‑1 (MQ=255)
                                                                                   ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttct                 >  2:1965994/1‑100 (MQ=255)
                                                                                    cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTCCTTAAACCATTCAAAAATACACTGACAGCCATCTTCTg                >  1:526813/1‑100 (MQ=255)
                                                                                    aagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCAAtcttct                 >  1:1946603/4‑99 (MQ=255)
                                                                                     aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  1:264269/100‑1 (MQ=255)
                                                                                     aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  1:2050585/100‑1 (MQ=255)
                                                                                      gaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa              <  2:1081384/99‑1 (MQ=255)
                                                                                      gaCCGGGAGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa              >  1:1081394/2‑100 (MQ=255)
                                                                                       aCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg             <  2:2203231/100‑1 (MQ=255)
                                                                                       aCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg             <  2:1721980/100‑1 (MQ=255)
                                                                                       aCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg             >  1:2548857/1‑100 (MQ=255)
                                                                                        ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa            >  1:449958/1‑100 (MQ=255)
                                                                                        ccGGGTGGTATTACCACCCGGTCTTTTATCATTACAGCGACTAATTCCAACCATACTCCTCCAACCATTCCAGCATACGCTCCCAGCCATCTCCTCCAGa            >  2:374877/1‑100 (MQ=255)
                                                                                         cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTGAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGAt           >  1:2528865/1‑100 (MQ=255)
                                                                                              ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc      >  1:1152129/1‑100 (MQ=255)
                                                                                               gtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa     >  2:1719023/1‑100 (MQ=255)
                                                                                                tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt    >  2:1006796/1‑100 (MQ=255)
                                                                                                naaTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt    >  2:168143/2‑100 (MQ=255)
                                                                                                  aTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGa  >  2:2275350/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGA  >  NC_000913/4015235‑4015432

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