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.

CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTG  >  NC_000913/4531019‑4531217
                                                                                                    |                                                                                                   
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGt                                                                                                     >  2:3155488/1‑100 (MQ=255)
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGt                                                                                                     >  1:1738975/1‑100 (MQ=255)
    aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCtt                                                                                                 >  1:1405999/1‑100 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCttt                                                                                                <  1:2327048/100‑1 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCttt                                                                                                <  2:513227/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                              >  2:656337/1‑100 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGc                                                                                           <  1:1010852/100‑1 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGc                                                                                           <  2:3518058/100‑1 (MQ=255)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGc                                                                                           <  1:1627689/100‑1 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGc                                                                                           <  1:2210957/100‑1 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGc                                                                                           >  1:1421559/1‑100 (MQ=25)
          gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGc                                                                                           <  1:267851/100‑1 (MQ=25)
           ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGca                                                                                          <  1:2304064/100‑1 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcac                                                                                         >  2:3138508/1‑100 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                         <  2:3837350/100‑3 (MQ=25)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                         <  2:1421550/100‑3 (MQ=25)
             accaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcaca                                                                                        >  2:1987180/1‑100 (MQ=255)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                       >  2:3383271/1‑100 (MQ=255)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccctg                                                                                       <  1:44808/100‑5 (MQ=21)
              ccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccctg                                                                                       <  2:1253173/100‑5 (MQ=25)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                      <  2:1738966/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                      <  2:2915166/100‑1 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  1:2521746/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  1:1027226/1‑100 (MQ=255)
                     aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                  >  2:333875/1‑98 (MQ=255)
                       cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGtt                                                                              >  2:3243117/1‑100 (MQ=255)
                        ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                             >  1:3826025/1‑100 (MQ=255)
                         gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  1:1690746/100‑1 (MQ=255)
                         gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  1:2365696/100‑1 (MQ=255)
                         gCGATATCAAGCTCCATACTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                            <  2:1113231/100‑1 (MQ=255)
                          cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGc                                                                           <  1:893602/100‑1 (MQ=255)
                             tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                        <  2:390154/100‑1 (MQ=255)
                               tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                      >  1:615246/1‑100 (MQ=255)
                                cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                     <  2:3300660/100‑1 (MQ=255)
                                   gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTg                                                                  >  1:2099025/1‑100 (MQ=255)
                                    cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgc                                                                 >  2:1500046/1‑100 (MQ=255)
                                    cTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgc                                                                 >  2:1534341/1‑100 (MQ=255)
                                        aTAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgccgcg                                                             >  2:2798514/1‑100 (MQ=255)
                                          aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                           <  2:1409917/100‑1 (MQ=255)
                                          aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                           >  2:3645862/1‑100 (MQ=255)
                                           gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                          >  2:663595/1‑100 (MQ=255)
                                            tACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGtt                                                         <  2:289458/100‑1 (MQ=255)
                                             attGATAGCATCCGCTACACCGAGCAGGACAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                        >  2:1306682/3‑100 (MQ=255)
                                             aCTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                        >  1:558190/1‑100 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                     <  2:1373228/100‑1 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                     >  1:329543/1‑100 (MQ=255)
                                                 atatCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                    <  2:1828816/100‑1 (MQ=255)
                                                  tAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACGGGATCGTTTGCCAGGAGTTGCGGCGGTCTATACa                                                    >  2:374859/1‑99 (MQ=255)
                                                   aGCATCCGCTACACCGAGCAGGGCAGTACAGC‑ATAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACaga                                                  >  2:700310/1‑100 (MQ=255)
                                                    gCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                 >  2:3191557/1‑100 (MQ=255)
                                                    gCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                 >  2:3463728/1‑100 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACagagn                                                <  2:1488803/100‑2 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                                >  2:3030349/1‑100 (MQ=255)
                                                      aTCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                               <  1:3296331/100‑1 (MQ=255)
                                                        ccGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGAc                                             <  2:2325375/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                           >  1:3108282/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACCGGc                                         >  2:3600965/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGc                                         <  1:2481097/100‑1 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                        <  1:2915484/100‑1 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                       >  1:2175590/1‑100 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                       >  2:3824667/1‑100 (MQ=255)
                                                                ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                     >  2:1687629/1‑100 (MQ=255)
                                                                 tctggtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCCCACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAcc                                    >  1:1372459/9‑100 (MQ=17)
                                                                    gCATGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                 >  2:734433/1‑100 (MQ=255)
                                                                    gCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                 >  2:1844347/1‑100 (MQ=255)
                                                                     cAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                <  1:2529134/100‑1 (MQ=255)
                                                                     cAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                <  1:1811309/100‑1 (MQ=255)
                                                                       gtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCATTTGCCAGGCGTTGCGGCGGGTCCTTCAGAGCGACAGCCGTACCGACCTg                              >  2:2548028/3‑100 (MQ=12)
                                                                       gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                              <  2:3090790/100‑1 (MQ=255)
                                                                       gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                              <  1:1083626/100‑1 (MQ=255)
                                                                        tgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  1:1918426/2‑100 (MQ=25)
                                                                        ggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  1:3044748/1‑100 (MQ=255)
                                                                         gCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                            <  1:497307/100‑1 (MQ=255)
                                                                         gCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                            >  2:2748188/1‑100 (MQ=25)
                                                                          cggtaaagcacTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:2677546/11‑100 (MQ=37)
                                                                          cagtaaagcacTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:3466815/11‑100 (MQ=37)
                                                                          cagtaaagcacTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           >  1:3466816/11‑100 (MQ=37)
                                                                          cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                           <  1:1204107/100‑1 (MQ=255)
                                                                           aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                          >  2:784038/1‑100 (MQ=255)
                                                                           aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                          <  1:1872574/100‑1 (MQ=255)
                                                                           aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGGTTGCCAGGAGTCGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                             >  2:120726/1‑97 (MQ=255)
                                                                            gTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  2:3117332/1‑100 (MQ=255)
                                                                            gTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  1:1548846/1‑100 (MQ=255)
                                                                            gTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                         >  1:1177755/1‑100 (MQ=255)
                                                                                  c‑cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATcc                   >  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)
                                                                                                    |                                                                                                   
CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTG  >  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: ‑