Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,903,571 A→G 61.7% 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,5710AG61.7% 106.8 / 134.8 133T164A (ACC→GCC) manYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base G (32/50);  total (53/80)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 134 total aligned reads displayed.

CAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903474‑1903670
                                                                                                 |                                                                                                   
cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGcc                                                                                                   <  1:3016003/100‑1 (MQ=255)
cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                                                   <  1:3008027/100‑1 (MQ=255)
cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                                                   <  2:314839/100‑1 (MQ=255)
  aCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAg                                                                                                 >  2:579159/1‑100 (MQ=255)
  aCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAg                                                                                                 >  2:4315142/1‑100 (MQ=255)
  aCCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAg                                                                                                 >  1:2982747/1‑100 (MQ=255)
    cTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCg                                                                                               <  1:2781050/100‑1 (MQ=255)
     tGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGTCGTCATCGTTGCGCTGTCGGTTGGTACCAGCGa                                                                                              <  1:752301/100‑1 (MQ=255)
      gACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGaa                                                                                             >  2:2827071/1‑100 (MQ=255)
      gACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGaa                                                                                             <  1:1186169/100‑1 (MQ=255)
      gACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGaa                                                                                             >  1:3090866/1‑100 (MQ=255)
       actgCGATGTCCTGGATGCAGGTTTCTTCTCTTTTCCTGTAAGCACTGCGTGCGGCTATTCCGGCCGTGATCGTTGCGCTGTCTGTTGGAGCCAGCGAAg                                                                                            <  1:958880/97‑1 (MQ=255)
       aCAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAg                                                                                            <  2:1950819/100‑1 (MQ=255)
        cAGTGATTTCCTGGATTCACGTTTCTTCTCTGTGCCTGCAAGCAATGCGTGTGTCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAg                                                                                            <  1:1979911/99‑1 (MQ=255)
        cAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGt                                                                                           >  1:1689173/1‑100 (MQ=255)
        cAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                                                           >  2:3156391/1‑100 (MQ=255)
          gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                         <  1:4065028/100‑1 (MQ=255)
          gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                                                         >  1:2756841/1‑100 (MQ=255)
          gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGGACCAGCGAAGTAc                                                                                         <  2:2965080/100‑1 (MQ=255)
           cGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGGGCCAGCGAAGTa                                                                                          >  1:1579353/1‑98 (MQ=255)
            gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAg                                                                                       <  2:476150/100‑1 (MQ=255)
            gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAg                                                                                       >  2:1814999/1‑100 (MQ=255)
              tttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCCTCGTTGCGCTGTCTGCTGGTGCCAGCGAAGTACAg                                                                                       >  1:3520938/1‑98 (MQ=255)
              tttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGaa                                                                                     <  2:4176737/100‑1 (MQ=255)
               ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGGGGGGGGCTATCCGGGCCGTCACCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAAc                                                                                    <  1:4407254/100‑1 (MQ=255)
               ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                                                    <  2:1356439/100‑1 (MQ=255)
               ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGCCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAAc                                                                                    >  1:2592516/1‑100 (MQ=255)
                 ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACAt                                                                                  >  2:548610/1‑100 (MQ=255)
                     gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTg                                                                              <  2:3427610/100‑1 (MQ=255)
                     gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTg                                                                              >  2:4202354/1‑100 (MQ=255)
                        ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAAt                                                                           >  2:4357031/1‑100 (MQ=255)
                         cACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATg                                                                          >  2:480769/1‑100 (MQ=255)
                         cACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATg                                                                          >  2:4203299/1‑100 (MQ=255)
                          aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                         <  2:5835/100‑1 (MQ=255)
                          aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                         <  2:1689180/100‑1 (MQ=255)
                           cGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCg                                                                        <  1:1411841/100‑1 (MQ=255)
                           cGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCACCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCg                                                                        <  1:998686/100‑1 (MQ=255)
                               tcttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTc                                                                    <  2:2451802/100‑1 (MQ=255)
                               tcttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTc                                                                    <  2:636760/100‑1 (MQ=255)
                                cttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTcc                                                                   <  1:2582752/100‑1 (MQ=255)
                                 ttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTccn                                                                  <  1:1665469/100‑2 (MQ=255)
                                   ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                                                >  2:4075096/1‑100 (MQ=255)
                                   ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                                                 >  2:3977955/1‑99 (MQ=255)
                                   ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                                                <  2:3149380/100‑1 (MQ=255)
                                    tctTTTCCTTCAAGCAATGCGTATGGCTATTCCGGCAGTCATCGTTGCGTTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGaa                                                               <  1:3380175/100‑1 (MQ=255)
                                    tctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGaa                                                               <  1:4202323/100‑1 (MQ=255)
                                       gTTCCTGCAAGCAATGCGTTTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtg                                                            <  1:2331022/100‑1 (MQ=255)
                                        ttCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtgg                                                           <  2:1734275/100‑1 (MQ=255)
                                         tCCTGCAAGCAATGCGTGTGGCTATTCTGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                          >  1:3435489/1‑100 (MQ=255)
                                          ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         >  1:3835066/1‑100 (MQ=255)
                                          ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         <  2:1897434/100‑1 (MQ=255)
                                          ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         <  2:1551070/100‑1 (MQ=255)
                                           cggcaagcaaTGCGTGTGGCTATTCCGGCCGTTATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        <  1:74831/98‑1 (MQ=255)
                                           cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         <  1:1000308/99‑1 (MQ=255)
                                           cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  1:2043034/1‑100 (MQ=255)
                                           cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  2:3407251/1‑100 (MQ=255)
                                            tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       >  2:1015049/1‑100 (MQ=255)
                                             gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      <  1:1255381/100‑1 (MQ=255)
                                             gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      <  2:3573536/100‑1 (MQ=255)
                                             gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      <  2:2135383/100‑1 (MQ=255)
                                                agcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                                                   >  2:3657110/1‑100 (MQ=255)
                                                  caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                                                 <  1:3786878/100‑1 (MQ=255)
                                                  caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGCAGTGGTGACCAATgg                                                 >  1:653044/1‑100 (MQ=255)
                                                     tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCt                                              <  2:3062345/100‑1 (MQ=255)
                                                           tgGCTCTTCCTGCCGTCCTCGGTGCGCTGCCTGTTGGTACCAGCGAAGTACAAAACATGCTGACTGCGATTCCGGAAGTGGTGAACAATGGTCTGAata                                         >  1:146530/1‑99 (MQ=255)
                                                            ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAatata                                       >  2:1988047/1‑99 (MQ=255)
                                                             gCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg                                      >  2:2199105/1‑100 (MQ=255)
                                                                aTTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCg                                   >  2:2363547/1‑100 (MQ=255)
                                                                   ccGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                                >  2:1275772/1‑100 (MQ=255)
                                                                   ccGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCACCGGTg                                >  2:839225/1‑100 (MQ=255)
                                                                    cGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                               <  2:352499/100‑1 (MQ=255)
                                                                    cGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                               <  2:3854409/100‑1 (MQ=255)
                                                                    cGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCGGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                               <  1:274619/100‑1 (MQ=255)
                                                                     ggCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                               <  1:1321020/99‑1 (MQ=255)
                                                                     ggCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc                              <  1:1688954/100‑1 (MQ=255)
                                                                     ggCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc                              >  1:1654421/1‑100 (MQ=255)
                                                                      gCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCa                             >  2:2192849/1‑100 (MQ=255)
                                                                        cGTCATCGTTGCGCTGTCGGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            <  2:4017693/99‑1 (MQ=255)
                                                                          tCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGAt                         <  1:4396167/100‑1 (MQ=255)
                                                                            aTCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCg                       <  1:323846/100‑1 (MQ=255)
                                                                              cGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtg                     <  2:1526221/100‑1 (MQ=255)
                                                                              cGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtg                     <  1:2241835/100‑1 (MQ=255)
                                                                              cGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtg                     <  2:739028/100‑1 (MQ=255)
                                                                               gTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtgg                    >  2:1985681/1‑100 (MQ=255)
                                                                               gTTGCGATGTCTGTTGGGGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGGGCGGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtgg                    <  1:4203268/100‑1 (MQ=255)
                                                                                ttGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  2:4318689/100‑1 (MQ=255)
                                                                                ttGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  1:3317143/100‑1 (MQ=255)
                                                                                 tGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggtg                  <  2:2627241/100‑1 (MQ=255)
                                                                                  gcgcTGTCTGTTGGTACCAGCGAAGTACAGAACAGGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATCATCGtggtgg                 <  1:3778728/100‑1 (MQ=255)
                                                                                   cgcTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCTCCGGTGGCATGATCGTggtggt                <  1:2560906/100‑1 (MQ=255)
                                                                                    gcTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggtt               >  1:3878940/1‑100 (MQ=255)
                                                                                     ctgtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttg              <  2:1760768/100‑1 (MQ=255)
                                                                                      tgtctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttgg             >  1:4044073/1‑100 (MQ=255)
                                                                                        tctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt           <  2:3920290/100‑1 (MQ=255)
                                                                                        tctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt           <  2:3640601/100‑1 (MQ=255)
                                                                                         ctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa          <  2:3998524/100‑1 (MQ=255)
                                                                                         ctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa          >  2:1825429/1‑100 (MQ=255)
                                                                                           gtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATg        <  1:4239696/100‑1 (MQ=255)
                                                                                              ggTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGa     >  2:1298631/1‑100 (MQ=255)
                                                                                                 aCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATgg  <  1:445352/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903474‑1903670

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: