Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,744 A→T 17.0% intergenic (‑109/+20) 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,7440AT17.0% 442.8 / 72.8 182intergenic (‑109/+20)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (92/57);  new base T (17/14);  total (109/73)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 186 total aligned reads displayed.

CTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCG  >  NC_000913/2980645‑2980843
                                                                                                   |                                                                                                   
cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                     >  1:2740467/1‑100 (MQ=255)
   aTGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGca                                                                                                  <  1:2333840/100‑1 (MQ=255)
    tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcac                                                                                                 >  2:2262483/1‑100 (MQ=255)
     gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                                >  1:374300/1‑100 (MQ=255)
     gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                                <  2:3134235/100‑1 (MQ=255)
        tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                             >  1:3076710/1‑100 (MQ=255)
        tgtTACCGCAACAAGGCTGAATAAGGAACGAAACAAACGCCTCAACGGCCAAGTGCCCAATCTAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATcc                                                                                             <  2:587294/100‑1 (MQ=255)
         gtTACCGCAACAAGGCTGAATAAGGAACGATATAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAATTGGCCGGATGTACAGCACATCCg                                                                                            <  2:2746080/100‑1 (MQ=255)
           tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  2:1697036/100‑1 (MQ=255)
           tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  1:3552261/100‑1 (MQ=255)
           tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                          <  1:3668426/100‑1 (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                                                                                        >  1:1468510/1‑100 (MQ=255)
             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)
                caacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTg                                                                                     >  1:2932067/1‑100 (MQ=255)
                caacaaGGCTGAATAAGGAACGAGACAAACGCATCAACGGCGAAGTGCCCAATCACCATTAACGAAAAGACGGCCGAGGGTACAGCCCATccagagcg                                                                                       >  2:2638621/1‑92 (MQ=255)
                 aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGccctttt                                                                                    <  2:3799219/100‑5 (MQ=38)
                  acaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGaa                                                                                   >  2:2785287/1‑100 (MQ=255)
                     aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                <  2:3450933/100‑1 (MQ=255)
                      ggCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCa                                                                               >  2:799878/1‑100 (MQ=255)
                       gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                              <  1:3251992/100‑1 (MQ=255)
                       gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                              <  2:27043/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)
                                   aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                  >  1:990160/1‑100 (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)
                                          aaaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACt                                                           >  1:240749/1‑100 (MQ=255)
                                             cGCCTCAAGGGCCAAGTGCCCAATCTCTATGAACGAAAAAAGGGCCGGATGGACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtct                                                        <  2:139084/100‑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)
                                                   aaCGGCCAAGTGCCCCATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTg                                                  <  1:2626123/100‑1 (MQ=255)
                                                    aCGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGc                                                 <  1:1918291/100‑1 (MQ=255)
                                                      ggCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCa                                               >  2:853476/1‑100 (MQ=255)
                                                          aaGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAg                                           <  1:4175684/100‑1 (MQ=255)
                                                           aGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAg                                           >  2:1465255/1‑99 (MQ=255)
                                                              gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                       >  1:1073452/1‑100 (MQ=255)
                                                              gCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTg                                       >  1:1962227/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)
                                                                 cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                      >  1:1802185/1‑98 (MQ=255)
                                                                    tctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATACGGCCCGAGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTCCCATCAGTTTGCGtt                                   >  2:1757152/1‑98 (MQ=255)
                                                                      tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                               <  1:801013/100‑1 (MQ=255)
                                                                       ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGata                              <  1:1902612/100‑1 (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:1605888/1‑100 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGCCATCAGACGCCGAGAGTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATg                            >  2:2474218/1‑97 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                          <  1:1262263/99‑1 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                         >  1:2108003/1‑100 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                         >  2:3745066/1‑100 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                         >  1:1715378/1‑100 (MQ=255)
                                                                            aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                         >  1:1362669/1‑100 (MQ=255)
                                                                             aCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTcc                        >  1:1770150/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                      >  2:1511160/1‑100 (MQ=255)
                                                                                aaaaaaGGGCCGTATGTACAGCACATCCGTCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGc                     <  1:3755633/100‑1 (MQ=255)
                                                                                 aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                    >  2:3299560/1‑100 (MQ=255)
                                                                                 aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                    <  1:3402099/100‑1 (MQ=255)
                                                                                   cacgggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTg                  >  2:909699/4‑100 (MQ=38)
                                                                                      gtgCCGGATGTACTGCACATCCGGCCCGTGAAATCAGACGCCGATTTGTCTCAACTTCTCGCCTGCCATCAGCTTGCGTTCGATATGTTCCAGCGTg                  >  2:3327248/3‑97 (MQ=255)
                                                                                      gtgCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa               >  1:2841982/3‑100 (MQ=255)
                                                                                      gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa               >  2:2565248/1‑100 (MQ=255)
                                                                                      gggCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACa               <  1:2996665/100‑1 (MQ=255)
                                                                                       ggCCGGATGTGCTGTACATCCGGCCCGTGACATCAGACGCCGTTATTTCTCAACTTCTCGCCTGCCATCAATCTTCATTCGATAAGTTCCACCGAGACa               >  2:237035/1‑99 (MQ=25)
                                                                                       ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt              >  2:3876603/1‑100 (MQ=39)
                                                                                        gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGCCGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATGTGTTCCAGGATGACAtt             >  2:508963/1‑100 (MQ=255)
                                                                                        gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt             >  2:135706/1‑100 (MQ=255)
                                                                                        gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt             >  1:2246637/1‑100 (MQ=255)
                                                                                            gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg         <  2:3262220/100‑1 (MQ=255)
                                                                                              tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt       <  1:4049309/100‑1 (MQ=255)
                                                                                              tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt       >  2:4035120/1‑100 (MQ=255)
                                                                                               gTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGtt      <  1:2122380/100‑1 (MQ=255)
                                                                                                tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGttt     >  1:2538300/1‑100 (MQ=255)
                                                                                                 aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTc    >  1:3162381/1‑100 (MQ=255)
                                                                                                  cAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc   >  1:51867/1‑100 (MQ=255)
                                                                                                  cAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc   >  2:3714241/1‑100 (MQ=255)
                                                                                                   aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc   >  2:1604571/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
CTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCG  >  NC_000913/2980645‑2980843

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: