Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 16.3% 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,3330AT16.3% 285.3 / 43.3 116intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (43/54);  new base T (9/10);  total (52/64)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 122 total aligned reads displayed.

GAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGA  >  NC_000913/4015236‑4015432
                                                                                                 |                                                                                                   
gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                   <  2:1157523/100‑1 (MQ=255)
  aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                 <  1:2313731/100‑1 (MQ=255)
     cgcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                              <  1:2487141/100‑1 (MQ=255)
        gcgGCTTACCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           <  2:2335110/100‑1 (MQ=255)
        gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           <  2:229791/100‑1 (MQ=255)
        gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           <  1:2207704/100‑1 (MQ=255)
        gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           <  1:2026138/100‑1 (MQ=255)
        gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                           <  1:1712504/100‑1 (MQ=255)
        gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                           <  1:1833349/100‑1 (MQ=255)
         cgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                          <  1:1716520/100‑1 (MQ=255)
         cgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                          <  1:703915/100‑1 (MQ=255)
           gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                        <  1:62152/100‑1 (MQ=255)
           gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                        >  1:1441553/1‑100 (MQ=255)
           gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCt                                                                                        <  1:727094/100‑1 (MQ=255)
               gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                    >  2:754381/1‑100 (MQ=255)
               gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCAAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                    <  1:1725112/100‑1 (MQ=255)
                ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                   >  1:1539211/1‑100 (MQ=255)
                ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                   <  1:487773/100‑1 (MQ=255)
                ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                   >  2:677810/1‑100 (MQ=255)
                ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                   >  1:1291310/1‑100 (MQ=255)
                 cAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                  <  1:1674010/100‑1 (MQ=255)
                  aGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                 <  1:195028/100‑1 (MQ=255)
                   gAACCCCCCGCGGCACTGGCGAACATGGTGCAGGCGGCGCACAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCAt                                                                                <  2:1907483/100‑1 (MQ=255)
                       cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                            <  2:2358094/100‑1 (MQ=255)
                        ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                           <  1:2150335/100‑1 (MQ=255)
                         cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGc                                                                          <  1:1963443/100‑1 (MQ=255)
                           cgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCg                                                                         >  2:938770/1‑99 (MQ=255)
                                cACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtct                                                                   >  2:501887/1‑100 (MQ=255)
                                  cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACttcttc                                                                 <  1:517930/100‑1 (MQ=255)
                                    ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTACAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACATCATccc                                                               >  2:339003/1‑100 (MQ=255)
                                    ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTccc                                                               <  1:1108194/100‑1 (MQ=255)
                                     gCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCcca                                                              <  1:458544/100‑1 (MQ=255)
                                      cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCcca                                                              >  2:2311659/1‑99 (MQ=255)
                                      cGAACATGGTGCAGGCGGCGCAGAACTTGCGCCGGGGGTAAAATCCAAACCGCGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccac                                                             <  2:2388399/100‑1 (MQ=255)
                                       gAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccacc                                                            <  1:591803/100‑1 (MQ=255)
                                           aTGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                        <  2:226016/100‑1 (MQ=255)
                                            tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACt                                                       >  2:110392/1‑100 (MQ=255)
                                            tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACt                                                       <  1:2537756/100‑1 (MQ=255)
                                             ggTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                      >  2:890737/1‑100 (MQ=255)
                                               tGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                    <  1:1456565/100‑1 (MQ=255)
                                               tGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCAAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                    >  1:1807751/1‑100 (MQ=255)
                                                   ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                >  2:372709/1‑100 (MQ=255)
                                                   ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                >  1:2015340/1‑100 (MQ=255)
                                                   ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                <  2:2124357/100‑1 (MQ=255)
                                                    gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                               <  2:2097212/100‑1 (MQ=255)
                                                     cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAc                                               <  2:2334070/99‑1 (MQ=255)
                                                     cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                              <  2:1390420/100‑1 (MQ=255)
                                                     cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                              <  2:1807737/100‑1 (MQ=255)
                                                     cggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                              >  1:2428021/1‑100 (MQ=255)
                                                       gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                            >  1:1854973/1‑100 (MQ=255)
                                                       gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                            <  1:1172948/100‑1 (MQ=255)
                                                       gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                            >  1:1854967/1‑100 (MQ=255)
                                                        cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAtt                                           <  1:890745/100‑1 (MQ=255)
                                                        cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAtt                                           <  2:1127033/100‑1 (MQ=255)
                                                        cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                            >  2:251784/1‑99 (MQ=255)
                                                        cgcAGAACTTGCGTCGGGGGAAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCAATACAACGAATTCTTCCCACCATACTGCTTAAACCAtt                                           >  2:1376943/1‑100 (MQ=255)
                                                           aGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTTTTTTTTTTTTTATTTCTTTTCACCATACTGCTTAAACCATTc                                          >  2:1197726/1‑98 (MQ=255)
                                                              aCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCa                                     <  1:2530281/100‑1 (MQ=255)
                                                              aCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCa                                     <  2:1688237/100‑1 (MQ=255)
                                                                  gCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACg                                 >  2:263816/1‑100 (MQ=255)
                                                                   cGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                                >  2:606327/1‑100 (MQ=255)
                                                                   cGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                                >  1:238047/1‑100 (MQ=255)
                                                                         gggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                          <  1:1044118/100‑1 (MQ=255)
                                                                          ggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagc                         <  1:656107/100‑1 (MQ=255)
                                                                           gaaaaaTCTAAACCGGGAGGTAATACCACCCGGTCATATATCAATACATAACCTTCATCCACCCATCATGCTTAAACCATTCCAGCAAAAGCTgccagc                         >  2:2351033/3‑99 (MQ=255)
                                                                           gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                        >  1:2243987/1‑100 (MQ=255)
                                                                           gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                        >  2:1693271/1‑100 (MQ=255)
                                                                               aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtct                    <  1:2085533/100‑1 (MQ=255)
                                                                                 tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtctcc                  <  2:2524732/100‑3 (MQ=255)
                                                                                   aagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                >  1:2519513/4‑100 (MQ=255)
                                                                                    aaaCCGGGTGGTAATACCGTCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  2:286273/100‑1 (MQ=255)
                                                                                    aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  1:2127918/100‑1 (MQ=255)
                                                                                    aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  1:900747/100‑1 (MQ=255)
                                                                                    aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATGACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               <  2:206524/100‑1 (MQ=255)
                                                                                     tccgatctggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa              <  2:1503060/93‑1 (MQ=255)
                                                                                     gaCCGGGGGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTTCTTCAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc               >  2:1134503/2‑99 (MQ=255)
                                                                                     aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa              >  1:908956/1‑100 (MQ=255)
                                                                                     aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCa              <  2:2456180/100‑1 (MQ=255)
                                                                                      aCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTTCCACCCATCTTCTGCa              >  2:2013910/1‑99 (MQ=255)
                                                                                      aCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAg             <  2:2155694/100‑1 (MQ=255)
                                                                                      aCCGGGAGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAATGCTTAAACCATTCCAGAATACGCTGCAAGCCATCTTCTg                >  2:1791953/1‑97 (MQ=255)
                                                                                       ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa            >  2:2357535/1‑100 (MQ=255)
                                                                                       ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa            >  1:469265/1‑100 (MQ=255)
                                                                                       ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa            >  2:2521341/1‑100 (MQ=255)
                                                                                       ccGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTCCAGa            >  2:11291/1‑100 (MQ=255)
                                                                                           gtggtAATACCCCCCGGTCTTTTCTCATTACATCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg        <  2:1650021/100‑1 (MQ=255)
                                                                                           gtggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg        <  2:806151/100‑1 (MQ=255)
                                                                                            tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg       <  2:408994/100‑1 (MQ=255)
                                                                                             ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCaagatcgga      >  1:1503071/1‑92 (MQ=255)
                                                                                             ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc      >  2:554718/1‑100 (MQ=255)
                                                                                             ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc      <  1:843710/100‑1 (MQ=255)
                                                                                              gtAATACCACCCGGTCTTTTCTCATTACAGTGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa     >  1:1465050/1‑100 (MQ=255)
                                                                                              gtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa     <  1:2192467/100‑1 (MQ=255)
                                                                                              gtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa     >  2:897657/1‑100 (MQ=255)
                                                                                               tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt    <  1:1111017/100‑1 (MQ=255)
                                                                                                aaTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCCTCCTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATg   >  2:1617107/1‑100 (MQ=255)
                                                                                                aaTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATg   <  2:2378396/100‑1 (MQ=255)
                                                                                                 aTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGa  >  1:179457/1‑100 (MQ=255)
                                                                                                 aTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGa  <  2:139403/100‑1 (MQ=255)
                                                                                                 aTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGa  >  1:2446466/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGA  >  NC_000913/4015236‑4015432

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: