Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 232,196 G→T 69.1% 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,1960GT69.1% 221.5 / 142.4 158G91C (GGC→TGC) yafEputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base T (58/51);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

GCCGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGT  >  NC_000913/232097‑232293
                                                                                                   |                                                                                                 
  cGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGc                                                                                                 <  1:851851/100‑1 (MQ=255)
  cGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGc                                                                                                 <  2:1296680/100‑1 (MQ=255)
   gATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                                >  1:1840108/1‑100 (MQ=255)
     tAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGt                                                                                              >  1:62886/1‑100 (MQ=255)
     tAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCgg                                                                                               <  2:3050986/99‑1 (MQ=255)
      aaCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGt                                                                                              <  1:1651281/99‑1 (MQ=255)
      aaCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTa                                                                                             <  1:2578339/100‑1 (MQ=255)
       aCGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCg                                                                                                >  1:3191300/1‑96 (MQ=255)
        cGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAgg                                                                                           <  2:245736/100‑1 (MQ=255)
         gCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGc                                                                                          >  2:3324202/1‑100 (MQ=255)
          cATTCGATATTTTAATCAGCCGGTCTTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGACAAGTGAATAGGATATTGAAACCTTGCGGTa                                                                                             <  1:2792142/96‑1 (MQ=255)
              tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                     >  1:382024/1‑100 (MQ=255)
              tGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGAtt                                                                                     <  2:1833141/100‑1 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                    >  2:812418/1‑100 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATtg                                                                                    <  2:383755/100‑1 (MQ=255)
               gATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATtg                                                                                    >  1:2131287/1‑100 (MQ=255)
                  atTGTTATCAGCCGTTATTCTGCCCATCACTTGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                 <  2:1554108/100‑1 (MQ=255)
                  atTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGagatcgg                                                                                 >  2:2213351/1‑94 (MQ=255)
                  atTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                 >  2:3205182/1‑100 (MQ=255)
                  atTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGa                                                                                 >  1:2270870/1‑100 (MQ=255)
                  atTGTGATCAGCCGGTATGCTGCCTATCACTGGCATGATGTTGGTGAAGCACTGCGAGAAGTGAATAGGAGATTGAAACCTGGCGGTAGGCTGATTGTGa                                                                                 <  1:2963860/100‑1 (MQ=255)
                   ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                                <  2:1145506/100‑1 (MQ=255)
                   ttGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                                >  1:3380413/1‑100 (MQ=255)
                    tGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGATGCAGCACTGCGCGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGAt                                                                                >  1:265125/1‑99 (MQ=255)
                       tATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGAc                                                                            >  2:3308422/1‑100 (MQ=255)
                         tCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGt                                                                          >  1:422296/1‑100 (MQ=255)
                          cAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTa                                                                         >  1:704418/1‑100 (MQ=255)
                           aGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGAGGTAGGCTGATTGTGATGGACGTa                                                                         >  1:2987846/1‑99 (MQ=255)
                            gCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGGTGGAGGTa                                                                         >  1:643094/1‑98 (MQ=255)
                             ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                      >  2:1484260/1‑100 (MQ=255)
                             ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                      <  2:3249780/100‑1 (MQ=255)
                             ccGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                      >  1:1039495/1‑100 (MQ=255)
                              cGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGGGATGGACGTAATGt                                                                     >  2:1842875/1‑100 (MQ=255)
                              cGTTATTCTGCCCATCACTGGCATGATGATGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATg                                                                      <  1:2043608/99‑1 (MQ=255)
                                  aTTCTGCCCATCACTGGCATGATGTTGGTGCGGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTcc                                                                 >  2:2993444/1‑100 (MQ=255)
                                  aTTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTcc                                                                 >  2:2559431/1‑100 (MQ=255)
                                     cTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCggg                                                              <  2:642229/100‑1 (MQ=255)
                                       gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                            <  2:723928/100‑1 (MQ=255)
                                       gCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                            >  1:980858/1‑100 (MQ=255)
                                        cccATCACTGGCATGATGTTGGTGCAGCCCTGCGAGCAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTc                                                            >  2:3236149/1‑99 (MQ=255)
                                          cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAcc                                                         >  2:1981412/1‑100 (MQ=255)
                                          cATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAcc                                                         >  1:2227489/1‑100 (MQ=255)
                                           aTCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCAccc                                                        >  1:1691852/1‑100 (MQ=255)
                                             cACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAg                                                      <  2:1100275/100‑1 (MQ=255)
                                              aCTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGt                                                     >  1:1011871/1‑100 (MQ=255)
                                               cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGt                                                     <  2:842730/99‑1 (MQ=255)
                                               cTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTg                                                    >  1:935540/1‑100 (MQ=255)
                                                    atgatgTTGGTGCAGCACTGCGAGAAGTGAATAAGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                                <  1:2878301/99‑1 (MQ=255)
                                                     tgatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGAc                                              >  1:2234071/1‑100 (MQ=255)
                                                      gatgTTGGTGCAGCACTGCGAGAAGTTGATAGGATATTGAAAGCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGcgcg                                                <  1:2907651/97‑1 (MQ=255)
                                                      gatgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACa                                             <  2:1537086/100‑1 (MQ=255)
                                                        tgTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATc                                           >  2:3160607/1‑100 (MQ=255)
                                                           tGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                        >  2:3112820/1‑100 (MQ=255)
                                                           tGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTgg                                        <  2:1808556/100‑1 (MQ=255)
                                                            ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                       <  1:1186850/100‑1 (MQ=255)
                                                            ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                       <  1:3353864/100‑1 (MQ=255)
                                                            ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGt                                       <  2:2227504/100‑1 (MQ=255)
                                                            ggTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCAGGTCACCCAGTGCGCGACATCt                                          >  1:2945464/1‑97 (MQ=255)
                                                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                    <  2:2108491/100‑1 (MQ=255)
                                                               gcagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTAc                                    <  2:1730225/100‑1 (MQ=255)
                                                                cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                   <  1:1175912/100‑1 (MQ=255)
                                                                cagcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACa                                   >  2:271416/1‑100 (MQ=255)
                                                                 agcaCTGCGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAg                                  >  2:1242272/1‑100 (MQ=255)
                                                                        cGAGAAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAg                           >  1:1468982/1‑100 (MQ=255)
                                                                           gaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                        <  2:1437611/100‑1 (MQ=255)
                                                                           gaAGTGAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAg                        <  1:1601890/100‑1 (MQ=255)
                                                                            aaGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGc                       >  2:2506213/1‑100 (MQ=255)
                                                                             aGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                      <  1:510548/100‑1 (MQ=255)
                                                                               tGAATAGGATATTGAACCCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCAtt                    <  1:3246372/100‑1 (MQ=255)
                                                                                gAATAGGGTATTGAAACCTGGCGGCAGGCTGATTGTGATGGACGTAATGTCTGCGGGTCAGCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCa                      >  1:1456724/1‑97 (MQ=255)
                                                                                gAATAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTa                   >  1:419940/1‑100 (MQ=255)
                                                                                 aaTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAc                  <  1:756874/100‑1 (MQ=255)
                                                                                  aTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCATCCAGTGCGCGACATCTGGTTACAGGCGGTAGAAGCATTAcg                 >  2:2888327/1‑100 (MQ=255)
                                                                                  aTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                 >  2:2907172/1‑100 (MQ=255)
                                                                                  aTAGGATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcg                 >  1:3068465/1‑100 (MQ=255)
                                                                                   tAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTAcgc                <  1:1871423/100‑1 (MQ=255)
                                                                                     ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa              <  2:2865976/100‑1 (MQ=255)
                                                                                     ggATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa              <  1:3091162/100‑1 (MQ=255)
                                                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa              >  1:1192323/1‑100 (MQ=255)
                                                                                     ggATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGa              >  1:229817/1‑100 (MQ=255)
                                                                                      gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt             <  2:1537795/100‑1 (MQ=255)
                                                                                      gATATTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGAt             <  1:1597807/100‑1 (MQ=255)
                                                                                       atatTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATa            >  1:90267/1‑100 (MQ=255)
                                                                                        tatTGCAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGGCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTGCGCGATAc           >  1:1996880/1‑100 (MQ=255)
                                                                                         atTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc          >  2:1142224/1‑100 (MQ=255)
                                                                                         atTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc          >  1:3273525/1‑100 (MQ=255)
                                                                                         atTGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc          <  2:1011874/100‑1 (MQ=255)
                                                                                         atTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATAcc          <  1:1957665/100‑1 (MQ=255)
                                                                                          ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct         <  2:117378/100‑1 (MQ=255)
                                                                                          ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct         >  2:1344563/1‑100 (MQ=255)
                                                                                          ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct         >  1:2120697/1‑100 (MQ=255)
                                                                                          ttGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACct         >  2:2102996/1‑100 (MQ=255)
                                                                                           tGAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctc        >  1:2602554/1‑100 (MQ=255)
                                                                                            gAAACCTTGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct       >  1:1757929/1‑100 (MQ=255)
                                                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTGACACAGTGCGCGACATCTGGTTACAGACGGTAGACGCATTGCGCGATACct         >  1:401253/1‑98 (MQ=255)
                                                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct       >  2:2347734/1‑100 (MQ=255)
                                                                                            gAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACctct       >  2:932019/1‑100 (MQ=255)
                                                                                              agccTGGCGGTAGGCTGATTGTGATGGACGTAATGTCGCCGGGTCACCCAGTGCGCGACATCTGGTTACGGACGGTAGAAGCATTACGCGATACCTCTCa     <  1:1336446/98‑1 (MQ=255)
                                                                                              aaCCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCa     <  1:115094/100‑1 (MQ=255)
                                                                                                 cTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGt  >  2:1477983/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GCCGATAACGCATTTGATATTGTTATCAGCCGTTATTCTGCCCATCACTGGCATGATGTTGGTGCAGCACTGCGAGAAGTGAATAGGATATTGAAACCTGGCGGTAGGCTGATTGTGATGGACGTAATGTCTCCGGGTCACCCAGTGCGCGACATCTGGTTACAGACGGTAGAAGCATTACGCGATACCTCTCACGT  >  NC_000913/232097‑232293

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: