Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,072,171 C→A 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,1710CA95.5% 302.6 / 3.8 88V161F (GTC→TTC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (30/54);  total (32/56)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCC  >  NC_000913/4072076‑4072266
                                                                                               |                                                                                               
gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                             <  1:2063541/100‑1 (MQ=255)
 tAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCt                                                                                            <  2:1568066/100‑1 (MQ=255)
  agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                           <  1:1107747/100‑1 (MQ=255)
   gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                          <  1:1352255/100‑1 (MQ=255)
    agGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCaa                                                                                         <  2:847829/100‑1 (MQ=255)
     ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                        <  2:2183119/100‑1 (MQ=255)
     ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                        >  2:291333/1‑100 (MQ=255)
     ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                        >  2:569507/1‑100 (MQ=255)
      gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                       <  2:1405275/100‑1 (MQ=255)
      gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTCATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                       <  1:1387596/100‑1 (MQ=255)
        cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAatgat                                                                                     <  2:261565/100‑1 (MQ=255)
        cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                     >  1:174805/1‑100 (MQ=255)
         aGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgatg                                                                                    >  2:1112099/1‑100 (MQ=255)
          gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGCTGAGCGCGTGTTCGCGATCGAATTTATCGCCACGACGCAGTGGGCGGACGACTGGCTCAAtgctgc                                                                                   >  1:56784/1‑100 (MQ=255)
          gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                   <  2:188253/100‑1 (MQ=255)
           cgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCt                                                                                  <  1:1913587/100‑1 (MQ=255)
             ccACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTa                                                                                >  2:1327420/1‑100 (MQ=255)
              cACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAc                                                                               <  1:1112096/100‑1 (MQ=255)
               aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                              <  1:2161557/100‑1 (MQ=255)
               aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                              <  1:574075/100‑1 (MQ=255)
                cTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCa                                                                             <  2:1557155/100‑1 (MQ=255)
                 tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCan                                                                            <  1:553748/100‑2 (MQ=255)
                 tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAg                                                                            <  1:1256719/100‑1 (MQ=255)
                 tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAg                                                                            >  2:246187/1‑100 (MQ=255)
                  gTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGa                                                                           >  1:2032027/1‑100 (MQ=255)
                     gCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCg                                                                        >  1:2593415/1‑100 (MQ=255)
                       ccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtg                                                                      >  2:2081634/1‑100 (MQ=255)
                        cAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtgt                                                                     >  2:2569129/1‑100 (MQ=255)
                         aGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                                    <  1:2211702/100‑1 (MQ=255)
                         aGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGtgtg                                                                    <  1:4706/100‑1 (MQ=255)
                             ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                                >  1:589944/1‑100 (MQ=255)
                              tttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTgg                                                               <  2:2470149/100‑1 (MQ=255)
                               ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                              <  1:258620/100‑1 (MQ=255)
                               ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                              <  2:2438221/100‑1 (MQ=255)
                                tGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGca                                                             >  2:1229989/1‑100 (MQ=255)
                                 ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                            <  2:2593433/100‑1 (MQ=255)
                                 ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                            <  1:2441284/100‑1 (MQ=255)
                                 ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                            <  1:747552/100‑1 (MQ=255)
                                   cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                          >  2:2279979/1‑100 (MQ=255)
                                   cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                          >  2:2218071/1‑100 (MQ=255)
                                         cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                    <  1:2115694/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                   <  1:2533988/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                   <  2:1426116/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                   <  2:2283242/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                   <  1:104607/100‑1 (MQ=255)
                                           atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                  <  1:894801/100‑1 (MQ=255)
                                           atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                  <  1:2447701/100‑1 (MQ=255)
                                           atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                  >  2:115580/1‑100 (MQ=255)
                                           atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                  <  1:682271/100‑1 (MQ=255)
                                             gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                                <  2:1873835/100‑1 (MQ=255)
                                             gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                                <  2:174805/100‑1 (MQ=255)
                                              atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                               >  2:209928/1‑100 (MQ=255)
                                              atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                               <  2:284238/100‑1 (MQ=255)
                                                cGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGCGGCACAGTTCGTTGAActct                                             <  1:1585850/100‑1 (MQ=255)
                                                 gCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCtt                                            >  2:1727560/1‑100 (MQ=255)
                                                 gCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCtt                                            <  2:589945/100‑1 (MQ=255)
                                                  cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCttt                                           <  2:1336375/100‑1 (MQ=255)
                                                  cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCttt                                           >  1:2427901/1‑100 (MQ=255)
                                                    tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAGCAGACCGTGTGAGTGGCACAGTTCTTTGAACTCTTTCa                                         <  1:1605718/100‑1 (MQ=255)
                                                    tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                         <  1:1875007/100‑1 (MQ=255)
                                                    tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                         >  2:226543/1‑100 (MQ=255)
                                                    tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                         >  2:199265/1‑100 (MQ=255)
                                                     gTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAc                                        <  2:2573979/100‑1 (MQ=255)
                                                     gTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAc                                        <  2:657678/100‑1 (MQ=255)
                                                     gTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAc                                        <  1:1327414/100‑1 (MQ=255)
                                                      ttCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAcc                                       <  2:2534865/100‑1 (MQ=255)
                                                         gcgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCGCAGTTCGTTGAACTCTTTCACCata                                    >  1:369256/1‑100 (MQ=255)
                                                           gATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATc                                  >  1:761035/1‑100 (MQ=255)
                                                            aTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATcc                                 >  2:116361/1‑100 (MQ=255)
                                                            aTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATcc                                 <  2:959523/100‑1 (MQ=255)
                                                             tCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCGATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCa                                >  1:1423101/1‑100 (MQ=255)
                                                              cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                               >  2:1067571/1‑100 (MQ=255)
                                                                aaTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGAc                             <  1:1038184/100‑1 (MQ=255)
                                                                 aTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACg                            <  2:1059315/100‑1 (MQ=255)
                                                                       cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                      <  1:1996292/100‑1 (MQ=255)
                                                                       cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                      >  1:705408/1‑100 (MQ=255)
                                                                       cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                      <  1:2569111/100‑1 (MQ=255)
                                                                         cgacgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                    >  2:1758573/3‑100 (MQ=255)
                                                                          cACGACGCGGTGGGCGGACGAATGGCTCAATGAAGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAATTCTTTCACCATATCCAGACGTTGCTgcgc                   >  1:1564654/1‑100 (MQ=255)
                                                                          cACGACGCGGTGGGCGGACGAATGGCCCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcgc                   <  1:1827356/100‑1 (MQ=255)
                                                                           acgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCa                  <  1:2589061/100‑1 (MQ=255)
                                                                               cgcgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCtt              <  1:465535/100‑1 (MQ=255)
                                                                               cgcgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGGGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCtt              <  1:465501/100‑1 (MQ=255)
                                                                                gcgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtc             <  2:1032894/100‑1 (MQ=255)
                                                                                      ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc        <  1:712807/99‑1 (MQ=255)
                                                                                      ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACt       >  2:2379416/1‑100 (MQ=255)
                                                                                       gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTa      >  1:2319370/1‑100 (MQ=255)
                                                                                         ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg    >  2:67801/1‑100 (MQ=255)
                                                                                           acgaaTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGcc  <  1:1229985/100‑1 (MQ=255)
                                                                                               |                                                                                               
GTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCC  >  NC_000913/4072076‑4072266

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