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% 455.9 / 83.4 163coding (534/1122 nt)sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (62/75);  new base . (18/7);  total (80/83)
Fisher's exact test for biased strand distribution p-value = 1.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 164 total aligned reads displayed.

CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTT  >  NC_000913/4531019‑4531216
                                                                                                    |                                                                                                  
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGt                                                                                                    >  1:217330/1‑100 (MQ=255)
 tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTg                                                                                                   >  1:1492631/1‑100 (MQ=255)
  gTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGc                                                                                                  >  2:2190603/1‑100 (MQ=255)
   tACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCt                                                                                                 >  2:500743/1‑100 (MQ=255)
    aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCGGCTACACCGAGCAGATCAGTACAGCGTAGACGATCGTCTAGTGCt                                                                                                 >  2:758571/1‑99 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATGCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCttt                                                                                               <  1:498610/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                             <  2:840610/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                             <  1:151662/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                             <  2:2791597/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTAc                                                                                             >  1:1620433/1‑100 (MQ=255)
        aaGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACt                                                                                            <  2:1513340/100‑1 (MQ=255)
         aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTg                                                                                           <  2:130238/100‑1 (MQ=255)
           ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGcc                                                                                         <  1:2607509/100‑2 (MQ=25)
            caccaGATAAACTGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                        <  1:180129/100‑3 (MQ=25)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcac                                                                                        <  1:2160571/100‑1 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                        <  1:2135366/100‑3 (MQ=25)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:2190617/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:2349949/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  2:1071246/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:285817/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:3715665/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctgc                                                                                     <  1:2685257/100‑6 (MQ=21)
                aGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCa                                                                                    >  2:2980446/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGTTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                 >  1:3223737/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                 <  1:3747672/100‑1 (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)
                          cGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGc                                                                          <  1:2901355/100‑1 (MQ=255)
                           gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGTTTTACTGCACACCAGATCGTTTGcc                                                                         <  1:531045/100‑1 (MQ=255)
                            atatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCa                                                                        <  1:1735989/100‑1 (MQ=255)
                             tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                       <  2:2455449/100‑1 (MQ=255)
                             tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                       <  2:95938/100‑1 (MQ=255)
                              atCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                      <  1:1073764/100‑1 (MQ=255)
                               tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                     >  2:3480726/1‑100 (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)
                                 aaGCTCCATAGTACTAATAGCATCCGCTACACCGAGCAGTGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                       <  2:125465/96‑1 (MQ=255)
                                   gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTg                                                                 >  1:3809030/1‑100 (MQ=255)
                                     tCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcg                                                               >  1:2417977/1‑100 (MQ=255)
                                      ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                              >  1:1059420/1‑100 (MQ=255)
                                      ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                              >  1:1913575/1‑100 (MQ=255)
                                         tAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGcggcgg                                                           <  2:1862183/100‑1 (MQ=255)
                                           gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                         <  2:1891345/100‑1 (MQ=255)
                                              cTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTtat                                                      <  2:1210635/100‑1 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                    >  2:632329/1‑100 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                    >  1:2958196/1‑100 (MQ=255)
                                                    gCATCCGCCACCTCGAGCAGGGGAGTAAAGACAGGACGAGCGTCTAGTGCTTTACTGCACACCACATCCGATGCCAGGAGGGGCAGCGGATtaca                                                     >  2:2504307/1‑93 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                               <  2:724471/100‑1 (MQ=255)
                                                      aTCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                              >  2:572840/1‑100 (MQ=255)
                                                       tCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                             >  2:387445/1‑100 (MQ=255)
                                                         cGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACa                                           <  2:3085623/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                          <  2:1888651/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGCTGCGGCGGGTTATACAGAGTGACAg                                          >  2:3535651/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                         >  2:2259440/1‑99 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                       <  1:798815/100‑1 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                       <  1:1317159/100‑1 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                      >  1:2650610/1‑100 (MQ=255)
                                                               nccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTa                                     >  1:3801490/2‑100 (MQ=255)
                                                                ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    <  1:3034302/100‑1 (MQ=255)
                                                                atctggtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    >  1:1772890/10‑100 (MQ=18)
                                                                atctggtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    >  1:1161390/10‑100 (MQ=18)
                                                                 cGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCCGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGGAcc                                   >  2:2139736/1‑100 (MQ=255)
                                                                  gAGCAGGTCATTACAGCATAGAGGACGGTCTCGTGGTGTACTGCAAACCAGGTCGTTGGCGAGGAGGGGCGGCGGTTCATACAGAGTGACAGGCGTAcc                                   <  2:3595163/99‑1 (MQ=255)
                                                                    ggtgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                >  1:2076357/6‑100 (MQ=21)
                                                                    gCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                <  2:1939826/100‑1 (MQ=255)
                                                                     cAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                               <  1:709100/100‑1 (MQ=255)
                                                                      tgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                              >  2:335304/4‑100 (MQ=25)
                                                                      tgtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGATATACAGAGTGACAGGCGTACCGACCt                              >  2:3040626/4‑100 (MQ=21)
                                                                      aGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                              <  1:2137616/100‑1 (MQ=255)
                                                                       gtgCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                             >  1:228207/3‑100 (MQ=25)
                                                                        ggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGa                            >  1:2813129/1‑100 (MQ=255)
                                                                         gCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                           <  1:3121130/100‑1 (MQ=255)
                                                                         gCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                           >  1:3688124/1‑100 (MQ=25)
                                                                          cagtaaagcacTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          >  1:537452/11‑100 (MQ=37)
                                                                          cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          <  2:856930/100‑1 (MQ=255)
                                                                          cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          <  1:768831/100‑1 (MQ=255)
                                                                           aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                         <  2:2549346/100‑1 (MQ=255)
                                                                           aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                         <  2:203850/100‑1 (MQ=255)
                                                                            gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                        >  2:3732993/1‑100 (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                   >  2:599897/1‑100 (MQ=255)
                                                                                 gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                   <  1:280683/100‑1 (MQ=255)
                                                                                 gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGTACCGACCTGCATGCCCATc                   >  2:3239660/1‑100 (MQ=255)
                                                                                         cGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAg            >  2:2528609/1‑100 (MQ=255)
                                                                                         cGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAg            >  2:1706740/1‑100 (MQ=255)
                                                                                            tCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgt         >  1:1145598/1‑100 (MQ=255)
                                                                                              gTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtcg       >  2:1544606/1‑100 (MQ=255)
                                                                                                   gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtcttt  >  1:3770092/1‑100 (MQ=255)
                                                                                                   gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt   >  2:1432595/1‑99 (MQ=255)
                                                                                                    |                                                                                                  
CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTT  >  NC_000913/4531019‑4531216

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: