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.

ACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTT  >  NC_000913/4531023‑4531221
                                                                                                    |                                                                                                   
aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCtt                                                                                                     <  1:1486589/100‑1 (MQ=255)
 ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATGCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCttt                                                                                                    <  1:498610/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                                  <  1:928816/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                                  <  2:2791597/100‑1 (MQ=255)
   gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                                  >  1:1620433/1‑100 (MQ=255)
     aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTg                                                                                                <  2:130238/100‑1 (MQ=255)
      gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGc                                                                                               <  2:1233264/100‑1 (MQ=255)
      gCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGc                                                                                               <  2:2535119/100‑1 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGGTCGTCTAGTGCTTTACTGca                                                                                              <  1:2149729/100‑1 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGca                                                                                              >  1:2122812/1‑100 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATAGTCTAGTGCTTTACTGca                                                                                              >  2:615101/1‑100 (MQ=255)
       ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGcc                                                                                              <  1:2607509/100‑2 (MQ=25)
        caccaGATAAACTGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                             <  1:180129/100‑3 (MQ=25)
        caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcac                                                                                             >  2:1735979/1‑100 (MQ=255)
        caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                             <  1:2135366/100‑3 (MQ=25)
         accaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcaca                                                                                            >  2:470021/1‑100 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:285817/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:3715665/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:2190617/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  1:2349949/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                          <  2:3720674/100‑1 (MQ=255)
           caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctgc                                                                                          <  1:2685257/100‑6 (MQ=21)
               tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                      >  1:711559/1‑100 (MQ=255)
               tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                      <  1:1244724/100‑1 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCg                                                                                    <  2:222705/98‑1 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGtt                                                                                  >  2:2694003/1‑100 (MQ=255)
                   cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCCAGTGCTTTACTGCACACCAGATCGtt                                                                                  >  2:3352099/1‑100 (MQ=255)
                    ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                 <  1:2298796/100‑1 (MQ=255)
                     gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTg                                                                                <  2:2910922/100‑1 (MQ=255)
                      cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCTCACCAGATCGTTTGc                                                                               >  1:3757747/1‑100 (MQ=255)
                      cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGc                                                                               <  1:2901355/100‑1 (MQ=255)
                       gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGTTTTACTGCACACCAGATCGTTTGcc                                                                              <  1:531045/100‑1 (MQ=255)
                        atatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCa                                                                             <  1:1735989/100‑1 (MQ=255)
                         tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                            <  2:95938/100‑1 (MQ=255)
                         tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                            <  2:2455449/100‑1 (MQ=255)
                          atCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                           <  1:1073764/100‑1 (MQ=255)
                           tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                          >  1:1562396/1‑100 (MQ=255)
                            cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                           >  2:694045/1‑98 (MQ=255)
                            cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                         <  1:2865909/100‑1 (MQ=255)
                             aaGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGGCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                         <  1:2748377/99‑1 (MQ=255)
                                 tCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcg                                                                    >  1:2417977/1‑100 (MQ=255)
                                  ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                                   >  1:1059420/1‑100 (MQ=255)
                                       gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                              <  2:1891345/100‑1 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         >  1:2958196/1‑100 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         >  2:632329/1‑100 (MQ=255)
                                            gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                         <  1:424646/100‑1 (MQ=255)
                                             aTAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACa                                                        <  1:1505453/100‑1 (MQ=255)
                                                  aTCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                                   <  2:3249873/100‑1 (MQ=255)
                                                  aTCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                                   >  2:572840/1‑100 (MQ=255)
                                                   tCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                                  >  2:387445/1‑100 (MQ=255)
                                                   tCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                                  >  1:3602895/1‑100 (MQ=255)
                                                     cGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGAc                                                 >  1:1684221/1‑99 (MQ=255)
                                                     cGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACa                                                <  2:1172503/100‑1 (MQ=255)
                                                      gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                               <  2:1888651/100‑1 (MQ=255)
                                                      gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGCTGCGGCGGGTTATACAGAGTGACAg                                               >  2:3535651/1‑100 (MQ=255)
                                                       cTACACCGAGCAGGGCAGTACAGCCTAGACGATCGGCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                              <  2:1442356/100‑1 (MQ=255)
                                                        tACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                              >  2:2259440/1‑99 (MQ=255)
                                                          cacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                           >  1:2650610/1‑100 (MQ=255)
                                                           nccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTa                                          >  1:3801490/2‑100 (MQ=255)
                                                            ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                         >  2:812214/1‑100 (MQ=255)
                                                            ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                         <  1:3034302/100‑1 (MQ=255)
                                                             cGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCCGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGGAcc                                        >  2:2139736/1‑100 (MQ=255)
                                                                ggtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                     >  1:2707091/6‑100 (MQ=21)
                                                                gCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                     <  2:1939826/100‑1 (MQ=255)
                                                                 gtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    >  1:1232006/5‑100 (MQ=21)
                                                                 gtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    >  2:2126607/5‑100 (MQ=21)
                                                                 cAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGATTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                                    <  2:2030091/100‑1 (MQ=255)
                                                                  tgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGATATACAGAGTGACAGGCGTACCGACCt                                   >  2:3040626/4‑100 (MQ=21)
                                                                   gtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCTAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  >  1:177976/3‑100 (MQ=25)
                                                                   gtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  >  1:228207/3‑100 (MQ=25)
                                                                   gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                  <  2:3504674/100‑1 (MQ=255)
                                                                   gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGTTATACCGCGTGACAGGCGTACCGACCt                                   >  2:322538/1‑99 (MQ=255)
                                                                    tgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGATGCGGCGGGTTATACAGAGTGACAGGCGTACCCACCTGa                                 >  1:541797/2‑100 (MQ=25)
                                                                    ggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                                 >  1:2813129/1‑100 (MQ=255)
                                                                     gCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                <  1:3121130/100‑1 (MQ=255)
                                                                     gCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                >  1:3688124/1‑100 (MQ=25)
                                                                      cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                               <  2:856930/100‑1 (MQ=255)
                                                                      cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                               <  1:768831/100‑1 (MQ=255)
                                                                         tACAGCCTAGACGATCGTCTAGTGCTTTACTGCGCACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAGAGAGTGACAGGCGTACCGACCTGAATGcc                            >  2:2301004/1‑100 (MQ=255)
                                                                         tACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGcc                            <  1:2263247/100‑1 (MQ=255)
                                                                           cAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                          <  1:1316168/100‑1 (MQ=255)
                                                                           cAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCa                          <  2:1014038/100‑1 (MQ=255)
                                                                             gcacTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCAt                         >  2:2205632/4‑100 (MQ=255)
                                                                             gcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                        >  1:3067597/4‑100 (MQ=38)
                                                                             gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGCGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCAt                         >  2:2611406/1‑99 (MQ=255)
                                                                             gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                        <  1:280683/100‑1 (MQ=255)
                                                                              ccTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTAGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATcc                       >  2:3027608/1‑100 (MQ=255)
                                                                                 tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGt                     <  2:3426195/100‑1 (MQ=255)
                                                                                   gACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGtata                  >  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)
                                                                                                    |                                                                                                   
ACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTT  >  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: