Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 232,196 G→T 100% G91C (GGC→TGC)  yafE → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913232,1960GT99.1% 375.1 / ‑4.0 112G91C (GGC→TGC) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (61/50);  total (62/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 112 total aligned reads displayed.

CGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACG  >  NC_000913/232099‑232292
                                                                                                 |                                                                                                
cGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGc                                                                                                >  1:2212357/1‑100 (MQ=255)
 gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                               >  1:2194149/1‑100 (MQ=255)
  aTAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCgg                                                                                              <  1:1753274/100‑1 (MQ=255)
    aaCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTa                                                                                            >  1:272226/1‑100 (MQ=255)
     aCGCATTTGATCTTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCgg                                                                                              <  2:2551812/97‑1 (MQ=255)
     aCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAg                                                                                           <  1:1072691/100‑1 (MQ=255)
     aCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAg                                                                                           >  1:1631169/1‑100 (MQ=255)
       gCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGc                                                                                         >  1:2146275/1‑100 (MQ=255)
         aTTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTg                                                                                       >  1:2532839/1‑100 (MQ=255)
           ttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAt                                                                                     >  1:2264530/1‑100 (MQ=255)
           ttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAt                                                                                     <  1:2092264/100‑1 (MQ=255)
            tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                    <  1:648006/100‑1 (MQ=255)
            tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                    >  2:940585/1‑100 (MQ=255)
            tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                    <  1:1427701/100‑1 (MQ=255)
            tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGCAGGCTGAtt                                                                                    >  2:519590/1‑100 (MQ=255)
              atatTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGAGGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                   <  1:2256312/99‑1 (MQ=255)
                atTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGCAGGCTGATTGTGa                                                                                <  1:1739523/100‑1 (MQ=255)
                 ttGTTATCCGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               <  2:2115454/100‑1 (MQ=255)
                 ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               <  2:534445/100‑1 (MQ=255)
                 ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACGGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                               <  1:1830243/100‑1 (MQ=255)
                  tGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATg                                                                              >  1:866745/1‑100 (MQ=255)
                   gTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATgg                                                                             >  2:281762/1‑100 (MQ=255)
                   gTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATgg                                                                             <  1:2292018/100‑1 (MQ=255)
                     tATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGACTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGAc                                                                           >  1:1089488/1‑100 (MQ=255)
                      aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACg                                                                          >  2:2291350/1‑100 (MQ=255)
                      aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACg                                                                          <  1:767186/100‑1 (MQ=255)
                       tCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGt                                                                         >  2:727478/1‑100 (MQ=255)
                        cAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTa                                                                        >  2:2204062/1‑100 (MQ=255)
                         aGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTaa                                                                       >  2:4440/1‑100 (MQ=255)
                         aGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTaa                                                                       >  1:1216066/1‑100 (MQ=255)
                          gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      <  1:1550960/100‑1 (MQ=255)
                          gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      <  1:449802/100‑1 (MQ=255)
                          gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      <  2:800118/100‑1 (MQ=255)
                          gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      <  1:2096268/100‑1 (MQ=255)
                          gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAAt                                                                      >  1:1855949/1‑100 (MQ=255)
                           ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                     >  1:1638569/1‑100 (MQ=255)
                           ccGTTATTCTGCCCATAACAGTCAGGATGTTGGTGCAGCACTGCGAGAAGTGAATAGAATATTGAAACCTTGCGATAGGCTGATTGTGATGGACGGAATg                                                                     <  2:2329084/100‑1 (MQ=255)
                            cGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGt                                                                    >  1:2089785/1‑100 (MQ=255)
                               tATTCTGCCCATCACTAGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTTATTGTGATGGACGTAATGtctc                                                                 <  2:1855060/100‑1 (MQ=255)
                                 ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCg                                                               >  2:2473965/1‑100 (MQ=255)
                                 ttCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAAGATGATTGTGATGGACGTAATGTCTGCg                                                               >  2:2417147/1‑100 (MQ=255)
                                  tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCgg                                                              >  2:76681/1‑100 (MQ=255)
                                   cTGCGCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTCTGATGGACGTAATGTCTCCggg                                                             <  2:1489405/100‑1 (MQ=255)
                                   cTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCggg                                                             <  1:958995/100‑1 (MQ=255)
                                    tGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGt                                                            >  1:1143601/1‑100 (MQ=255)
                                     gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                           <  1:1753416/100‑1 (MQ=255)
                                     gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                           >  2:1576609/1‑100 (MQ=255)
                                     gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                           >  2:1997263/1‑100 (MQ=255)
                                     gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGGAGTGTGATTGTGATGGACGTAATGTCTCCGGGTc                                                           >  2:473270/1‑100 (MQ=255)
                                       ccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAc                                                         <  1:404466/100‑1 (MQ=255)
                                        cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAcc                                                        >  1:471185/1‑100 (MQ=255)
                                         aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAccc                                                       <  1:342799/100‑1 (MQ=255)
                                          tCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCa                                                      >  1:1158361/1‑100 (MQ=255)
                                             cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   >  2:2363558/1‑100 (MQ=255)
                                             cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                   <  2:1091321/100‑1 (MQ=255)
                                                gCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTgcgc                                                >  1:659262/1‑100 (MQ=255)
                                                  atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGa                                              >  1:832454/1‑100 (MQ=255)
                                                  atgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGa                                              <  1:1250133/100‑1 (MQ=255)
                                                   tgatgTTGGTGCAGCGCTGCCAGAAGTGAATAGGATTTTGAATCCTTGCGGTAGGCTGATTGTGATGGACGTAAGGGCTCCTGGTCACCCAg                                                     <  2:1911729/92‑1 (MQ=255)
                                                   tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGCTATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                             <  2:1351781/100‑1 (MQ=255)
                                                    gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCCCCCAGTGCGCGACa                                            <  2:1652986/100‑1 (MQ=255)
                                                      tgatgaagcagcaATGCAATATTTGAATATTTTATTGAAACCTTGCGGTAGGATGATTGTGATGCACGTAATGTATACGGGTCACCCAGTGCCCGACAt                                           >  2:526453/8‑99 (MQ=255)
                                                      tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                          >  2:1510138/1‑100 (MQ=255)
                                                      tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                          >  2:90147/1‑100 (MQ=255)
                                                        ttGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTg                                        <  1:4440/100‑1 (MQ=255)
                                                         tGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                       >  2:190398/1‑100 (MQ=255)
                                                          ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                      <  2:1433165/100‑1 (MQ=255)
                                                           gTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGtt                                     >  2:1165672/1‑100 (MQ=255)
                                                               agcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAg                                 >  1:181048/1‑100 (MQ=255)
                                                                  aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                              <  1:557186/100‑1 (MQ=255)
                                                                  aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                              <  1:602282/100‑1 (MQ=255)
                                                                   cTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGAc                               <  2:2505751/98‑1 (MQ=255)
                                                                    tGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                            <  2:1936056/100‑1 (MQ=255)
                                                                    tGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGt                            >  1:1460110/1‑100 (MQ=255)
                                                                     gCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                           <  1:533973/100‑1 (MQ=255)
                                                                     gCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                           <  2:377486/100‑1 (MQ=255)
                                                                     gCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                           <  1:1165680/100‑1 (MQ=255)
                                                                      cGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                          <  2:2267157/100‑1 (MQ=255)
                                                                         gaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                       <  2:1692144/100‑1 (MQ=255)
                                                                           aGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                     >  1:194986/1‑100 (MQ=255)
                                                                            gTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                    >  2:1112160/1‑100 (MQ=255)
                                                                             tGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                   >  1:1874233/1‑100 (MQ=255)
                                                                             tGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                   >  2:1384705/1‑100 (MQ=255)
                                                                             tGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                   >  1:382149/1‑100 (MQ=255)
                                                                              gAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                  >  2:766353/1‑100 (MQ=255)
                                                                              gAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                  >  1:35817/1‑100 (MQ=255)
                                                                                aTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                <  2:464470/100‑1 (MQ=255)
                                                                                aTAGGATATTGAAACCTTGCGGTAGGATGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                >  1:1555500/1‑100 (MQ=255)
                                                                                 tAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc               <  1:1643163/100‑1 (MQ=255)
                                                                                  aGGATATTGAAGCCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg              >  1:1449439/1‑100 (MQ=255)
                                                                                  aGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg              >  2:7002/1‑100 (MQ=255)
                                                                                   ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa             >  2:60544/1‑100 (MQ=255)
                                                                                   ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa             >  2:727144/1‑100 (MQ=255)
                                                                                     atatTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa           >  2:2145354/1‑100 (MQ=255)
                                                                                      tatTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAc          <  1:1946456/100‑1 (MQ=255)
                                                                                       atTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc         <  1:1058779/100‑1 (MQ=255)
                                                                                        ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGCTAcc         >  1:1039914/1‑99 (MQ=255)
                                                                                        ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct        >  1:2468120/1‑100 (MQ=255)
                                                                                         tGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc       >  1:2280557/1‑100 (MQ=255)
                                                                                              ccTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg  >  1:2530185/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACG  >  NC_000913/232099‑232292

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: