Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,742 A→G 16.4% intergenic (‑107/+22) ygeA ← / ← araE Asp/Glu_racemase family protein/arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,980,7420AG16.4% 417.3 / 66.7 183intergenic (‑107/+22)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (91/60);  new base G (16/14);  total (108/75)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 186 total aligned reads displayed.

TCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980643‑2980841
                                                                                                   |                                                                                                   
tcCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTa                                                                                                     >  2:3679131/1‑100 (MQ=255)
  cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                   >  1:2740467/1‑100 (MQ=255)
    gATGTGTGTTACCGCAACAAGGCTTAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                 <  1:1631262/100‑1 (MQ=255)
      tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                               <  1:572749/100‑1 (MQ=255)
      tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                               >  1:188217/1‑100 (MQ=255)
      tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                               >  2:2262483/1‑100 (MQ=255)
       gtgtgtTACCTCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCACTCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                    <  2:1321858/94‑1 (MQ=255)
       gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                              <  2:3134235/100‑1 (MQ=255)
       gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                              >  1:374300/1‑100 (MQ=255)
          tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           >  1:3076710/1‑100 (MQ=255)
          tgtTACCGCAACAAGGCTGAATAAGGAACGAAACAAACGCCTCAACGGCCAAGTGCCCAATCTAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATcc                                                                                           <  2:587294/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:3552261/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:3668426/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  1:2034639/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                        <  2:1697036/100‑1 (MQ=255)
             tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                        >  1:1526603/1‑100 (MQ=255)
             tACAGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCg                                                                                          <  2:1108679/98‑1 (MQ=255)
             gacagCAACAAGGCTGAATAACGAACGACAAAAAAGCCTCAACGGCCAAGTGCCCAATAAAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATCCGGc                                                                                        <  2:2297895/96‑1 (MQ=255)
              aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCTTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  2:557376/100‑1 (MQ=38)
              aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  1:3124486/100‑1 (MQ=39)
              aCCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  1:326362/100‑1 (MQ=39)
               ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGccc                                                                                      >  2:3050179/1‑100 (MQ=255)
                 gCATCAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAGAAGGGCCGGATGTACAGCACATCCGGcc                                                                                       >  2:3075252/1‑97 (MQ=255)
                 gCAATAATGCTGAATAAGGAACGAGACATACGCCTCAACGGGCAAGTTCCCAATTGCTATCAACGAAAAAAGGGCCGGATGTGGTGTACATCCGGcct                                                                                      <  2:3179424/98‑2 (MQ=255)
                 gCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTAGTAACGAAAAAAGCACCGGATGTACAGCACATccagacc                                                                                      >  2:2671404/1‑93 (MQ=255)
                 gCAACAAGGCTGAATAAGGAACGACACAAATGGCTCGACGGCCAAGTGCCCAATCAAAAAAAAAAAAAAAAGGGCCGGATGTGCTGTACATCCGGccct                                                                                     <  2:652649/99‑2 (MQ=16)
                 gCAACAAGGCTGAAGAAGGAACGACACAAACGCCTCAACGGCCAAGTGCCGAATCAAAAAAAAAGAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                       <  2:1894886/97‑1 (MQ=19)
                  caacaaGGCTGAATAAGGAACGAGACAAACGCCTCCACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   >  2:2668967/1‑100 (MQ=255)
                  caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                   >  1:2932067/1‑100 (MQ=255)
                   aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccctttt                                                                                  <  2:3799219/100‑5 (MQ=38)
                    acaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGcccttttt                                                                                 <  1:1799816/100‑6 (MQ=37)
                       aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                              <  2:3450933/100‑1 (MQ=255)
                        ggCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCa                                                                             >  2:799878/1‑100 (MQ=255)
                         gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                            <  1:3251992/100‑1 (MQ=255)
                         gCTGAATAAGGAACGAGACAAACGCCTCAACAGCCAAGTGCCCAATCTCTATAAACGAAAAAAGGGTCGGATGTACAGCACATCCGGCCCGTGAAAt                                                                               <  2:570341/97‑1 (MQ=255)
                          cTGAAAAAGAAACGAGACAAACGCCGCAACGGCCAATTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACCTCCGGCCCGTGCAATccg                                                                            <  2:3921958/99‑3 (MQ=255)
                            gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACg                                                                         >  1:4180509/1‑100 (MQ=255)
                             aaTAAGGATCTGGACAAACGTCTCAACGGCCAAGAGCCCAATCACTATTAACGAAAAAAGGGCCTGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                            <  2:1473354/96‑1 (MQ=255)
                                  ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCACTATAAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGata                                                                   <  2:2722037/100‑1 (MQ=255)
                                     aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                >  1:990160/1‑100 (MQ=255)
                                      cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                               <  1:2262488/100‑1 (MQ=255)
                                      cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                               <  2:1099800/100‑1 (MQ=255)
                                      cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                               <  1:4023351/100‑1 (MQ=255)
                                      cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                               >  1:1225377/1‑100 (MQ=255)
                                       gagaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtct                                                              >  2:2914748/1‑100 (MQ=255)
                                        agaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtctc                                                             <  1:1637540/100‑1 (MQ=255)
                                           cAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                          >  2:1254323/1‑100 (MQ=255)
                                             aaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                          <  2:1368599/98‑1 (MQ=255)
                                                gCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtctc                                                     >  2:2277272/1‑100 (MQ=255)
                                                 ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                    >  1:693954/1‑100 (MQ=255)
                                                  cTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGc                                                   >  2:106277/1‑100 (MQ=255)
                                                   tCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGcc                                                  >  1:3790134/1‑100 (MQ=255)
                                                    cAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                    <  2:1557933/97‑1 (MQ=255)
                                                        ggCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCa                                             >  2:853476/1‑100 (MQ=255)
                                                         gCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCAt                                            >  1:2586902/1‑100 (MQ=255)
                                                            aaGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAg                                         <  1:4175684/100‑1 (MQ=255)
                                                                gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                     >  1:1073452/1‑100 (MQ=255)
                                                                 cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                    <  2:671775/100‑1 (MQ=255)
                                                                 cccAATCTCTATAAACGAAAAAAGGGCCGGATCTCCCGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                     <  1:3989781/99‑1 (MQ=255)
                                                                  ccAATCTCTATTCACGAAAAAAGGGCCCGATGTACATCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCg                                   <  1:3761836/100‑1 (MQ=255)
                                                                      tctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATACGGCCCGAGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTCCCATCAGTTTGCGtt                                 >  2:1757152/1‑98 (MQ=255)
                                                                        tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                             <  1:801013/100‑1 (MQ=255)
                                                                         ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGata                            >  2:1547038/1‑100 (MQ=255)
                                                                           atgaccaaaaaaaCGGGCGTATGGACCGCACACCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCCCGCCTGCCATCAGTTTGCGTTCGat                             <  2:2553732/90‑1 (MQ=255)
                                                                           aTTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATg                          >  2:75429/1‑100 (MQ=255)
                                                                            ttAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGt                         <  1:3441451/100‑1 (MQ=255)
                                                                             tAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                        >  1:1340772/1‑100 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGCCATCAGACGCCGAGAGTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATg                          >  2:2474218/1‑97 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                        <  1:1262263/99‑1 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                       >  2:3745066/1‑100 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                       >  1:2108003/1‑100 (MQ=255)
                                                                              aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGAGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTGCCATCAGTTTTCGTTGGATATGTTc                       >  2:2582260/1‑100 (MQ=255)
                                                                                cGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGATTGCGTTCGATATGTTCCa                     >  1:872696/1‑100 (MQ=255)
                                                                                 gAAAAAAGGGCCTGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    >  1:2759808/1‑100 (MQ=255)
                                                                                 gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGGTATTTCTCAACATCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    <  2:1770398/100‑1 (MQ=255)
                                                                                 gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    <  1:1097453/100‑1 (MQ=255)
                                                                                 gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                    <  1:1768077/100‑1 (MQ=255)
                                                                                  aaaaaaGGGCCGTATGTACAGCACATCCGTCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGc                   <  1:3755633/100‑1 (MQ=255)
                                                                                   aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                  >  1:1382200/1‑100 (MQ=255)
                                                                                   aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                  >  2:3299560/1‑100 (MQ=255)
                                                                                   aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                  <  1:3402099/100‑1 (MQ=255)
                                                                                    aaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGt                 <  1:1936836/100‑1 (MQ=255)
                                                                                     cacgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                >  2:909699/4‑100 (MQ=38)
                                                                                       aGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc              <  2:1402151/100‑1 (MQ=255)
                                                                                        gtgCCGGATGTACTGCACATCCGGCCCGTGAAATCAGACGCCGATTTGTCTCAACTTCTCGCCTGCCATCAGCTTGCGTTCGATATGTTCCAGCGTg                >  2:3327248/3‑97 (MQ=255)
                                                                                        gtgCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  1:2841982/3‑100 (MQ=255)
                                                                                        gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             <  1:2996665/100‑1 (MQ=255)
                                                                                        gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  2:2565248/1‑100 (MQ=255)
                                                                                         ggCCGGATGTGCTGTACATCCGGCCCGTGACATCAGACGCCGTTATTTCTCAACTTCTCGCCTGCCATCAATCTTCATTCGATAAGTTCCACCGAGACa             >  2:237035/1‑99 (MQ=25)
                                                                                         ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa             >  1:1914142/1‑99 (MQ=39)
                                                                                         ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt            >  2:3876603/1‑100 (MQ=39)
                                                                                          gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGCCGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATGTGTTCCAGGATGACAtt           >  2:508963/1‑100 (MQ=255)
                                                                                          gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  1:2193352/1‑100 (MQ=255)
                                                                                          gCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt           >  2:275794/1‑100 (MQ=255)
                                                                                              gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg       <  2:3262220/100‑1 (MQ=255)
                                                                                                tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt     >  2:4035120/1‑100 (MQ=255)
                                                                                                tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt     <  1:4049309/100‑1 (MQ=255)
                                                                                                  tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGttt   >  1:2398037/1‑100 (MQ=255)
                                                                                                  tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTAGGttt   >  1:2441262/1‑100 (MQ=255)
                                                                                                   aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc  >  1:3162381/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTC  >  NC_000913/2980643‑2980841

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