Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,903,571 A→G 36.1% T164A (ACC→GCC)  manY → mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,903,5710AG36.1% 113.2 / 121.3 119T164A (ACC→GCC) manYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (32/44);  new base G (22/21);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 122 total aligned reads displayed.

GGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903472‑1903670
                                                                                                   |                                                                                                   
ggCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTa                                                                                                     <  2:1637628/100‑1 (MQ=255)
  cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGcc                                                                                                   <  1:1694451/100‑1 (MQ=255)
  cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                                                   >  2:2192252/1‑100 (MQ=255)
  cAACCTGACAGCGATCTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                                                   <  1:2807669/100‑1 (MQ=255)
   aaCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCa                                                                                                  <  1:2651980/100‑1 (MQ=255)
    aCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAg                                                                                                 <  1:2141243/100‑1 (MQ=255)
    aCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAg                                                                                                 <  1:1907175/100‑1 (MQ=255)
     ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGc                                                                                                <  1:2518725/100‑1 (MQ=255)
     ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGc                                                                                                <  2:2291931/100‑1 (MQ=255)
       tGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGa                                                                                              <  1:2579204/100‑1 (MQ=255)
         aCAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAg                                                                                            <  1:1537692/100‑1 (MQ=255)
         aCAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAg                                                                                            <  2:727500/100‑1 (MQ=255)
            gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTAc                                                                                         >  2:783224/1‑100 (MQ=255)
             cGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTAc                                                                                         >  2:2904509/1‑99 (MQ=255)
             cGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACa                                                                                        <  1:2782966/100‑1 (MQ=255)
             cGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACa                                                                                        <  2:2817351/100‑1 (MQ=255)
               aTTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGa                                                                                      <  2:3302453/100‑1 (MQ=255)
               aTTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGa                                                                                      <  1:2043269/100‑1 (MQ=255)
                 ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                                                    <  1:3050799/100‑1 (MQ=255)
                 ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                                                    <  1:2443190/100‑1 (MQ=255)
                 ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                                                    <  2:2363344/100‑1 (MQ=255)
                   ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACAt                                                                                  <  2:2489856/100‑1 (MQ=255)
                   ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACAt                                                                                  >  1:3363058/1‑100 (MQ=255)
                    cTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATg                                                                                 >  1:3412581/1‑100 (MQ=255)
                     tGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGc                                                                                <  1:2228322/100‑1 (MQ=255)
                       gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTg                                                                              <  2:2699223/100‑1 (MQ=255)
                       gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                                              <  1:710193/100‑1 (MQ=255)
                       gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                                              >  1:2195744/1‑100 (MQ=255)
                       gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                                              <  1:1499305/100‑1 (MQ=255)
                            aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGc                                                                         >  2:2219355/1‑100 (MQ=255)
                            aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGc                                                                         >  2:1332452/1‑100 (MQ=255)
                            aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGc                                                                         <  2:1857386/100‑1 (MQ=255)
                            aCGTTTCGTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                         <  1:1860015/100‑1 (MQ=255)
                               tttCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGAt                                                                      <  1:785393/100‑1 (MQ=255)
                                   ttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCg                                                                  <  2:1893605/100‑1 (MQ=255)
                                    tctcTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                                                 >  1:171905/1‑100 (MQ=255)
                                    tctcTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                                                 >  2:377617/1‑100 (MQ=255)
                                     ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                                                <  2:1744065/100‑1 (MQ=255)
                                         gTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAggg                                                            >  2:3436617/1‑98 (MQ=255)
                                         gTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAg                                                              >  2:3037594/1‑98 (MQ=255)
                                            ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         >  1:214177/1‑100 (MQ=255)
                                             cTGCAAGCAATGCGTGTGGCTCTTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  1:2191402/1‑100 (MQ=255)
                                             cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  1:1250640/1‑100 (MQ=255)
                                             cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  2:1153935/1‑100 (MQ=255)
                                              tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       >  2:2721578/1‑100 (MQ=255)
                                              tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       >  1:2634092/1‑100 (MQ=255)
                                              tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGCATGCGATTCCGGAAGTGGTGAc                                                       <  1:2007045/100‑1 (MQ=255)
                                              tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       >  1:2069054/1‑100 (MQ=255)
                                              tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       <  1:248334/100‑1 (MQ=255)
                                               gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      >  2:3156723/1‑100 (MQ=255)
                                               gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      >  1:1559378/1‑100 (MQ=255)
                                                 aagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCaa                                                    >  2:1448204/1‑100 (MQ=255)
                                                  agcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                                                   >  1:1487866/1‑100 (MQ=255)
                                                   gcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATg                                                  >  2:1228876/1‑100 (MQ=255)
                                                    caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                                                 <  1:798207/100‑1 (MQ=255)
                                                    caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                                                 <  1:438208/100‑1 (MQ=255)
                                                      aTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                               <  1:672898/100‑1 (MQ=255)
                                                        gCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTg                                             <  2:356936/100‑1 (MQ=255)
                                                         cGTGTGTCTATTGCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGa                                            <  2:3232828/100‑1 (MQ=255)
                                                           tgtgGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat                                          >  2:50110/1‑100 (MQ=255)
                                                           tgtgGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat                                          >  1:435019/1‑100 (MQ=255)
                                                              ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                       >  2:1194014/1‑100 (MQ=255)
                                                              ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                       <  1:2561889/100‑1 (MQ=255)
                                                              ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                       >  2:121824/1‑100 (MQ=255)
                                                               gCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg                                      <  1:1958275/100‑1 (MQ=255)
                                                               gCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg                                      <  2:3424898/100‑1 (MQ=255)
                                                                cTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGc                                     <  2:3008663/100‑1 (MQ=255)
                                                                 tATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc                                    <  1:3436589/100‑1 (MQ=255)
                                                                 tATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc                                    <  2:1096078/100‑1 (MQ=255)
                                                                   ttCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg                                  >  2:340725/1‑100 (MQ=255)
                                                                   ttCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg                                  >  1:418109/1‑100 (MQ=255)
                                                                     ccGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                                <  2:2195761/100‑1 (MQ=255)
                                                                       ggCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATTCTTAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                                <  1:177983/98‑1 (MQ=255)
                                                                         ccGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            <  1:102478/100‑1 (MQ=255)
                                                                         ccGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            <  1:2219337/100‑1 (MQ=255)
                                                                          cggtaTCGTGGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATg                           <  1:2735136/96‑1 (MQ=255)
                                                                          cGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATg                           >  2:208308/1‑100 (MQ=255)
                                                                          cGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGAGTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATg                           >  1:1547948/1‑100 (MQ=255)
                                                                            tCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACGAATGGTCTGAATATCGCCGGTGGCATGAt                         >  1:2075555/1‑100 (MQ=255)
                                                                            tCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGAt                         <  1:166804/100‑1 (MQ=255)
                                                                             cATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATc                        >  2:2359925/1‑100 (MQ=255)
                                                                             cATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATc                        >  2:2481541/1‑100 (MQ=255)
                                                                              aTCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCg                       <  1:2479119/100‑1 (MQ=255)
                                                                              aTCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCg                       <  1:3156699/100‑1 (MQ=255)
                                                                                 gTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtgg                    <  1:809829/100‑1 (MQ=255)
                                                                                  ttGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  2:887642/100‑1 (MQ=255)
                                                                                  ttGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  2:1404948/100‑1 (MQ=255)
                                                                                    gcgcTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGtggtgg                 >  1:3290564/1‑100 (MQ=255)
                                                                                         tcttgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  2:3384244/97‑1 (MQ=255)
                                                                                         gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  2:1345849/100‑1 (MQ=255)
                                                                                            tgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggtaga         >  1:3384219/1‑97 (MQ=255)
                                                                                             gtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATg        >  1:3407495/1‑100 (MQ=255)
                                                                                             gtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATg        >  2:511486/1‑100 (MQ=255)
                                                                                              ttGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGc       <  1:217665/100‑1 (MQ=255)
                                                                                               tGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCg      >  2:824964/1‑100 (MQ=255)
                                                                                                 gTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGAt    >  1:1171155/1‑100 (MQ=255)
                                                                                                  tACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATg   >  1:1077658/1‑100 (MQ=255)
                                                                                                  tACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAAGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATg   <  2:1373226/100‑1 (MQ=255)
                                                                                                   aCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATgg  >  2:3073501/1‑100 (MQ=255)
                                                                                                   aCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATgg  <  1:1153924/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903472‑1903670

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