Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,585 Δ1 bp 100% coding (629/1173 nt) malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,5820G.53.7% 140.9 / 220.8 95coding (626/1173 nt)malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base . (27/24);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 102 total aligned reads displayed.

GAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGG  >  NC_000913/1701484‑1701681
                                                                                                  |                                                                                                   
gAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCAGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgg                                                                                                    >  1:2916415/1‑100 (MQ=255)
  aaGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggg                                                                                                  <  1:1696662/100‑1 (MQ=255)
   aGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcg                                                                                                 <  1:2467336/100‑3 (MQ=255)
    gtgtGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcg                                                                                                 <  2:3464150/99‑3 (MQ=255)
     tgtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcgag                                                                                               >  2:3212821/1‑96 (MQ=255)
     tgtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTgggcgag                                                                                               <  2:2881870/100‑5 (MQ=255)
      gtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGag                                                                                              <  1:721367/100‑1 (MQ=255)
      gtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGag                                                                                              <  2:2883827/100‑1 (MQ=255)
      gtgGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGagc                                                                                             <  1:2543149/100‑1 (MQ=255)
         gACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGagcag                                                                                           >  1:1853022/1‑100 (MQ=255)
         gACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGagcagc                                                                                          <  1:90003/100‑1 (MQ=255)
           cGTGCGATGTGCTGGAGTTCATGGCTGACCTGTGCGAGCGTGAGGGTGTGCGGGTTATTTCCGATGAAATACATGTGGATATGGTTTGGGGCGag                                                                                              <  2:1471504/95‑1 (MQ=255)
           cGTGCGACGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTAGTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCg                                                                                        >  2:2894609/1‑100 (MQ=255)
            gTGCGATGGGCTGGAGATCATGGCTGACCTGTGGGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGc                                                                                       <  2:1307165/100‑1 (MQ=255)
            gTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGc                                                                                       >  1:1743964/1‑100 (MQ=255)
            gTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGc                                                                                       >  2:3174995/1‑100 (MQ=255)
             tGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCCGCCGCa                                                                                      >  2:670175/1‑100 (MQ=255)
             tGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCa                                                                                      <  1:295013/100‑1 (MQ=255)
              gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCa                                                                                      <  1:1191715/100‑1 (MQ=255)
              gCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCa                                                                                      <  2:2158534/100‑1 (MQ=255)
               cGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATGTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCat                                                                                     <  2:1479159/100‑1 (MQ=255)
                 aTGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCatat                                                                                   >  1:3671961/1‑100 (MQ=255)
                 aTGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATAtt                                                                                  >  1:844490/1‑100 (MQ=255)
                   gAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTc                                                                                 <  1:605865/100‑1 (MQ=255)
                    aGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTccc                                                                               <  1:2688828/100‑1 (MQ=255)
                     gCTGGAGATCATGGCTGCCCTGTGCGAGCGTCATGGTTTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTccc                                                                               <  2:3590979/99‑1 (MQ=255)
                       tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTg                                                                             <  2:1888618/100‑1 (MQ=255)
                       tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTgg                                                                            <  1:1575047/100‑1 (MQ=255)
                       tGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTgg                                                                            <  1:3812460/100‑1 (MQ=255)
                           gaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAg                                                                          <  2:829504/99‑1 (MQ=255)
                           gaTCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTa                                                                        >  1:182351/1‑100 (MQ=255)
                               aTGGCTGACCTGTGCGAGCGTCATGGTGGGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAAtgt                                                                    <  2:1528824/100‑1 (MQ=255)
                                tGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAAtgtg                                                                   <  1:2875415/100‑1 (MQ=255)
                                 ggCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                                  >  2:2987796/1‑100 (MQ=255)
                                 ggCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                                  <  2:2996851/100‑1 (MQ=255)
                                  gCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                                  >  1:1523550/1‑100 (MQ=255)
                                  gCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGc                                                                 <  1:2191532/100‑1 (MQ=255)
                                  gCTGACCTGTGCGAGCGTCACGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTgg                                                                  >  2:2523921/1‑100 (MQ=255)
                                    tGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTc                                                               >  1:3160036/1‑100 (MQ=255)
                                      aCCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTcgc                                                             >  1:1687171/1‑100 (MQ=255)
                                       ccTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTcgc                                                             >  1:726066/1‑100 (MQ=255)
                                         tgtgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGga                                                          >  2:3202092/1‑100 (MQ=255)
                                           tgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGga                                                          <  1:761355/99‑1 (MQ=255)
                                           tgCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgaga                                                        <  2:3720946/100‑1 (MQ=255)
                                            gCGAGCGTCATGGTGTGCGGGTTATTTCCGCTGAAATCCATATGGATGTGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGgaga                                                        >  2:14665/1‑100 (MQ=255)
                                            gCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCgg                                                           <  2:2122006/97‑1 (MQ=255)
                                            gCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGAc                                                       <  1:1161726/100‑1 (MQ=255)
                                              gAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACt                                                      <  1:1822397/100‑1 (MQ=255)
                                              gAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTg                                                     >  1:3794337/1‑100 (MQ=255)
                                                  gTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGc                                                  <  1:680318/100‑1 (MQ=255)
                                                  gTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCg                                                 >  1:541347/1‑100 (MQ=255)
                                                   tCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCg                                                 >  1:3887686/1‑100 (MQ=255)
                                                      tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                             >  1:3558051/1‑100 (MQ=255)
                                                      tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                             >  2:1861429/1‑100 (MQ=255)
                                                      tGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                             <  2:3046567/100‑1 (MQ=255)
                                                       ggTGTGCGGGTTATTTCCGATGAAATCCATATGGGTATGGGTTGGGGCGAGCAGCCGCATATTCGCTGGAGTGATGTGGCTCGCGGAGACTGGGCGTTGc                                             >  2:4102180/1‑100 (MQ=255)
                                                       ggTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGc                                             <  2:1275362/100‑1 (MQ=255)
                                                        gtgtGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGGGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCt                                            <  1:2498908/100‑1 (MQ=255)
                                                        gtgtGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTa                                           >  1:2488732/1‑100 (MQ=255)
                                                         tgtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTaa                                          <  1:1923197/100‑1 (MQ=255)
                                                          gtgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGCGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTaa                                          >  2:1343348/1‑100 (MQ=255)
                                                           tgCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACg                                        >  2:148345/1‑100 (MQ=255)
                                                             cGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGt                                       >  1:965680/1‑100 (MQ=255)
                                                              gggTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCg                                     <  1:3087498/100‑1 (MQ=255)
                                                              gggTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCg                                     <  2:1976859/100‑1 (MQ=255)
                                                               ggTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCg                                     >  2:285887/1‑100 (MQ=255)
                                                               ggTTATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCgg                                    >  2:923843/1‑100 (MQ=255)
                                                                 ttATTTCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTAGCGGAGACTGGGCGTTGCTAACGTCGGGc                                  >  2:3687378/1‑100 (MQ=255)
                                                                  tATTTCCGATGAAATCCATATGGCTCTGGTTTGGGGCGCGCGGCCGCATAGTCACCGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGc                                  <  2:2466916/100‑1 (MQ=255)
                                                                   aTTTCCGATGAAATCGATATGGATATGGCTCGGGGCGCGCAGCCACAGAGTCCCTGGAGTAATGTGGCTCCCGGAGACTGGGCGTCTCTAACGTCgagt                                  <  2:813685/99‑4 (MQ=255)
                                                                    tttCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTc                                <  2:587865/100‑1 (MQ=255)
                                                                    tttCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAAAGTc                                      >  2:1783830/1‑93 (MQ=255)
                                                                     ttCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                              <  1:903461/100‑1 (MQ=255)
                                                                     ttCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGa                              <  2:2964308/100‑1 (MQ=255)
                                                                      tCCGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaa                             >  1:1624411/1‑100 (MQ=255)
                                                                       ccGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaaa                            <  2:132920/100‑1 (MQ=255)
                                                                       ccGATGAAATCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGAATGGGCGTTGCTAACGTCGGGCTCGaaa                            >  2:3402854/1‑100 (MQ=255)
                                                                         gATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGaaaa                           <  2:1613761/100‑1 (MQ=255)
                                                                             aaaTCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGttt                       >  2:1522150/1‑100 (MQ=255)
                                                                             aaaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTc                      <  1:1861442/100‑1 (MQ=255)
                                                                              aaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGGCTGGGCGTTGCTGACGTCGGGCTCGAAAAGTTTc                      >  2:3824155/1‑99 (MQ=255)
                                                                              aaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCa                     >  1:2154789/1‑100 (MQ=255)
                                                                              aaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCa                     <  1:724612/100‑1 (MQ=255)
                                                                              aaTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCa                     >  1:560677/1‑100 (MQ=255)
                                                                               aTCCATATGGATATGGTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                    <  1:3694209/100‑1 (MQ=255)
                                                                                    tctggatatggTTTGGGGCGAGCAGCCGCATGTTCCCTGGAGAAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATAtt                >  1:2794784/3‑100 (MQ=255)
                                                                                    tatggatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATAtt                >  2:65821/1‑100 (MQ=255)
                                                                                     atggatatggTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTcc              >  1:3589583/1‑100 (MQ=255)
                                                                                      tggatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTcc              <  1:45979/100‑1 (MQ=255)
                                                                                       ggatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTccc             >  2:3627653/1‑100 (MQ=255)
                                                                                       ggatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGGAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCGCGAAAAGTTTCAATAtttc              >  2:3790595/1‑97 (MQ=255)
                                                                                        gatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCg            <  2:4095563/100‑1 (MQ=255)
                                                                                        gatatggTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTcccac           >  2:2516038/1‑98 (MQ=255)
                                                                                          tatggTTTGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGCCTGGGCGTTGCTAACGTCGGGCTCGACAAGTTTCAATACTCCCGccc         >  2:1965296/1‑100 (MQ=255)
                                                                                           atggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGcc          <  2:865527/99‑1 (MQ=255)
                                                                                            tggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCt        >  1:2358189/1‑100 (MQ=255)
                                                                                            tggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCt        <  1:3274975/100‑1 (MQ=255)
                                                                                             ggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGGAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGcc          <  2:2924459/97‑1 (MQ=255)
                                                                                               tttGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGAc     >  1:3685935/1‑100 (MQ=255)
                                                                                                  ggggCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCgg  >  2:1492790/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GAAAGTGTGGACGTGCGATGAGCTGGAGATCATGGCTGACCTGTGCGAGCGTCATGGTGTGCGGGTTATTTCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGG  >  NC_000913/1701484‑1701681

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: