Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,746 C→T 16.7% intergenic (‑111/+18) 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,7460CT16.7% 431.8 / 64.4 178intergenic (‑111/+18)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base C (93/55);  new base T (16/14);  total (109/69)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-02

Only 100 of 188 total aligned reads displayed.

GATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGG  >  NC_000913/2980647‑2980844
                                                                                                   |                                                                                                  
gATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGc                                                                                                    <  2:2463847/100‑1 (MQ=255)
 aTGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGca                                                                                                   <  1:2333840/100‑1 (MQ=255)
   gtgtgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGcaca                                                                                                 >  1:374300/1‑100 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                              >  1:3076710/1‑100 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                              <  2:569740/100‑1 (MQ=255)
      tgtTACCGCAACAAGGCTGAATAAGGAACGAAACAAACGCCTCAACGGCCAAGTGCCCAATCTAAAAAAAAAAAAAAAGGGCCGGATGTACAGCACATcc                                                                                              <  2:587294/100‑1 (MQ=255)
       gtTACCGCAACAAGGCTGAATAAGGAACGATATAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAATTGGCCGGATGTACAGCACATCCg                                                                                             <  2:2746080/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGTCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                           <  1:3996755/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                           <  1:3668426/100‑1 (MQ=255)
         tACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCGGc                                                                                           <  2:1697036/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:326362/100‑1 (MQ=39)
           ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGccc                                                                                         >  2:3050179/1‑100 (MQ=255)
           ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGccc                                                                                         >  1:1468510/1‑100 (MQ=255)
           ccGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATAACAATAAAAAAAAAAAGGGCCGGATGTGCTGTACATCCGGcc                                                                                          <  2:3663842/99‑1 (MQ=25)
             gCATCAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAGAAGGGCCGGATGTACAGCACATCCGGcc                                                                                          >  2:3075252/1‑97 (MQ=255)
             gCAATAATGCTGAATAAGGAACGAGACATACGCCTCAACGGGCAAGTTCCCAATTGCTATCAACGAAAAAAGGGCCGGATGTGGTGTACATCCGGcct                                                                                         <  2:3179424/98‑2 (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)
               aacaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGa                                                                                     >  1:376766/1‑100 (MQ=255)
                acaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGaa                                                                                    >  2:2785287/1‑100 (MQ=255)
                 caaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGaaa                                                                                   <  1:445900/100‑1 (MQ=255)
                   aGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATc                                                                                 <  2:3450933/100‑1 (MQ=255)
                    ggCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCa                                                                                >  2:799878/1‑100 (MQ=255)
                    ggCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCa                                                                                >  1:3919381/1‑100 (MQ=255)
                     gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                               <  1:3251992/100‑1 (MQ=255)
                     gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAg                                                                               <  2:27043/100‑1 (MQ=255)
                      cTGAAAAAGAAACGAGACAAACGCCGCAACGGCCAATTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACCTCCGGCCCGTGCAATccg                                                                               <  2:3921958/99‑3 (MQ=255)
                          aTAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGcc                                                                          <  1:3124450/100‑1 (MQ=255)
                                 aCGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATAttt                                                                   >  1:990160/1‑100 (MQ=255)
                                  cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                                  <  2:1099800/100‑1 (MQ=255)
                                  cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTtc                                                                  <  1:2262488/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)
                                     gaCAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCa                                                               >  1:408195/1‑100 (MQ=255)
                                       cAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAAc                                                             >  2:1254323/1‑100 (MQ=255)
                                        aaaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACt                                                            >  1:240749/1‑100 (MQ=255)
                                           cGCCTCAAGGGCCAAGTGCCCAATCTCTATGAACGAAAAAAGGGCCGGATGGACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTtct                                                         <  2:139084/100‑1 (MQ=255)
                                             ccTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                       >  1:693954/1‑100 (MQ=255)
                                              cTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCg                                                       >  2:1662140/1‑99 (MQ=255)
                                              cTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGc                                                      >  2:106277/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)
                                                        aaGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAg                                            <  1:3757647/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)
                                                             cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                       >  1:3153051/1‑100 (MQ=255)
                                                               cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGc                                       >  1:1802185/1‑98 (MQ=255)
                                                                  tctcTATTAACGAAAAAAGGGCCGGATGTACAGCACATACGGCCCGAGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTCCCATCAGTTTGCGtt                                    >  2:1757152/1‑98 (MQ=255)
                                                                    tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                                <  2:3739894/100‑1 (MQ=255)
                                                                    tctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGat                                <  1:801013/100‑1 (MQ=255)
                                                                     ctATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGata                               <  1:1902612/100‑1 (MQ=255)
                                                                      tATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGatat                              >  1:3439104/1‑100 (MQ=255)
                                                                       aTTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATg                             >  2:75429/1‑100 (MQ=255)
                                                                         tAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                           >  1:1605888/1‑100 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                           <  1:1262263/99‑1 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                          >  1:2108003/1‑100 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                          >  2:2515449/1‑100 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTc                          >  1:1362669/1‑100 (MQ=255)
                                                                          aaCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGAGAAATCAGACGCCGATATTTCTCCACTTCTCGCCTGCCATCAGTTTTCGTTGGATATGTTc                          >  2:2582260/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                       <  1:1768077/100‑1 (MQ=255)
                                                                             gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                       >  2:1511160/1‑100 (MQ=255)
                                                                               aaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                     >  2:3299560/1‑100 (MQ=255)
                                                                               aaaaaGGGCCGGATGGACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCg                     <  2:2269622/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                <  1:2996665/100‑1 (MQ=255)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGACATCAGACGCCGTTATTTCTCAACTTCTCGCCTGCCATCAATCTTCATTCGATAAGTTCCACCGAGACa                >  2:237035/1‑99 (MQ=25)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt               >  2:2560480/1‑100 (MQ=39)
                                                                                     ggCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAt               >  2:3876603/1‑100 (MQ=39)
                                                                                      gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGCCGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATGTGTTCCAGGATGACAtt              >  2:508963/1‑100 (MQ=255)
                                                                                      gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt              >  2:135706/1‑100 (MQ=255)
                                                                                      gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACAtt              >  1:2193352/1‑100 (MQ=255)
                                                                                          gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTg          <  2:3262220/100‑1 (MQ=255)
                                                                                            tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt        <  1:4035147/100‑1 (MQ=255)
                                                                                            tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt        >  2:4035120/1‑100 (MQ=255)
                                                                                            tGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGt        <  1:3466757/100‑1 (MQ=255)
                                                                                             gTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGtt       >  1:3565292/1‑100 (MQ=255)
                                                                                              tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGttt      >  1:2538300/1‑100 (MQ=255)
                                                                                                catcacaTCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTcc    <  2:2932057/97‑1 (MQ=255)
                                                                                                  gCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCgg  >  1:527516/1‑100 (MQ=255)
                                                                                                  gCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCgg  <  1:3350003/100‑1 (MQ=255)
                                                                                                  gCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCgg  >  2:1845094/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGG  >  NC_000913/2980647‑2980844

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: