Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 232,196 G→T 24.6% 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,1960GT24.6% 371.0 / 122.3 188G91C (GGC→TGC) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (70/71);  new base T (18/28);  total (88/100)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 193 total aligned reads displayed.

GCCGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAC  >  NC_000913/232097‑232295
                                                                                                   |                                                                                                   
gCCGATAACGCATTTGATATTGTTACCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAACAGGATATTGAAACCTg                                                                                                     >  1:2416350/1‑100 (MQ=255)
 ccGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTgg                                                                                                    >  1:3312844/1‑100 (MQ=255)
  cGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGc                                                                                                   <  2:3948404/100‑1 (MQ=255)
   gATGACGCATTTGCTACTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                                  <  2:2859108/100‑1 (MQ=255)
   gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                                  <  1:3911834/100‑1 (MQ=255)
   gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                                  <  1:2157244/100‑1 (MQ=255)
   gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCg                                                                                                  <  1:2113648/100‑1 (MQ=255)
    aTAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAATCCTGGCg                                                                                                  <  1:1226312/99‑1 (MQ=255)
    aTAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCgg                                                                                                 <  2:1577748/100‑1 (MQ=255)
      aaCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTa                                                                                               >  1:3738887/1‑100 (MQ=255)
       aCGCATTTGATATTGTTATCAGCCGTTATTCTGCTCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAg                                                                                              >  1:1783405/1‑100 (MQ=255)
       aCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAg                                                                                              >  2:2264829/1‑100 (MQ=255)
        cGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAgg                                                                                             >  1:1713379/1‑100 (MQ=255)
         gCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGc                                                                                            <  2:525597/100‑1 (MQ=255)
         gCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGc                                                                                            <  1:3459731/100‑1 (MQ=255)
         gCATTGGATATTGTGATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGACGGTAGGc                                                                                            <  2:3361789/100‑1 (MQ=255)
           aTTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTg                                                                                          >  1:3919759/1‑100 (MQ=255)
             ttGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGAt                                                                                        >  2:2335871/1‑100 (MQ=255)
              tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGGAGGCTGAtt                                                                                       >  2:3479034/1‑100 (MQ=255)
              tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGAtt                                                                                       >  2:1062628/1‑100 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                      <  1:1005332/100‑1 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                      >  2:3369645/1‑100 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATtg                                                                                      <  1:343156/100‑1 (MQ=255)
                   ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                                  <  2:2401874/100‑1 (MQ=255)
                   ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGAt                                                                                  <  2:2745866/100‑1 (MQ=255)
                    tGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGa                                                                                   >  2:631001/1‑98 (MQ=255)
                     gTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATgg                                                                                <  1:706299/100‑1 (MQ=255)
                       tATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGAc                                                                              <  1:1004503/100‑1 (MQ=255)
                       tATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGAc                                                                              <  2:1995788/100‑1 (MQ=255)
                       tATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGAc                                                                              <  1:3420613/100‑1 (MQ=255)
                        aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACg                                                                             <  2:3276824/100‑1 (MQ=255)
                        aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGAc                                                                              >  2:3715389/1‑99 (MQ=255)
                        aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACg                                                                             <  2:2122821/100‑1 (MQ=255)
                        aTCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACg                                                                             <  2:2103253/100‑1 (MQ=255)
                           aGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTaa                                                                          >  1:2473818/1‑100 (MQ=255)
                            gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAAt                                                                         >  2:1439346/1‑100 (MQ=255)
                              cGTTATTCTGCCCATCACTTGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATg                                                                        <  2:968291/99‑1 (MQ=255)
                              cGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGt                                                                       >  2:1656602/1‑100 (MQ=255)
                                    tCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCgg                                                                 <  2:3309507/100‑1 (MQ=255)
                                     cTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCggg                                                                <  1:1161463/100‑1 (MQ=255)
                                        cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCa                                                             <  1:3001761/100‑1 (MQ=255)
                                        cccATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCa                                                             >  1:2193321/1‑100 (MQ=255)
                                             cACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAg                                                        <  2:3663040/100‑1 (MQ=255)
                                             cACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAg                                                        >  1:266026/1‑100 (MQ=255)
                                                 ggCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTgcg                                                    <  2:3217181/100‑1 (MQ=255)
                                                   cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                                  <  1:2903713/100‑1 (MQ=255)
                                                   cATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                                  <  1:971262/100‑1 (MQ=255)
                                                      gatgTTGGTGCAGCACTGCGAGAAGTTAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                               <  2:2378209/100‑1 (MQ=255)
                                                      gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                               >  1:2534915/1‑100 (MQ=255)
                                                           tGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                          >  2:750272/1‑100 (MQ=255)
                                                            ggTGCAGCACTGGGAGAAGTGAATAGGATATTGAAGCCTGGCGGTAGGCTGAGTGGGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                         <  2:2967197/100‑1 (MQ=255)
                                                             gTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGtt                                        <  2:1205982/100‑1 (MQ=255)
                                                              tGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTa                                       <  1:3558411/100‑1 (MQ=255)
                                                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                      >  2:2516684/1‑100 (MQ=255)
                                                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                      >  1:310514/1‑100 (MQ=255)
                                                                cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                     >  1:1378871/1‑100 (MQ=255)
                                                                  gcaCTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGa                                   >  1:3981795/1‑100 (MQ=255)
                                                                    aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                 <  1:2945118/100‑1 (MQ=255)
                                                                    aCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACg                                 <  1:3377046/100‑1 (MQ=255)
                                                                     cTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACgg                                <  1:3973891/100‑1 (MQ=255)
                                                                      tGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCAAGTGCGCGACATCTGGTTACAGACGGt                               <  1:3962966/100‑1 (MQ=255)
                                                                       ncGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGANGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                              >  1:2023735/2‑100 (MQ=255)
                                                                       gCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTa                              <  1:1542930/100‑1 (MQ=255)
                                                                        cGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGTCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                             <  2:1782294/100‑1 (MQ=255)
                                                                          agaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGGAGaa                           <  1:7350/100‑1 (MQ=255)
                                                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                          <  1:3476039/100‑1 (MQ=255)
                                                                           gaAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                          >  1:3185869/1‑100 (MQ=255)
                                                                             aGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                        <  2:1120329/100‑1 (MQ=255)
                                                                             aGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                        >  1:3331474/1‑100 (MQ=255)
                                                                              gTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                       >  1:1422066/1‑100 (MQ=255)
                                                                              gTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAt                       >  1:1826083/1‑100 (MQ=255)
                                                                               tGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                      >  1:3306471/1‑100 (MQ=255)
                                                                                gaacaGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                     <  2:2416336/96‑1 (MQ=255)
                                                                                gAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                     >  1:531858/1‑100 (MQ=255)
                                                                                   tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                  <  1:801587/100‑1 (MQ=255)
                                                                                   tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                  >  1:1911479/1‑100 (MQ=255)
                                                                                    aGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg                 <  2:1263601/100‑1 (MQ=255)
                                                                                    aGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg                 <  1:3222047/100‑1 (MQ=255)
                                                                                    aGGATATTGAAACCATGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgcg                 <  2:2811229/100‑1 (MQ=255)
                                                                                     ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                >  2:3879655/1‑100 (MQ=255)
                                                                                     ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                >  1:1869962/1‑100 (MQ=255)
                                                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa                >  1:3932193/1‑100 (MQ=255)
                                                                                      gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt               <  2:3031406/100‑1 (MQ=255)
                                                                                         atTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc            <  1:1648169/100‑1 (MQ=255)
                                                                                         atTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc            <  1:359093/100‑1 (MQ=255)
                                                                                         atTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc            >  1:2062831/1‑100 (MQ=255)
                                                                                          ttGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct           >  1:491255/1‑100 (MQ=255)
                                                                                           tGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc          >  2:1967491/1‑100 (MQ=255)
                                                                                           tGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc          <  1:838704/100‑1 (MQ=255)
                                                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct         >  2:3807533/1‑100 (MQ=255)
                                                                                             aatccTGGCGGTAGGCTGATTGTGATGGACGTAGTGTCTCCGGGTCACCCAGCGCGCGACATCTGGTTACAGCCGGTAGAAGCATTACGCGATACCtctc        <  1:4106056/97‑1 (MQ=255)
                                                                                             aaaCCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCtctc        <  2:14781/100‑1 (MQ=255)
                                                                                              atccTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATATGGTTACAGACGGTAGAAGCATTACGCGATACCtctc        <  2:2215516/97‑1 (MQ=255)
                                                                                              aaCCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa       <  2:2562410/100‑1 (MQ=255)
                                                                                              aaCCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa       <  2:190204/100‑1 (MQ=255)
                                                                                               aCCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCAc      <  2:3885028/100‑1 (MQ=255)
                                                                                                ccTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTGCAGACGGTAGAAGCATTACGCGATACCTCTCAc      <  2:2473804/99‑1 (MQ=255)
                                                                                                ccTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg     >  2:895999/1‑100 (MQ=255)
                                                                                                ccTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACg     >  1:3380364/1‑100 (MQ=255)
                                                                                                 cTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGt    >  2:3376652/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCCGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGTAC  >  NC_000913/232097‑232295

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