Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,531,122 T→G 15.9% K177T (AAA→ACA)  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,1220TG15.9% 434.3 / 67.3 161K177T (AAA→ACA) sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (62/73);  new base G (18/8);  total (80/81)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 168 total aligned reads displayed.

ACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTT  >  NC_000913/4531023‑4531221
                                                                                                    |                                                                                                   
aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCtt                                                                                                     <  1:1486589/100‑1 (MQ=255)
 ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATGCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCttt                                                                                                    <  1:498610/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                                  <  1:928816/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                                  <  2:2791597/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                                  >  1:1620433/1‑100 (MQ=255)
     aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTg                                                                                                <  2:130238/100‑1 (MQ=255)
      gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGc                                                                                               <  2:1233264/100‑1 (MQ=255)
      gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGc                                                                                               <  2:2535119/100‑1 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGGTCGTCTAGTGCTTTACTGca                                                                                              <  1:2149729/100‑1 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGca                                                                                              >  1:2122812/1‑100 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATAGTCTAGTGCTTTACTGca                                                                                              >  2:615101/1‑100 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGcc                                                                                              <  1:2607509/100‑2 (MQ=25)
        caccaGATAAACTGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                             <  1:180129/100‑3 (MQ=25)
        caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcac                                                                                             >  2:1735979/1‑100 (MQ=255)
        caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                             <  1:2135366/100‑3 (MQ=25)
         accaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcaca                                                                                            >  2:470021/1‑100 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:285817/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:3715665/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:2190617/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:2349949/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  2:3720674/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccctgc                                                                                          <  1:2685257/100‑6 (MQ=21)
               tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                      >  1:711559/1‑100 (MQ=255)
               tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                      <  1:1244724/100‑1 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCg                                                                                    <  2:222705/98‑1 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGtt                                                                                  >  2:2694003/1‑100 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCCAGTGCTTTACTGCACACCAGATCGtt                                                                                  >  2:3352099/1‑100 (MQ=255)
                    ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                 <  1:2298796/100‑1 (MQ=255)
                     gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                                <  2:2910922/100‑1 (MQ=255)
                      cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCTCACCAGATCGTTTGc                                                                               >  1:3757747/1‑100 (MQ=255)
                      cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGc                                                                               <  1:2901355/100‑1 (MQ=255)
                       gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGTTTTACTGCACACCAGATCGTTTGcc                                                                              <  1:531045/100‑1 (MQ=255)
                        atatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCa                                                                             <  1:1735989/100‑1 (MQ=255)
                         tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                            <  2:95938/100‑1 (MQ=255)
                         tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                            <  2:2455449/100‑1 (MQ=255)
                          atCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                           <  1:1073764/100‑1 (MQ=255)
                           tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                          >  1:1562396/1‑100 (MQ=255)
                            cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                           >  2:694045/1‑98 (MQ=255)
                            cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                         <  1:2865909/100‑1 (MQ=255)
                             aaGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGGCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                         <  1:2748377/99‑1 (MQ=255)
                                 tCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcg                                                                    >  1:2417977/1‑100 (MQ=255)
                                  ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                                   >  1:1059420/1‑100 (MQ=255)
                                       gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                              <  2:1891345/100‑1 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         >  1:2958196/1‑100 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         >  2:632329/1‑100 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         <  1:424646/100‑1 (MQ=255)
                                             aTAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                        <  1:1505453/100‑1 (MQ=255)
                                                  aTCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                                   <  2:3249873/100‑1 (MQ=255)
                                                  aTCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                                   >  2:572840/1‑100 (MQ=255)
                                                   tCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                                  >  2:387445/1‑100 (MQ=255)
                                                   tCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                                  >  1:3602895/1‑100 (MQ=255)
                                                     cGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGAc                                                 >  1:1684221/1‑99 (MQ=255)
                                                     cGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACa                                                <  2:1172503/100‑1 (MQ=255)
                                                      gCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                               <  2:1888651/100‑1 (MQ=255)
                                                      gCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGCTGCGGCGGGTTATACAGAGTGACAg                                               >  2:3535651/1‑100 (MQ=255)
                                                       cTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGGCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                              <  2:1442356/100‑1 (MQ=255)
                                                        tACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                              >  2:2259440/1‑99 (MQ=255)
                                                          cacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                           >  1:2650610/1‑100 (MQ=255)
                                                           nccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTa                                          >  1:3801490/2‑100 (MQ=255)
                                                            ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                         >  2:812214/1‑100 (MQ=255)
                                                            ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                         <  1:3034302/100‑1 (MQ=255)
                                                             cGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCCGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGGAcc                                        >  2:2139736/1‑100 (MQ=255)
                                                                ggtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                     >  1:2707091/6‑100 (MQ=21)
                                                                gCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                     <  2:1939826/100‑1 (MQ=255)
                                                                 gtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    >  1:1232006/5‑100 (MQ=21)
                                                                 gtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    >  2:2126607/5‑100 (MQ=21)
                                                                 cAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGATTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    <  2:2030091/100‑1 (MQ=255)
                                                                  tgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGATATACAGAGTGACAGGCGTACCGACCt                                   >  2:3040626/4‑100 (MQ=21)
                                                                   gtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCTAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  >  1:177976/3‑100 (MQ=25)
                                                                   gtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  >  1:228207/3‑100 (MQ=25)
                                                                   gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  <  2:3504674/100‑1 (MQ=255)
                                                                   gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGTTATACCGCGTGACAGGCGTACCGACCt                                   >  2:322538/1‑99 (MQ=255)
                                                                    tgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGATGCGGCGGGTTATACAGAGTGACAGGCGTACCCACCTGa                                 >  1:541797/2‑100 (MQ=25)
                                                                    ggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                                 >  1:2813129/1‑100 (MQ=255)
                                                                     gCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                <  1:3121130/100‑1 (MQ=255)
                                                                     gCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                >  1:3688124/1‑100 (MQ=25)
                                                                      cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                               <  2:856930/100‑1 (MQ=255)
                                                                      cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                               <  1:768831/100‑1 (MQ=255)
                                                                         tACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCGCACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAGAGAGTGACAGGCGTACCGACCTGAATGcc                            >  2:2301004/1‑100 (MQ=255)
                                                                         tACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGcc                            <  1:2263247/100‑1 (MQ=255)
                                                                           cAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                          <  1:1316168/100‑1 (MQ=255)
                                                                           cAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                          <  2:1014038/100‑1 (MQ=255)
                                                                             gcacTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCAt                         >  2:2205632/4‑100 (MQ=255)
                                                                             gcacTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                        >  1:3067597/4‑100 (MQ=38)
                                                                             gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGCGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCAt                         >  2:2611406/1‑99 (MQ=255)
                                                                             gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                        <  1:280683/100‑1 (MQ=255)
                                                                              c‑cTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTAGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATcc                       >  2:3027608/1‑100 (MQ=255)
                                                                                 tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGt                     <  2:3426195/100‑1 (MQ=255)
                                                                                   gACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGtata                  >  1:1958605/1‑100 (MQ=255)
                                                                                     cGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAg                 >  2:2528609/1‑100 (MQ=255)
                                                                                     cGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAg                 >  2:1706740/1‑100 (MQ=255)
                                                                                       aTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcg               >  2:743783/1‑100 (MQ=255)
                                                                                        tCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgt              >  1:1145598/1‑100 (MQ=255)
                                                                                          gTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtcg            >  2:1544606/1‑100 (MQ=255)
                                                                                              aGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt        <  1:2171730/100‑1 (MQ=255)
                                                                                               gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtcttt       >  1:3770092/1‑100 (MQ=255)
                                                                                               gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt        >  2:1432595/1‑99 (MQ=255)
                                                                                                  cTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtc    >  1:2556620/1‑100 (MQ=255)
                                                                                                   tttACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtct   >  1:989109/1‑100 (MQ=255)
                                                                                                    ttACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtctt  >  2:1241883/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
ACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTT  >  NC_000913/4531023‑4531221

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