Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,531,118 Δ1 bp 15.3% coding (534/1122 nt) sgcX ← putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,531,1180T.15.3% 432.2 / 76.7 150coding (534/1122 nt)sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (68/59);  new base . (14/9);  total (82/68)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-02

Only 100 of 159 total aligned reads displayed.

CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTG  >  NC_000913/4531019‑4531217
                                                                                                    |                                                                                                   
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGt                                                                                                     >  2:3155488/1‑100 (MQ=255)
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGt                                                                                                     >  1:1738975/1‑100 (MQ=255)
    aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCtt                                                                                                 >  1:1405999/1‑100 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCttt                                                                                                <  1:2327048/100‑1 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCttt                                                                                                <  2:513227/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                              >  2:656337/1‑100 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGc                                                                                           <  1:1010852/100‑1 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGc                                                                                           <  2:3518058/100‑1 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGc                                                                                           <  1:1627689/100‑1 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGc                                                                                           <  1:2210957/100‑1 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGc                                                                                           >  1:1421559/1‑100 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGc                                                                                           <  1:267851/100‑1 (MQ=25)
           ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGca                                                                                          <  1:2304064/100‑1 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcac                                                                                         >  2:3138508/1‑100 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                         <  2:3837350/100‑3 (MQ=25)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                         <  2:1421550/100‑3 (MQ=25)
             accaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcaca                                                                                        >  2:1987180/1‑100 (MQ=255)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                       >  2:3383271/1‑100 (MQ=255)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctg                                                                                       <  1:44808/100‑5 (MQ=21)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctg                                                                                       <  2:1253173/100‑5 (MQ=25)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                      <  2:1738966/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                      <  2:2915166/100‑1 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  1:2521746/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  1:1027226/1‑100 (MQ=255)
                     aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  2:333875/1‑98 (MQ=255)
                       cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGtt                                                                              >  2:3243117/1‑100 (MQ=255)
                        ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                             >  1:3826025/1‑100 (MQ=255)
                         gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  1:1690746/100‑1 (MQ=255)
                         gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  1:2365696/100‑1 (MQ=255)
                         gCGATATCAAGCTCCATACTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  2:1113231/100‑1 (MQ=255)
                          cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGc                                                                           <  1:893602/100‑1 (MQ=255)
                             tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                        <  2:390154/100‑1 (MQ=255)
                               tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                      >  1:615246/1‑100 (MQ=255)
                                cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                     <  2:3300660/100‑1 (MQ=255)
                                   gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTg                                                                  >  1:2099025/1‑100 (MQ=255)
                                    cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgc                                                                 >  2:1500046/1‑100 (MQ=255)
                                    cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgc                                                                 >  2:1534341/1‑100 (MQ=255)
                                        aTAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgccgcg                                                             >  2:2798514/1‑100 (MQ=255)
                                          aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                           <  2:1409917/100‑1 (MQ=255)
                                          aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                           >  2:3645862/1‑100 (MQ=255)
                                           gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                          >  2:663595/1‑100 (MQ=255)
                                            tACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGtt                                                         <  2:289458/100‑1 (MQ=255)
                                             attGATAGCATCCGCTACACCGAGCAGGACAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                        >  2:1306682/3‑100 (MQ=255)
                                             aCTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                        >  1:558190/1‑100 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                     <  2:1373228/100‑1 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                     >  1:329543/1‑100 (MQ=255)
                                                 atatCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                    <  2:1828816/100‑1 (MQ=255)
                                                  tAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACGGGATCGTTTGCCAGGAGTTGCGGCGGTCTATACa                                                    >  2:374859/1‑99 (MQ=255)
                                                   aGCATCCGCTACACCGAGCAGGGCAGTACAGCATAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACaga                                                  >  2:700310/1‑100 (MQ=255)
                                                    gCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                 >  2:3191557/1‑100 (MQ=255)
                                                    gCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                 >  2:3463728/1‑100 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagagn                                                <  2:1488803/100‑2 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                                >  2:3030349/1‑100 (MQ=255)
                                                      aTCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                               <  1:3296331/100‑1 (MQ=255)
                                                        ccGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGAc                                             <  2:2325375/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                           >  1:3108282/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACCGGc                                         >  2:3600965/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGc                                         <  1:2481097/100‑1 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                        <  1:2915484/100‑1 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                       >  1:2175590/1‑100 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                       >  2:3824667/1‑100 (MQ=255)
                                                                ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                     >  2:1687629/1‑100 (MQ=255)
                                                                 tctggtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCCCACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAcc                                    >  1:1372459/9‑100 (MQ=17)
                                                                    gCATGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                 >  2:734433/1‑100 (MQ=255)
                                                                    gCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                 >  2:1844347/1‑100 (MQ=255)
                                                                     cAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                <  1:2529134/100‑1 (MQ=255)
                                                                     cAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                <  1:1811309/100‑1 (MQ=255)
                                                                       gtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCATTTGCCAGGCGTTGCGGCGGGTCCTTCAGAGCGACAGCCGTACCGACCTg                              >  2:2548028/3‑100 (MQ=12)
                                                                       gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                              <  2:3090790/100‑1 (MQ=255)
                                                                       gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                              <  1:1083626/100‑1 (MQ=255)
                                                                        tgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  1:1918426/2‑100 (MQ=25)
                                                                        ggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  1:3044748/1‑100 (MQ=255)
                                                                         gCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                            <  1:497307/100‑1 (MQ=255)
                                                                         gCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                            >  2:2748188/1‑100 (MQ=25)
                                                                          cggtaaagcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:2677546/11‑100 (MQ=37)
                                                                          cagtaaagcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:3466815/11‑100 (MQ=37)
                                                                          cagtaaagcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:3466816/11‑100 (MQ=37)
                                                                          cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           <  1:1204107/100‑1 (MQ=255)
                                                                           aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                          >  2:784038/1‑100 (MQ=255)
                                                                           aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                          <  1:1872574/100‑1 (MQ=255)
                                                                           aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGGTTGCCAGGAGTCGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  2:120726/1‑97 (MQ=255)
                                                                            gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  2:3117332/1‑100 (MQ=255)
                                                                            gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  1:1548846/1‑100 (MQ=255)
                                                                            gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  1:1177755/1‑100 (MQ=255)
                                                                                  ccTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATcc                   >  1:532617/1‑100 (MQ=255)
                                                                                    cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCg                  <  1:3562877/100‑1 (MQ=255)
                                                                                    cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCg                  <  1:3446087/100‑1 (MQ=255)
                                                                                    cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCg                  >  1:2513182/1‑100 (MQ=255)
                                                                                      aGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGta                >  2:1076877/1‑100 (MQ=255)
                                                                                          gATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGc            >  2:1530327/1‑100 (MQ=255)
                                                                                           aTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcg           >  2:3619161/1‑100 (MQ=255)
                                                                                              gTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtcg        <  1:668318/100‑1 (MQ=255)
                                                                                              gTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtcg        >  1:3633454/1‑100 (MQ=255)
                                                                                                cTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCgtcgtc      >  1:3950115/1‑100 (MQ=255)
                                                                                                cTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCCGAATGCCCATCCGTATAGCgtcgtc      >  2:3604849/1‑100 (MQ=255)
                                                                                                 tAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtct     >  2:226996/1‑100 (MQ=255)
                                                                                                  aGTGCTTTACTGCACACCAGATCGTTTGCCGGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt    <  2:2825454/100‑1 (MQ=255)
                                                                                                    tGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttg  <  2:558187/100‑1 (MQ=255)
                                                                                                    tGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttg  <  1:3684379/100‑1 (MQ=255)
                                                                                                    tGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttg  <  1:1070194/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTG  >  NC_000913/4531019‑4531217

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: