Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,310,468 T→G 41.1% E171D (GAA→GAC alsA ← D‑allose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,310,4680TG41.1% 88.2 / 197.6 160E171D (GAA→GACalsAD‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (55/39);  new base G (38/28);  total (93/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

TTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGA  >  NC_000913/4310370‑4310567
                                                                                                  |                                                                                                   
ttATGCAAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGtt                                                                                                    <  2:788020/100‑1 (MQ=255)
   tGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCAt                                                                                                 >  2:838415/1‑100 (MQ=255)
    gCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATc                                                                                                >  2:684628/1‑100 (MQ=255)
     cGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCTACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATcc                                                                                               >  2:3005401/1‑100 (MQ=255)
     cGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATcc                                                                                               <  2:1595543/100‑1 (MQ=255)
     cGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATcc                                                                                               >  2:3224022/1‑100 (MQ=255)
      gagaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCa                                                                                              <  1:1816607/100‑1 (MQ=255)
       agaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCat                                                                                             >  1:10360/1‑100 (MQ=255)
       agaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCat                                                                                             <  2:3549199/100‑1 (MQ=255)
        gaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatg                                                                                            <  1:613787/100‑1 (MQ=255)
         atatAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatga                                                                                           <  2:3603819/100‑1 (MQ=255)
          tctaGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGATACAGATAGTCCACCTCTTTATTGGTGAGTGAAGAGGTGGGGTCAGCCatga                                                                                           >  1:3721904/3‑99 (MQ=255)
          tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatgat                                                                                          >  2:2414804/1‑100 (MQ=255)
          tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCatgat                                                                                          <  1:3902925/100‑1 (MQ=255)
           ataGACGATGGCCGTACCCTCTTTACGCACCTGCTTCATGCTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGTAGGTGGGGTCATCCatgat                                                                                          >  1:3824943/1‑99 (MQ=255)
            taGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatgatga                                                                                        >  1:4157506/1‑100 (MQ=255)
            taGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCatgatg                                                                                         >  1:3256998/1‑99 (MQ=255)
             aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatgatgat                                                                                       >  1:2386546/1‑100 (MQ=255)
             aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatgatgat                                                                                       <  2:473353/100‑1 (MQ=255)
             aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCatgatgat                                                                                       >  2:1098366/1‑100 (MQ=255)
              gACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGCGGTGGGTTCATCCatgatga                                                                                        >  1:3352096/1‑98 (MQ=255)
                   tGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACttt                                                                                 >  1:2292268/1‑100 (MQ=255)
                     gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTgg                                                                               <  1:720832/100‑1 (MQ=255)
                     gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTgg                                                                               <  1:3022896/100‑1 (MQ=255)
                     gCCGTACCCTCTTTACGCAACTAATTCATGATCAGAAACAGATAGTCCACCTCTTTATTTGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTgg                                                                               >  1:309239/1‑100 (MQ=255)
                      ccGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGc                                                                              >  2:3982853/1‑100 (MQ=255)
                          aCCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCGCCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCg                                                                          <  1:2688130/100‑1 (MQ=255)
                           cccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGa                                                                         >  2:3472981/1‑100 (MQ=255)
                            ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAg                                                                        <  1:1076007/100‑1 (MQ=255)
                            catctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAg                                                                        <  1:2580687/98‑1 (MQ=255)
                               ctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGc                                                                       <  2:497339/98‑1 (MQ=255)
                                 ttACGCAACTGATTCATGATCAGAATCAGATAGTCCACCTCTTTATTGGTTAGTGTGTAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATc                                                                    <  1:4258416/99‑1 (MQ=255)
                                 ttACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCa                                                                   >  1:1683114/1‑100 (MQ=255)
                                   aCGCAACTGATTCATGTTCAGAAACAGATAGTGCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCTATGATGATGACTTTGGCATCGAGCATCAg                                                                  <  1:272875/99‑1 (MQ=255)
                                    cGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCg                                                                <  2:2493631/100‑1 (MQ=255)
                                    cGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCg                                                                <  2:365181/100‑1 (MQ=255)
                                    cGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCg                                                                <  2:2368672/100‑1 (MQ=255)
                                      cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGCTGACTTTGGCATCTAGCATCAGCGtt                                                              >  1:3972630/1‑100 (MQ=255)
                                      cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGAGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGtt                                                              >  1:1281254/1‑100 (MQ=255)
                                         cTGATTCATGATCAGAAACAGATAGTCCACCTCGTTATTGGTGAGTGTGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGttt                                                             <  1:1912198/98‑1 (MQ=255)
                                          tGATTCATGATCAGAAACAGATAGTCCTCCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTgg                                                          <  2:928206/100‑1 (MQ=255)
                                          tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTgg                                                          <  2:2396210/100‑1 (MQ=255)
                                          tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTgg                                                          <  2:1733205/100‑1 (MQ=255)
                                           gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGc                                                         >  2:4106461/1‑100 (MQ=255)
                                           gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGc                                                         >  2:362100/1‑100 (MQ=255)
                                           gATTCATGATCAGAAACAGATAGTCCACCTCTTCATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTgg                                                          <  2:4183980/99‑1 (MQ=255)
                                            aTTCATGAGCAGAAACAGATAGTCCACCTCTTTATTTGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCa                                                        <  1:2361778/100‑1 (MQ=255)
                                             ttCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCaa                                                       >  2:2411455/1‑100 (MQ=255)
                                              tCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAt                                                      >  1:2155801/1‑100 (MQ=255)
                                              tCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAt                                                      >  1:957477/1‑100 (MQ=255)
                                               cctGATCAGAAACGGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAtt                                                     >  2:2902580/3‑100 (MQ=255)
                                               cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAtt                                                     <  2:2922404/100‑1 (MQ=255)
                                               cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAtt                                                     >  2:3525332/1‑100 (MQ=255)
                                                  gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                  >  2:1168102/1‑100 (MQ=255)
                                                  gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                  >  2:4244064/1‑100 (MQ=255)
                                                  gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                  >  1:2609274/1‑100 (MQ=255)
                                                   aTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTa                                                 >  1:3484418/1‑100 (MQ=255)
                                                   aTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATAGAGCATCAGCGTTTTGGCAATTTCTa                                                 >  1:2192845/1‑100 (MQ=255)
                                                    tCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAg                                                <  2:1432309/100‑1 (MQ=255)
                                                      aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGcn                                              <  1:4243771/100‑2 (MQ=255)
                                                      aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCa                                              >  1:2110095/1‑100 (MQ=255)
                                                        aaaCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATc                                            <  1:4201464/100‑1 (MQ=255)
                                                          aCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTg                                          >  1:3249049/1‑100 (MQ=255)
                                                              aTAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtg                                      <  2:341254/100‑1 (MQ=255)
                                                               tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgt                                     >  2:60567/1‑100 (MQ=255)
                                                               tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgt                                     >  1:1307339/1‑100 (MQ=255)
                                                                 gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTgg                                   >  2:14608/1‑100 (MQ=255)
                                                                  tCCACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGc                                  <  2:1977428/100‑1 (MQ=255)
                                                                   ccACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                 >  2:3294748/1‑100 (MQ=255)
                                                                   ccACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                 >  1:451350/1‑100 (MQ=255)
                                                                   ccACCTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                 >  2:1712002/1‑100 (MQ=255)
                                                                      ccTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                              >  1:600529/1‑100 (MQ=255)
                                                                      ccTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                              <  2:3534465/100‑1 (MQ=255)
                                                                      ccTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                              <  1:3803337/100‑1 (MQ=255)
                                                                      ccTCTTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                              >  2:1456528/1‑100 (MQ=255)
                                                                       ctctTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                             <  1:1731935/100‑1 (MQ=255)
                                                                       ctctTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                             <  2:3738012/100‑1 (MQ=255)
                                                                        tctTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatag                            <  2:716215/100‑1 (MQ=255)
                                                                         ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatag                            <  2:1033288/99‑1 (MQ=255)
                                                                         ctTTATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgataga                           >  1:2980728/1‑100 (MQ=255)
                                                                          tttATTGGTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagat                          <  2:3452962/100‑1 (MQ=255)
                                                                           ttATTGGTGAGTGAGTAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagata                         <  2:4217868/100‑1 (MQ=255)
                                                                             aTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaaa                       >  1:1376159/1‑100 (MQ=255)
                                                                                ggTGAGTGAGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                    >  2:4190170/1‑100 (MQ=255)
                                                                                 gtgagtgaGGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCg                   <  1:4043999/100‑1 (MQ=255)
                                                                                  tgagtgagGAGGTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGc                  >  1:2227154/1‑100 (MQ=255)
                                                                                     gtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCAc               >  1:2168690/1‑100 (MQ=255)
                                                                                     gtgagGAGGAGGGGTCATACATGATGATGACTTTGGCATCGAGAATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTAGCTAc               >  1:483042/1‑100 (MQ=255)
                                                                                     atgagGAGCTGGCGTCATCTATGATGATGACATTGGAATCGAGCATTAGCGTTTTGGCAATTTCTAGCATCTTCTTGTTGCTGATAGATAAATTCGc                  >  1:787735/2‑97 (MQ=255)
                                                                                       gtggaggTGGGTTAATCCATGATGATGACTTGGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt             >  2:2876321/3‑100 (MQ=255)
                                                                                        aggatgTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt             <  1:912308/99‑1 (MQ=255)
                                                                                         ggaggTGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc           >  2:435161/1‑100 (MQ=255)
                                                                                           aggTGGGTTTATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtctc         <  1:1042894/100‑1 (MQ=255)
                                                                                              tGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc      >  1:180345/1‑100 (MQ=255)
                                                                                              tGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc      >  1:3394973/1‑100 (MQ=255)
                                                                                              tGGGGTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc      >  1:3119505/1‑100 (MQ=255)
                                                                                               gggTTCATCCATGATGATGACTTTGGAATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCt     >  2:2597164/1‑100 (MQ=255)
                                                                                                 gTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAg   >  2:3745297/1‑100 (MQ=255)
                                                                                                  ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGa  <  2:2561371/100‑1 (MQ=255)
                                                                                                  ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGa  >  1:814324/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGA  >  NC_000913/4310370‑4310567

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: