Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,531,122 T→G 9.5% 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,1220TG 9.5% 298.8 / 17.7 111K177T (AAA→ACA) sgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base T (55/45);  new base G (4/7);  total (59/52)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 112 total aligned reads displayed.

CACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGC  >  NC_000913/4530976‑4531271
                                                                                                                                                   |                                                                                                                                                     
 acGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACATC‑CTAGACGATCGTCTAGTGCTTTACt                                                                                                                                                   >  2:898065/1‑150 (MQ=255)
   ggCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGc                                                                                                                                                 <  1:1317050/150‑1 (MQ=39)
    gcttAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGcc                                                                                                                                                <  1:139041/147‑2 (MQ=38)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcac                                                                                                                                               >  1:819732/1‑150 (MQ=255)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGcac                                                                                                                                               <  1:465569/150‑1 (MQ=255)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                                                                               <  2:1425080/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                                                                               <  2:1311781/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                                                                               <  2:890641/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccc                                                                                                                                               <  1:1711898/150‑3 (MQ=39)
       tAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                                                                             >  2:1488342/1‑150 (MQ=255)
       tAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATAGCACGCTCCATAGTACTGATACCATCCGCTACACCGCGCATGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                                                                             <  1:951731/150‑1 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            <  2:1276787/150‑1 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            >  1:1251712/1‑150 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            >  1:1005827/1‑150 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGccctgc                                                                                                                                            <  1:111596/150‑6 (MQ=38)
          tACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAg                                                                                                                                          >  1:1091550/1‑150 (MQ=255)
            cGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCTTCTAGTGCTGTACTGCACACCAGag                                                                                                                                        >  1:664427/1‑149 (MQ=255)
             gggAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATc                                                                                                                                       <  2:1413652/150‑1 (MQ=255)
                 aCAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                                                                   <  2:819732/150‑1 (MQ=255)
                 aCAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                                                                   <  1:1343670/150‑1 (MQ=255)
                    aTGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGcc                                                                                                                                >  1:405543/1‑150 (MQ=255)
                      gCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                                                                              <  1:693119/150‑1 (MQ=255)
                         gcgGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                                                                           <  2:891575/150‑1 (MQ=255)
                           ggATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGtt                                                                                                                         >  2:670932/1‑150 (MQ=255)
                               atTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                                                                                     >  1:562460/1‑150 (MQ=255)
                                ttAACTTCTTCCTGTACCGACGCCACCAGCTAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggc                                                                                                                    <  1:1246994/150‑1 (MQ=255)
                                 tAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggcg                                                                                                                   <  1:1194263/150‑1 (MQ=255)
                                 tAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggcg                                                                                                                   >  2:313953/1‑150 (MQ=255)
                                  ataTTCTTCCTGTACCGAAGCCACCAGATAAACGGCTATATCAAGCTCCATAGTACTGTTAGCATCCGCTACACCGAGCAGGGCAGTACAGG‑CTAGAGGATCGTGTAGTGCTTTACTGCAGACAAGATCGTTTGCCAGGAGTTGcggcgg                                                                                                                  <  1:1229372/148‑1 (MQ=255)
                                   aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                                 <  1:1604621/150‑1 (MQ=255)
                                   aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                                 <  2:583985/150‑1 (MQ=255)
                                     ttcttcCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGtt                                                                                                               <  2:230809/150‑1 (MQ=255)
                                      tcttcCTGTACCGAAGCCACCAGATAAACGGAGATATCAAGCTCCATAGTACTTAAAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                                                                              <  2:1065657/150‑1 (MQ=255)
                                           cTGTACCGAAGCCACCAGATAAACGGCGATATCACGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATActg                                                                                                         >  2:323785/1‑148 (MQ=255)
                                            tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACaga                                                                                                        >  2:1588974/1‑150 (MQ=255)
                                             gTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACCGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                                                                       >  2:137042/1‑150 (MQ=255)
                                              tACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                                                                                      <  2:619456/150‑1 (MQ=255)
                                               aCCGAATCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACTCCGAGCATAGCAGTACAGC‑CTAGACTATCGTCTAGTGCTTTACTGCACACCATATCGTTTGCAAGGAGATGCGTCGGGTTATACagagag                                                                                                     >  2:21311/1‑148 (MQ=255)
                                                   aaGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGGTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                                                                                 <  1:1609608/150‑1 (MQ=255)
                                                    aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                                                                                <  1:662888/150‑1 (MQ=255)
                                                      ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                                                                              <  2:1562445/150‑1 (MQ=255)
                                                       caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                                                                             >  2:200768/1‑150 (MQ=255)
                                                           aGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACATAGTGACAGGCGTACCg                                                                                         >  2:1040785/1‑150 (MQ=255)
                                                               aaaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                     <  1:401257/150‑1 (MQ=255)
                                                               aaaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                     >  2:928810/1‑150 (MQ=255)
                                                                aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCATGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  1:1311263/150‑1 (MQ=255)
                                                                aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  1:1308889/150‑1 (MQ=255)
                                                                aaCGGCGAAATCCAGCTCCATACTACTCATAGCAGCCGCTACACCGAGCAGGGCAGTACTGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  2:1073980/150‑1 (MQ=255)
                                                                  cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                                                                  >  1:1154169/1‑150 (MQ=255)
                                                                  cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                                                                  >  1:1420365/1‑150 (MQ=255)
                                                                   ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTAAAACGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCATATCGTTTGCCAGGAGTTGATGCGGGTTATACAGAGTGACAGGCGTACAGACCTGAAt                                                                                 >  2:720771/1‑150 (MQ=255)
                                                                    gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                                                                                <  1:922216/150‑1 (MQ=255)
                                                                      gATATCAAGCTCCATAGTACTGATAGCATCCGCTACCCCGAGCAGGGCAGTACAGC‑CTAGACGCTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGCCCGGAATGcc                                                                              >  2:1020941/1‑150 (MQ=255)
                                                                      gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGcc                                                                              >  1:885111/1‑150 (MQ=255)
                                                                        tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGTCAGGCGTACCGACCTGAATGCCCa                                                                            >  2:887090/1‑150 (MQ=255)
                                                                        tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGACTGCCCa                                                                            >  1:1182471/1‑150 (MQ=255)
                                                                           cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTGTGCCAGGAGTTGCGGCGGGTTATACCGAGTGACAGGCGTACCGCCCTGAATGCCCATcc                                                                         >  2:1190051/1‑150 (MQ=255)
                                                                            aaGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCg                                                                        <  1:1352767/150‑1 (MQ=255)
                                                                            aaGCTCCAAAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCCTCCg                                                                        >  1:1128272/1‑150 (MQ=255)
                                                                              gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGta                                                                      >  2:1078956/1‑150 (MQ=255)
                                                                                   aTAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcg                                                                 <  1:1488342/150‑1 (MQ=255)
                                                                                     aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtc                                                               <  1:200768/150‑1 (MQ=255)
                                                                                          tGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt                                                          <  1:1443736/150‑1 (MQ=255)
                                                                                           gATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGATTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGGATAGCGTCgtcttt                                                         >  1:1403237/1‑150 (MQ=255)
                                                                                            aTAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtcgttg                                                        >  1:1623601/1‑150 (MQ=255)
                                                                                              aGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtc                                                      <  2:927976/150‑1 (MQ=255)
                                                                                                cATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtctt                                                    <  1:1286510/150‑1 (MQ=255)
                                                                                                cATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtctt                                                    <  1:1726086/150‑1 (MQ=255)
                                                                                                  tCCGCTAAACCGAGCAGGGCAGTACAGC‑ATAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCACATCCGTATAGCGTCgtcgttgtagttt                                                  >  1:603567/1‑145 (MQ=255)
                                                                                                   ccGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTgg                                                 >  1:1603822/1‑150 (MQ=255)
                                                                                                    cGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTGTGCCTGGAGTTGCGGCGGTTTATACAGAGTGACAGGCGTACCGACCTGAAGGCCCATCCGTATAGCGTCGTCtttggctttggc                                                >  2:863938/1‑150 (MQ=255)
                                                                                                         cacCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATCCCCCTCCGTATAGCGTCGTGTTTGTCTTTGGCACCaa                                           >  2:832650/1‑150 (MQ=255)
                                                                                                           ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAata                                         >  2:39475/1‑150 (MQ=255)
                                                                                                           ccGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAata                                         <  1:1078956/150‑1 (MQ=255)
                                                                                                                cAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAAt                                    >  2:1304604/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATcc                                  >  2:1721334/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATcc                                  >  2:841458/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCAACAATATCAAAcc                                  >  1:1417257/1‑150 (MQ=255)
                                                                                                                     cAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcaca                               <  1:175823/150‑1 (MQ=255)
                                                                                                                     cAGTAAAGCACTAGACGATCGTCTAG‑GCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcaca                               >  1:1251166/1‑150 (MQ=39)
                                                                                                                      aGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAg                              <  2:1498321/150‑1 (MQ=255)
                                                                                                                       gTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAAAATCAATCAACAGt                             >  2:1583073/1‑150 (MQ=255)
                                                                                                                            gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                        >  2:1147990/1‑150 (MQ=255)
                                                                                                                            gC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGAGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                        <  1:898065/150‑1 (MQ=255)
                                                                                                                                tAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGAc                     >  2:932123/1‑150 (MQ=255)
                                                                                                                                  gACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAg                   >  1:996837/1‑150 (MQ=255)
                                                                                                                                        cGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTg             <  1:670932/150‑1 (MQ=255)
                                                                                                                                        cGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTATCGACAGAAGGTg             <  2:1561204/150‑1 (MQ=255)
                                                                                                                                          tCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGAc           <  1:39475/150‑1 (MQ=255)
                                                                                                                                          tCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGAc           <  1:1304604/150‑1 (MQ=255)
                                                                                                                                          tCTAGTGCTTTACTGCACACCAGATCGTTTGCCACGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGAc           <  2:1091550/150‑1 (MQ=255)
                                                                                                                                           cTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACt          >  1:441201/1‑150 (MQ=255)
                                                                                                                                            tAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTGTTTGGCCCCAATTTGAATCCCCAGTTTGTCGACAGAAGGTGGCtg         >  2:357435/1‑150 (MQ=255)
                                                                                                                                            tAGTGCTTTACTGCACACCAGATCGTTTGACAGGAGTTCCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGACCATCCGTATAGCGTCGTCTTTGTCTTTGGCTCCAAAATCAATCCACAGATTGTCGACATAAGTTGACtt         >  1:269158/1‑149 (MQ=255)
                                                                                                                                             aGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgc        <  2:1005827/150‑1 (MQ=255)
                                                                                                                                              gTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGCCCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgcg       >  2:46621/1‑150 (MQ=255)
                                                                                                                                              aggCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGCGTTAAACAGAGTGACAGGCGTACCGACCAGAATGACCATCCGTATAGCGTCGTCTTTTTCTTTGGCACCAATATAAAAACACAGTTTGTCCACAGAAGGAGACagcg       >  2:485250/3‑150 (MQ=255)
                                                                                                                                                gCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgcgtg     >  1:139062/1‑150 (MQ=255)
                                                                                                                                                 cTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCCCCAATATCAATCCACAGTTTGTCGACAGAAGGTGACtgcgtgc    >  2:700223/1‑150 (MQ=255)
                                                                                                                                                   ttACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTgcgc  <  2:1233220/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGC‑CTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGC  >  NC_000913/4530976‑4531271

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