Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,171 C→A 100% V161F (GTC→TTC)  yihT ← 6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,1710CA99.4% 569.3 / ‑4.2 167V161F (GTC→TTC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (78/86);  total (79/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CTTTGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAC  >  NC_000913/4072072‑4072260
                                                                                                   |                                                                                         
cTTTGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGCTCGCTTGTTCGCGATCGAATTTATCGACACGACGCGGTGGGCGgacgac                                                                                           >  2:327959/1‑100 (MQ=255)
  ttGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATg                                                                                         <  1:4014669/100‑1 (MQ=255)
    gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                       >  1:1136095/1‑100 (MQ=255)
    gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                       <  2:3067059/100‑1 (MQ=255)
    gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                       <  2:3712417/100‑1 (MQ=255)
     tAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCt                                                                                      >  1:1074067/1‑100 (MQ=255)
      agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                     >  2:748191/1‑100 (MQ=255)
       gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                    >  1:20091/1‑100 (MQ=255)
       gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                    >  1:3888259/1‑100 (MQ=255)
       gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                    >  2:4311659/1‑100 (MQ=255)
        agGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCaa                                                                                   <  2:3432777/100‑1 (MQ=255)
        agGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCaa                                                                                   <  1:1343584/100‑1 (MQ=255)
         ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                  >  1:3288/1‑100 (MQ=255)
         ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                  <  2:2107801/100‑1 (MQ=255)
          gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                 <  2:3237273/100‑1 (MQ=255)
          gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                 <  1:995871/100‑1 (MQ=255)
            cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                               >  2:2590046/1‑100 (MQ=255)
            cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                               >  2:2027939/1‑100 (MQ=255)
             aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGGGGTGGGCGGACGAATGGCTCAatgatg                                                                              >  1:1593452/1‑100 (MQ=255)
             aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgatg                                                                              <  2:2358379/100‑1 (MQ=255)
              gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgatg                                                                              <  1:4166316/99‑1 (MQ=255)
              gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                             <  2:3726647/100‑1 (MQ=255)
                gcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCtt                                                                           <  2:1344066/100‑1 (MQ=255)
                gcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCtt                                                                           >  2:77289/1‑100 (MQ=255)
                gcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCtt                                                                           >  1:3118158/1‑100 (MQ=255)
                gcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCtt                                                                           <  1:2990188/100‑1 (MQ=255)
                 ccACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCTCTTTCGATTTTTTCGCCACGACGCTGTGTTCTTTTGATTGGCTCAatgatgt                                                                             <  1:1612929/97‑2 (MQ=255)
                 ccACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTa                                                                          >  1:3273489/1‑100 (MQ=255)
                  cACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAc                                                                         >  2:3050567/1‑100 (MQ=255)
                   aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                        <  2:1055873/100‑1 (MQ=255)
                      gTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGa                                                                     >  2:1554027/1‑100 (MQ=255)
                         gCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCg                                                                  >  2:2436121/1‑100 (MQ=255)
                         gCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCg                                                                  >  1:3481106/1‑100 (MQ=255)
                         gCCCAGCTCTTTGGCGGCATAGGTGATCGCTTGTTCGCGATCGACTTTATCGCCACGACGCGTCGGGCGGACGAATGGCTCAATGTGGGTGACCAGACCg                                                                  >  1:715762/1‑100 (MQ=255)
                             agttctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                              <  1:1800257/97‑1 (MQ=255)
                             aggtctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTAGCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                              <  1:357416/97‑1 (MQ=255)
                             aGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                              <  2:176657/100‑1 (MQ=255)
                             aGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                              <  1:2295145/100‑1 (MQ=255)
                              gCTCTTTGGCGGCATCCATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCg                                                                  >  1:1394827/1‑95 (MQ=255)
                                 ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                          <  2:4043313/100‑1 (MQ=255)
                                  tttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTgg                                                         <  1:122362/100‑1 (MQ=255)
                                   ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                        <  2:1243645/100‑1 (MQ=255)
                                   ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                        <  2:118286/100‑1 (MQ=255)
                                   ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                        <  1:291182/100‑1 (MQ=255)
                                     ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                      >  2:3199092/1‑100 (MQ=255)
                                     ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                      <  1:1656223/100‑1 (MQ=255)
                                     ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                      >  2:555907/1‑100 (MQ=255)
                                     ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                      <  2:3410840/100‑1 (MQ=255)
                                      gcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcaca                                                     >  2:1184392/1‑100 (MQ=255)
                                        ggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGt                                                   >  1:1633739/1‑100 (MQ=255)
                                          cATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                 >  2:3971021/1‑100 (MQ=255)
                                          cATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                 >  2:1613603/1‑100 (MQ=255)
                                           aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                >  2:2712231/1‑100 (MQ=255)
                                           aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                >  1:4222544/1‑100 (MQ=255)
                                           aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                <  2:1482889/100‑1 (MQ=255)
                                           aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                >  2:1145584/1‑100 (MQ=255)
                                              gatgatCGCTTGTTCGCGATCGAATTTATCGCCTCGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                             <  1:3713359/100‑1 (MQ=255)
                                              gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                             >  1:2790853/1‑100 (MQ=255)
                                               atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                            <  2:3850918/100‑1 (MQ=255)
                                                tgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGaa                                           <  2:3037952/100‑1 (MQ=255)
                                                 gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                          <  2:2749413/100‑1 (MQ=255)
                                                      cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCtt                                      >  2:2236804/1‑99 (MQ=255)
                                                      cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGGGAGTGGCACAGTTCGTTGAACTCttt                                     >  1:2803660/1‑100 (MQ=255)
                                                       ttGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTc                                    <  1:2590033/100‑1 (MQ=255)
                                                        tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                   <  2:1092606/100‑1 (MQ=255)
                                                           tagcgATCGAATTAATCGCCACGACGCGGTGCGCGGACGAATAGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                >  2:2970115/3‑100 (MQ=255)
                                                           tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                >  2:1048135/1‑100 (MQ=255)
                                                                aTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATcc                           <  2:2727979/100‑1 (MQ=255)
                                                                aTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATcc                           >  1:4220258/1‑100 (MQ=255)
                                                                 tCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCa                          <  1:2838632/100‑1 (MQ=255)
                                                                  cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                         <  1:217906/100‑1 (MQ=255)
                                                                  cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                         >  2:3789114/1‑100 (MQ=255)
                                                                  cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                         <  2:1118951/100‑1 (MQ=255)
                                                                   gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCa                          >  1:3954820/1‑98 (MQ=255)
                                                                   gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGa                        <  2:3266378/100‑1 (MQ=255)
                                                                    aaTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGAc                       <  1:568118/100‑1 (MQ=255)
                                                                        tATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtg                   >  2:1581774/1‑100 (MQ=255)
                                                                        tATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtg                   >  2:626834/1‑100 (MQ=255)
                                                                         aTCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgc                  <  1:3626688/100‑1 (MQ=255)
                                                                         aTCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgc                  >  2:1682670/1‑100 (MQ=255)
                                                                         aTCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgc                  <  1:442663/100‑1 (MQ=255)
                                                                            gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCCCCATATCAACACGTtgctga               >  2:1903581/1‑99 (MQ=255)
                                                                            gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc               >  1:720197/1‑100 (MQ=255)
                                                                            gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                <  1:2447516/99‑1 (MQ=255)
                                                                             ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg              >  1:1711229/1‑100 (MQ=255)
                                                                              cACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcgc             <  2:3825260/100‑1 (MQ=255)
                                                                              cACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcgc             >  2:1033863/1‑100 (MQ=255)
                                                                                cgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATAGCCAGACGTTGCTGCGCAt           >  1:1839888/1‑100 (MQ=255)
                                                                                 gacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATc          <  2:1281499/100‑1 (MQ=255)
                                                                                  acgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCt         <  1:925494/100‑1 (MQ=255)
                                                                                  acgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCt         >  2:4016686/1‑100 (MQ=255)
                                                                                     cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca      <  1:62780/100‑1 (MQ=255)
                                                                                     cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca      <  1:531090/100‑1 (MQ=255)
                                                                                     cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca      <  1:1243324/100‑1 (MQ=255)
                                                                                     cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca      <  1:1844531/100‑1 (MQ=255)
                                                                                      ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat     >  1:2069493/1‑100 (MQ=255)
                                                                                      ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat     >  1:1769935/1‑100 (MQ=255)
                                                                                         gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc  >  1:2700100/1‑100 (MQ=255)
                                                                                          ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc  <  1:1524543/99‑1 (MQ=255)
                                                                                          ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc  <  1:1025094/99‑1 (MQ=255)
                                                                                                   |                                                                                         
CTTTGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAC  >  NC_000913/4072072‑4072260

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