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,1710CA98.8% 516.0 / ‑3.2 153V161F (GTC→TTC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (76/74);  total (79/74)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 158 total aligned reads displayed.

TTGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCCACA  >  NC_000913/4072074‑4072269
                                                                                                 |                                                                                                  
ttGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATg                                                                                                  <  2:460385/100‑1 (MQ=255)
ttGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTGCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATg                                                                                                  <  2:954348/100‑1 (MQ=255)
 tGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATgg                                                                                                 >  2:579106/1‑100 (MQ=255)
  gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                                >  1:2271314/1‑100 (MQ=255)
  gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGCACGAATGGc                                                                                                >  2:1010901/1‑100 (MQ=255)
   tAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGCTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCt                                                                                               >  2:4065577/1‑100 (MQ=255)
    agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                              >  1:1083684/1‑100 (MQ=255)
     gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                             >  1:1673290/1‑100 (MQ=255)
       ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                           >  2:1210498/1‑100 (MQ=255)
       ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                           >  2:3768373/1‑100 (MQ=255)
       ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                           <  2:3445247/100‑1 (MQ=255)
        gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                          >  1:1981523/1‑100 (MQ=255)
        gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                          >  1:2671544/1‑100 (MQ=255)
          cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                        >  1:1613857/1‑100 (MQ=255)
          cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                        >  1:406157/1‑100 (MQ=255)
          cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                        >  1:547940/1‑100 (MQ=255)
          cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                        >  1:2317296/1‑100 (MQ=255)
           aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGTACGAGTGGCTCAatgatg                                                                                       >  1:2716458/1‑100 (MQ=255)
           aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgatg                                                                                       <  1:2854630/100‑1 (MQ=255)
           aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgatg                                                                                       >  1:176871/1‑100 (MQ=255)
            gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                      <  2:1126943/100‑1 (MQ=255)
            gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                      <  1:3838320/100‑1 (MQ=255)
            gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                      <  1:2744448/100‑1 (MQ=255)
            gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                      <  2:2025319/100‑1 (MQ=255)
             cgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCt                                                                                     >  2:340041/1‑100 (MQ=255)
               ccACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTa                                                                                   >  1:43200/1‑100 (MQ=255)
                cACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAc                                                                                  <  2:1292927/100‑1 (MQ=255)
                 aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                                 <  2:3909955/100‑1 (MQ=255)
                 aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGATTa                                                                                   >  1:8847/1‑98 (MQ=255)
                   tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAg                                                                               <  2:2900590/100‑1 (MQ=255)
                   tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAg                                                                               <  2:1880056/100‑1 (MQ=255)
                    gTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGa                                                                              <  1:3768397/100‑1 (MQ=255)
                      cGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGAcc                                                                            >  2:67475/1‑100 (MQ=255)
                        cccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgt                                                                          >  1:2271890/1‑100 (MQ=255)
                         ccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtg                                                                         >  1:2310294/1‑100 (MQ=255)
                          cAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtgt                                                                        >  2:1666790/1‑100 (MQ=255)
                           aGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                                       <  1:2736334/100‑1 (MQ=255)
                            gCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGa                                                                      >  1:3046329/1‑100 (MQ=255)
                             ctctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAg                                                                     <  1:751912/100‑1 (MQ=255)
                             ctctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAg                                                                     >  1:1584192/1‑100 (MQ=255)
                               ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                                   <  1:403634/100‑1 (MQ=255)
                                 ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                                 <  2:176871/100‑1 (MQ=255)
                                 ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                                 <  2:2411673/100‑1 (MQ=255)
                                    gcggcAGCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGAGGCTTAGCAGACCGAGAGAGTGGcaca                                                              >  2:2219274/1‑100 (MQ=255)
                                        cATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                          >  1:3068496/1‑100 (MQ=255)
                                        cATCGATGATCGCTTGATCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGGATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                          >  2:2686433/1‑100 (MQ=255)
                                          tCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGGCCGTGTGAGTGGCACAGTTCGt                                                        <  2:1438299/100‑1 (MQ=255)
                                          tCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGt                                                        >  1:3759947/1‑100 (MQ=255)
                                            gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                      >  2:2300702/1‑100 (MQ=255)
                                            gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                      <  1:3145053/100‑1 (MQ=255)
                                             atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                     <  2:597518/100‑1 (MQ=255)
                                               gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                                   <  1:3408287/100‑1 (MQ=255)
                                               gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                                   <  2:2099462/100‑1 (MQ=255)
                                                atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                                  <  2:3733438/100‑1 (MQ=255)
                                                atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                                  <  2:95766/100‑1 (MQ=255)
                                                  cGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAActct                                                >  1:1232619/1‑100 (MQ=255)
                                                    cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCttt                                              <  1:1381362/100‑1 (MQ=255)
                                                      tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                            <  1:884809/100‑1 (MQ=255)
                                                      tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                            <  2:424694/100‑1 (MQ=255)
                                                       gTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAc                                           <  1:2566621/100‑1 (MQ=255)
                                                         tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAcc                                          >  2:3057097/1‑99 (MQ=255)
                                                         tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                         <  1:3854356/100‑1 (MQ=255)
                                                         tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                         <  2:78447/100‑1 (MQ=255)
                                                          cgcgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCat                                        >  1:394650/1‑100 (MQ=255)
                                                          cgcgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCat                                        <  2:2667799/100‑1 (MQ=255)
                                                          cgcgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCat                                        >  1:35696/1‑100 (MQ=255)
                                                            cgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCatat                                      >  1:2001267/1‑100 (MQ=255)
                                                              aTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATcc                                    >  1:3156591/1‑100 (MQ=255)
                                                                  aaTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGAc                                <  2:969387/100‑1 (MQ=255)
                                                                  aaTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGAc                                >  1:1063268/1‑100 (MQ=255)
                                                                      tATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtg                            >  1:994825/1‑100 (MQ=255)
                                                                       aTCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtg                            >  2:3919963/1‑99 (MQ=255)
                                                                       aTCACCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTTAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgc                           >  1:801272/1‑100 (MQ=255)
                                                                          gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                        >  1:3956095/1‑100 (MQ=255)
                                                                          gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                        >  1:474754/1‑100 (MQ=255)
                                                                           ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                       <  1:2383324/100‑1 (MQ=255)
                                                                            cACGACGCGGTGGGCGGACGAATGGCTTAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                       <  2:3539586/99‑1 (MQ=255)
                                                                            cACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcgc                      >  1:3928585/1‑100 (MQ=255)
                                                                             acgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCa                     >  1:441772/1‑100 (MQ=255)
                                                                                acgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCt                  <  2:2449532/100‑1 (MQ=255)
                                                                                acgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGGTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCt                  <  2:3883107/100‑1 (MQ=255)
                                                                                    ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat              >  1:723692/1‑100 (MQ=255)
                                                                                    ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat              >  2:225917/1‑100 (MQ=255)
                                                                                    ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat              >  1:1865798/1‑100 (MQ=255)
                                                                                     gTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatc             <  2:1083684/100‑1 (MQ=255)
                                                                                      tGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatca            <  1:1841665/100‑1 (MQ=255)
                                                                                      tGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatca            >  2:4113069/1‑100 (MQ=255)
                                                                                       gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCCTCAc           >  1:417667/1‑100 (MQ=255)
                                                                                       gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc           >  1:2775215/1‑100 (MQ=255)
                                                                                       gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc           >  2:3544179/1‑100 (MQ=255)
                                                                                        ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACt          >  1:3347012/1‑100 (MQ=255)
                                                                                         gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatcn            <  2:795891/97‑2 (MQ=255)
                                                                                         gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTa         >  1:1051947/1‑100 (MQ=255)
                                                                                          cGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTAc        >  2:1704125/1‑100 (MQ=255)
                                                                                           ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg       <  2:2519432/100‑1 (MQ=255)
                                                                                           ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg       >  1:2103134/1‑100 (MQ=255)
                                                                                           ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg       >  2:1654353/1‑100 (MQ=255)
                                                                                             acgaaTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGcc     <  2:430296/100‑1 (MQ=255)
                                                                                             acgaaTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGcc     >  2:3961577/1‑100 (MQ=255)
                                                                                                acTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCcaca  <  2:645118/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TTGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCCACA  >  NC_000913/4072074‑4072269

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