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.

CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTT  >  NC_000913/4531019‑4531216
                                                                                                    |                                                                                                  
cTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGt                                                                                                    >  1:217330/1‑100 (MQ=255)
 tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTg                                                                                                   >  1:1492631/1‑100 (MQ=255)
  gTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGc                                                                                                  >  2:2190603/1‑100 (MQ=255)
   tACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCt                                                                                                 >  2:500743/1‑100 (MQ=255)
    aCCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCGGCTACACCGAGCAGATCAGTACAGC‑GTAGACGATCGTCTAGTGCt                                                                                                 >  2:758571/1‑99 (MQ=255)
     ccGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATGCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCttt                                                                                               <  1:498610/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                             <  2:840610/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                             <  1:151662/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                             <  2:2791597/100‑1 (MQ=255)
       gAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTAc                                                                                             >  1:1620433/1‑100 (MQ=255)
        aaGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACt                                                                                            <  2:1513340/100‑1 (MQ=255)
         aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTg                                                                                           <  2:130238/100‑1 (MQ=255)
           ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGcc                                                                                         <  1:2607509/100‑2 (MQ=25)
            caccaGATAAACTGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                        <  1:180129/100‑3 (MQ=25)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcac                                                                                        <  1:2160571/100‑1 (MQ=255)
            caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                        <  1:2135366/100‑3 (MQ=25)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:2190617/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:2349949/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  2:1071246/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:285817/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                     <  1:3715665/100‑1 (MQ=255)
               caGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccctgc                                                                                     <  1:2685257/100‑6 (MQ=21)
                aGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCa                                                                                    >  2:2980446/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGTTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                 >  1:3223737/1‑100 (MQ=255)
                   tAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAt                                                                                 <  1:3747672/100‑1 (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)
                          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:2455449/100‑1 (MQ=255)
                             tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                       <  2:95938/100‑1 (MQ=255)
                              atCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAgg                                                                      <  1:1073764/100‑1 (MQ=255)
                               tCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGa                                                                     >  2:3480726/1‑100 (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)
                                 aaGCTCCATAGTACTAATAGCATCCGCTACACCGAGCAGTGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                       <  2:125465/96‑1 (MQ=255)
                                   gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTg                                                                 >  1:3809030/1‑100 (MQ=255)
                                     tCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcg                                                               >  1:2417977/1‑100 (MQ=255)
                                      ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                              >  1:1059420/1‑100 (MQ=255)
                                      ccATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                              >  1:1913575/1‑100 (MQ=255)
                                         tAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGcggcgg                                                           <  2:1862183/100‑1 (MQ=255)
                                           gTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGt                                                         <  2:1891345/100‑1 (MQ=255)
                                              cTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTtat                                                      <  2:1210635/100‑1 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                    >  2:632329/1‑100 (MQ=255)
                                                gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATAc                                                    >  1:2958196/1‑100 (MQ=255)
                                                    gCATCCGCCACCTCGAGCAGGGGAGTAAAGA‑CAGGACGAGCGTCTAGTGCTTTACTGCACACCACATCCGATGCCAGGAGGGGCAGCGGATtaca                                                     >  2:2504307/1‑93 (MQ=255)
                                                     cATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                               <  2:724471/100‑1 (MQ=255)
                                                      aTCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTg                                              >  2:572840/1‑100 (MQ=255)
                                                       tCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGa                                             >  2:387445/1‑100 (MQ=255)
                                                         cGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACa                                           <  2:3085623/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                          <  2:1888651/100‑1 (MQ=255)
                                                          gCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGCTGCGGCGGGTTATACAGAGTGACAg                                          >  2:3535651/1‑100 (MQ=255)
                                                            tACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                         >  2:2259440/1‑99 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                       <  1:798815/100‑1 (MQ=255)
                                                             acacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                       <  1:1317159/100‑1 (MQ=255)
                                                              cacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                      >  1:2650610/1‑100 (MQ=255)
                                                               nccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTa                                     >  1:3801490/2‑100 (MQ=255)
                                                                ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    <  1:3034302/100‑1 (MQ=255)
                                                                atctggtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    >  1:1772890/10‑100 (MQ=18)
                                                                atctggtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTAc                                    >  1:1161390/10‑100 (MQ=18)
                                                                 cGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCCGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGGAcc                                   >  2:2139736/1‑100 (MQ=255)
                                                                  gAGCAGGTCATTACAGC‑ATAGAGGACGGTCTCGTGGTGTACTGCAAACCAGGTCGTTGGCGAGGAGGGGCGGCGGTTCATACAGAGTGACAGGCGTAcc                                   <  2:3595163/99‑1 (MQ=255)
                                                                    ggtgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                >  1:2076357/6‑100 (MQ=21)
                                                                    gCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAc                                <  2:1939826/100‑1 (MQ=255)
                                                                     cAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGAcc                               <  1:709100/100‑1 (MQ=255)
                                                                      tgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                              >  2:335304/4‑100 (MQ=25)
                                                                      tgtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTCGCGGCGGGATATACAGAGTGACAGGCGTACCGACCt                              >  2:3040626/4‑100 (MQ=21)
                                                                      aGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                              <  1:2137616/100‑1 (MQ=255)
                                                                       gtgCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                             >  1:228207/3‑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)
                                                                          cagtaaagcacTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          >  1:537452/11‑100 (MQ=37)
                                                                          cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          <  2:856930/100‑1 (MQ=255)
                                                                          cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAAt                          <  1:768831/100‑1 (MQ=255)
                                                                           aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                         <  2:2549346/100‑1 (MQ=255)
                                                                           aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                         <  2:203850/100‑1 (MQ=255)
                                                                            gTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGc                        >  2:3732993/1‑100 (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                   >  2:599897/1‑100 (MQ=255)
                                                                                 gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATc                   <  1:280683/100‑1 (MQ=255)
                                                                                 gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGGGACAGGCGTACCGACCTGCATGCCCATc                   >  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)
                                                                                                    |                                                                                                  
CTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTT  >  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: ‑