Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,116 T→G Y122S (TAT→TCT)  thiB ← thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,1160TG94.1% 406.5 / 12.8 136Y122S (TAT→TCT) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (63/65);  total (64/72)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 137 total aligned reads displayed.

TAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTG  >  NC_000913/75017‑75211
                                                                                                   |                                                                                               
tAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAt                                                                                                 >  2:582628/1‑100 (MQ=255)
   aTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAAt                                                                                              >  1:3642139/1‑100 (MQ=255)
    tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATc                                                                                             <  2:1553225/100‑1 (MQ=255)
    tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATc                                                                                             <  1:1539374/100‑1 (MQ=255)
     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCa                                                                                            <  2:974805/100‑1 (MQ=255)
     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCa                                                                                            <  2:40277/100‑1 (MQ=255)
     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCa                                                                                            >  1:681451/1‑100 (MQ=255)
       cccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCaaa                                                                                          >  2:1489879/1‑100 (MQ=255)
       cccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCGTGGTGGCGGGTTTTTCAGTGTGGGCTCGGCATAAACGAAGGCTAAGGATCCAGAATCaaa                                                                                          >  2:3023841/1‑100 (MQ=255)
        ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAAc                                                                                         <  1:1441808/100‑1 (MQ=255)
        ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAAc                                                                                         <  1:714478/100‑1 (MQ=255)
        ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAAc                                                                                         <  1:2179689/100‑1 (MQ=255)
         cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACg                                                                                        <  1:1615404/100‑1 (MQ=255)
             aGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTAc                                                                                    >  1:3660013/1‑100 (MQ=255)
              gTTTTGATCGCTCTCAACCAGTTCTATCAGGCTTTGTGGCGGGGGTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACg                                                                                   >  1:1299935/1‑100 (MQ=255)
                tttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGaa                                                                                 >  1:2084028/1‑100 (MQ=255)
                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAg                                                                               <  2:1460680/100‑1 (MQ=255)
                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAg                                                                               >  2:2009250/1‑100 (MQ=255)
                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAg                                                                               <  1:3411649/100‑1 (MQ=255)
                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAg                                                                               <  1:2631744/100‑1 (MQ=255)
                   gATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAgt                                                                              >  1:3441411/1‑100 (MQ=255)
                    aTCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAgtg                                                                             >  2:3921254/1‑100 (MQ=255)
                        ctctCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCat                                                                         <  1:502410/100‑1 (MQ=255)
                            cAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGGCcttatt                                                                     >  1:135840/1‑100 (MQ=255)
                            cAACCAGTTCTATCAGGCTTTGTGGCGGGTTTTTCAGGTTGTTCTTGTCATACAAGACGGCGAAGTAGCCAGCATCAAACGGTACGAAAGCGTCattatt                                                                     >  2:3444295/1‑100 (MQ=255)
                             aaCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTc                                                                    <  2:116216/100‑1 (MQ=255)
                             aaCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTc                                                                    >  2:1295383/1‑100 (MQ=255)
                              aCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTcc                                                                   >  2:2870009/1‑100 (MQ=255)
                              aCCAGTTCTTTCACGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATCAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTcc                                                                   <  1:3681373/100‑1 (MQ=255)
                                cAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAg                                                                 >  1:94005/1‑100 (MQ=255)
                                cAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAg                                                                 <  2:405102/100‑1 (MQ=255)
                                  gTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGACAGTGTCATTATTCCAgcc                                                               >  1:3824209/1‑100 (MQ=255)
                                     cTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAgccggc                                                            <  2:3706068/100‑3 (MQ=255)
                                     cTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAgccgcc                                                            >  2:3968972/1‑100 (MQ=255)
                                     cTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAgccgcc                                                            >  2:2586198/1‑100 (MQ=255)
                                     cTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAgccgcc                                                            >  1:3645799/1‑100 (MQ=255)
                                       ttCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGccgccgn                                                          <  2:1110041/100‑2 (MQ=255)
                                       ttCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCgg                                                          <  1:2592189/100‑1 (MQ=255)
                                       ttCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCgg                                                          <  2:830268/100‑1 (MQ=255)
                                         cAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGa                                                        >  2:108572/1‑100 (MQ=255)
                                           ggCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGGGAAGTAGCCAGAATCAAACGGTACGGAAGTGTCATTATTCCAGCCGCCGGGAAc                                                      >  2:2787358/1‑100 (MQ=255)
                                           ggCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAAc                                                      <  1:55034/100‑1 (MQ=255)
                                           ggCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAAc                                                      >  2:2474573/1‑100 (MQ=255)
                                            gCTTTGTGGCGGGTTTTTCAGTTTGTTCTGGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACg                                                     <  1:236414/100‑1 (MQ=255)
                                             cTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGt                                                    >  1:2503470/1‑100 (MQ=255)
                                             cTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGt                                                    >  1:2375775/1‑100 (MQ=255)
                                             cTTGGTGGCGGGTTTTTCAGTTTGTTCTTGGCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGt                                                    <  1:2663294/100‑1 (MQ=255)
                                                 gtgGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAAc                                                >  2:1466687/1‑100 (MQ=255)
                                                  tggtggGTTTCTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCAGCCGCCGGGAACGTTAACgg                                              >  2:3677036/1‑100 (MQ=255)
                                                   ggCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACgg                                              >  1:1253274/1‑100 (MQ=255)
                                                    gCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACggn                                             <  1:3248931/100‑2 (MQ=255)
                                                     cGGGTTTTTCAGTTTTTTCTCGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCa                                            >  2:1156617/1‑100 (MQ=255)
                                                     cGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCa                                            >  1:2894146/1‑100 (MQ=255)
                                                       ggTTTTTCAGTTTGTTCTGGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCAta                                          <  2:3621139/100‑2 (MQ=255)
                                                        gTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc                                         >  2:2457097/1‑100 (MQ=255)
                                                         tttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCg                                        <  1:154033/100‑1 (MQ=255)
                                                          ttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGc                                       >  1:3216524/1‑100 (MQ=255)
                                                            ttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTg                                     >  2:571004/1‑100 (MQ=255)
                                                             tCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGTCACAAGCGGTACGAATGTGTCATTTTTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGc                                    <  1:1424261/100‑1 (MQ=255)
                                                             tCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGc                                    <  2:3342329/100‑1 (MQ=255)
                                                              cAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTCTTCCAGCCGCCGGGAACGTTAACGGCATCCGCGGcc                                   >  2:3618373/1‑100 (MQ=255)
                                                              cAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGcc                                   <  2:2855551/100‑1 (MQ=255)
                                                              cAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGcc                                   <  2:2051664/100‑1 (MQ=255)
                                                                gTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCcac                                 <  1:1477947/100‑1 (MQ=255)
                                                                gTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCcac                                 <  1:2513225/100‑1 (MQ=255)
                                                                  ttGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCacac                               <  2:945293/100‑1 (MQ=255)
                                                                   tGTTCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACAcc                              <  1:1479687/100‑1 (MQ=255)
                                                                      tCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCt                           <  1:870472/100‑1 (MQ=255)
                                                                      tCTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCt                           >  2:45230/1‑100 (MQ=255)
                                                                       cTTGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCtt                          >  2:877155/1‑100 (MQ=255)
                                                                        ttGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCttt                         >  2:680160/1‑100 (MQ=255)
                                                                         tGTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCtttt                        >  1:2796970/1‑100 (MQ=255)
                                                                          gTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTg                       <  2:1627638/100‑1 (MQ=255)
                                                                          gTCATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTg                       >  1:1344365/1‑100 (MQ=255)
                                                                            cATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGc                     <  1:3337751/100‑1 (MQ=255)
                                                                            cATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGc                     >  2:1871297/1‑100 (MQ=255)
                                                                            cATAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGc                     >  1:654547/1‑100 (MQ=255)
                                                                             aTAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCGTTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCa                    <  2:3262047/100‑1 (MQ=255)
                                                                             aTAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCa                    >  1:81175/1‑100 (MQ=255)
                                                                             aTAAACGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCa                    >  2:3172410/1‑100 (MQ=255)
                                                                                aaCGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAAc                 <  2:91737/100‑1 (MQ=255)
                                                                                aaCGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAAc                 >  2:1699970/1‑100 (MQ=255)
                                                                                  cGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg               <  1:1154936/100‑1 (MQ=255)
                                                                                  cGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg               <  1:2495458/100‑1 (MQ=255)
                                                                                  cGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg               <  1:2748944/100‑1 (MQ=255)
                                                                                  cGAAGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAg               <  1:2811910/100‑1 (MQ=255)
                                                                                    aaGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTc             <  1:448733/100‑1 (MQ=255)
                                                                                     aGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc            >  1:359820/1‑100 (MQ=255)
                                                                                     aGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc            <  1:721553/100‑1 (MQ=255)
                                                                                     aGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc            <  1:2141736/100‑1 (MQ=255)
                                                                                     aGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc            >  1:853559/1‑100 (MQ=255)
                                                                                     aGGCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc            >  1:1405117/1‑100 (MQ=255)
                                                                                       gCGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCgg          <  2:3589886/100‑1 (MQ=255)
                                                                                        cGAAGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGt         >  2:2803935/1‑100 (MQ=255)
                                                                                          aagaagCAAGAATCAAACGGGACGAAAGTGTCTTTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGttt       <  2:3429738/100‑1 (MQ=255)
                                                                                          aaGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGttt       <  2:1553580/100‑1 (MQ=255)
                                                                                          aaGTAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGttt       >  1:924151/1‑100 (MQ=255)
                                                                                             tAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTAc    >  1:2886627/1‑100 (MQ=255)
                                                                                             tAGCCAGAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTAc    <  1:1087252/100‑1 (MQ=255)
                                                                                               gCCAGAATCAAACGGTACGAAAGTGTCATTATTCCGGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTg  <  2:446411/100‑1 (MQ=255)
                                                                                                   |                                                                                               
TAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTG  >  NC_000913/75017‑75211

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: