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.

CACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGC  >  NC_000913/4530976‑4531271
                                                                                                                                                   |                                                                                                                                                     
 acGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACATCCTAGACGATCGTCTAGTGCTTTACt                                                                                                                                                   >  2:898065/1‑150 (MQ=255)
   ggCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGc                                                                                                                                                 <  1:1317050/150‑1 (MQ=39)
    gcttAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGcc                                                                                                                                                <  1:139041/147‑2 (MQ=38)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcac                                                                                                                                               >  1:819732/1‑150 (MQ=255)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGcac                                                                                                                                               <  1:465569/150‑1 (MQ=255)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                                                                               <  2:1425080/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                                                                               <  2:1311781/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                                                                               <  2:890641/150‑3 (MQ=39)
     cGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccc                                                                                                                                               <  1:1711898/150‑3 (MQ=39)
       tAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                                                                             >  2:1488342/1‑150 (MQ=255)
       tAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATAGCACGCTCCATAGTACTGATACCATCCGCTACACCGCGCATGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCacac                                                                                                                                             <  1:951731/150‑1 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            <  2:1276787/150‑1 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            >  1:1251712/1‑150 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACAcc                                                                                                                                            >  1:1005827/1‑150 (MQ=255)
        aaTACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGccctgc                                                                                                                                            <  1:111596/150‑6 (MQ=38)
          tACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAg                                                                                                                                          >  1:1091550/1‑150 (MQ=255)
            cGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCTTCTAGTGCTGTACTGCACACCAGag                                                                                                                                        >  1:664427/1‑149 (MQ=255)
             gggAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATc                                                                                                                                       <  2:1413652/150‑1 (MQ=255)
                 aCAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                                                                   <  2:819732/150‑1 (MQ=255)
                 aCAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGttt                                                                                                                                   <  1:1343670/150‑1 (MQ=255)
                    aTGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGcc                                                                                                                                >  1:405543/1‑150 (MQ=255)
                      gCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAg                                                                                                                              <  1:693119/150‑1 (MQ=255)
                         gcgGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAg                                                                                                                           <  2:891575/150‑1 (MQ=255)
                           ggATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGtt                                                                                                                         >  2:670932/1‑150 (MQ=255)
                               atTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcgg                                                                                                                     >  1:562460/1‑150 (MQ=255)
                                ttAACTTCTTCCTGTACCGACGCCACCAGCTAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggc                                                                                                                    <  1:1246994/150‑1 (MQ=255)
                                 tAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggcg                                                                                                                   <  1:1194263/150‑1 (MQ=255)
                                 tAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTgcggcg                                                                                                                   >  2:313953/1‑150 (MQ=255)
                                  ataTTCTTCCTGTACCGAAGCCACCAGATAAACGGCTATATCAAGCTCCATAGTACTGTTAGCATCCGCTACACCGAGCAGGGCAGTACAGGCTAGAGGATCGTGTAGTGCTTTACTGCAGACAAGATCGTTTGCCAGGAGTTGcggcgg                                                                                                                  <  1:1229372/148‑1 (MQ=255)
                                   aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                                 <  1:1604621/150‑1 (MQ=255)
                                   aaTTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCggg                                                                                                                 <  2:583985/150‑1 (MQ=255)
                                     ttcttcCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGtt                                                                                                               <  2:230809/150‑1 (MQ=255)
                                      tcttcCTGTACCGAAGCCACCAGATAAACGGAGATATCAAGCTCCATAGTACTTAAAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTta                                                                                                              <  2:1065657/150‑1 (MQ=255)
                                           cTGTACCGAAGCCACCAGATAAACGGCGATATCACGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATActg                                                                                                         >  2:323785/1‑148 (MQ=255)
                                            tGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACaga                                                                                                        >  2:1588974/1‑150 (MQ=255)
                                             gTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACCGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATACagag                                                                                                       >  2:137042/1‑150 (MQ=255)
                                              tACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGt                                                                                                      <  2:619456/150‑1 (MQ=255)
                                               aCCGAATCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACTCCGAGCATAGCAGTACAGCCTAGACTATCGTCTAGTGCTTTACTGCACACCATATCGTTTGCAAGGAGATGCGTCGGGTTATACagagag                                                                                                     >  2:21311/1‑148 (MQ=255)
                                                   aaGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGGTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAg                                                                                                 <  1:1609608/150‑1 (MQ=255)
                                                    aGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAgg                                                                                                <  1:662888/150‑1 (MQ=255)
                                                      ccaccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCg                                                                                              <  2:1562445/150‑1 (MQ=255)
                                                       caccaGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGt                                                                                             >  2:200768/1‑150 (MQ=255)
                                                           aGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACATAGTGACAGGCGTACCg                                                                                         >  2:1040785/1‑150 (MQ=255)
                                                               aaaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                     <  1:401257/150‑1 (MQ=255)
                                                               aaaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCt                                                                                     >  2:928810/1‑150 (MQ=255)
                                                                aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCATGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  1:1311263/150‑1 (MQ=255)
                                                                aaCGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  1:1308889/150‑1 (MQ=255)
                                                                aaCGGCGAAATCCAGCTCCATACTACTCATAGCAGCCGCTACACCGAGCAGGGCAGTACTGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGGTATACAGAGTGACAGGCGTACCGACCTg                                                                                    <  2:1073980/150‑1 (MQ=255)
                                                                  cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                                                                  >  1:1154169/1‑150 (MQ=255)
                                                                  cGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGaa                                                                                  >  1:1420365/1‑150 (MQ=255)
                                                                   ggCGATATCAAGCTCCATAGTACTGATAGCATCCGCTAAAACGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCATATCGTTTGCCAGGAGTTGATGCGGGTTATACAGAGTGACAGGCGTACAGACCTGAAt                                                                                 >  2:720771/1‑150 (MQ=255)
                                                                    gCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATg                                                                                <  1:922216/150‑1 (MQ=255)
                                                                      gATATCAAGCTCCATAGTACTGATAGCATCCGCTACCCCGAGCAGGGCAGTACAGCCTAGACGCTCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGCCCGGAATGcc                                                                              >  2:1020941/1‑150 (MQ=255)
                                                                      gATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGcc                                                                              >  1:885111/1‑150 (MQ=255)
                                                                        tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGTCAGGCGTACCGACCTGAATGCCCa                                                                            >  2:887090/1‑150 (MQ=255)
                                                                        tatCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGACTGCCCa                                                                            >  1:1182471/1‑150 (MQ=255)
                                                                           cAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTGTGCCAGGAGTTGCGGCGGGTTATACCGAGTGACAGGCGTACCGCCCTGAATGCCCATcc                                                                         >  2:1190051/1‑150 (MQ=255)
                                                                            aaGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCg                                                                        <  1:1352767/150‑1 (MQ=255)
                                                                            aaGCTCCAAAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCCTCCg                                                                        >  1:1128272/1‑150 (MQ=255)
                                                                              gCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGta                                                                      >  2:1078956/1‑150 (MQ=255)
                                                                                   aTAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcg                                                                 <  1:1488342/150‑1 (MQ=255)
                                                                                     aGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGcgtc                                                               <  1:200768/150‑1 (MQ=255)
                                                                                          tGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctt                                                          <  1:1443736/150‑1 (MQ=255)
                                                                                           gATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGATTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGGATAGCGTCgtcttt                                                         >  1:1403237/1‑150 (MQ=255)
                                                                                            aTAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtcgttg                                                        >  1:1623601/1‑150 (MQ=255)
                                                                                              aGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtc                                                      <  2:927976/150‑1 (MQ=255)
                                                                                                cATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtctt                                                    <  1:1286510/150‑1 (MQ=255)
                                                                                                cATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCgtctttgtctt                                                    <  1:1726086/150‑1 (MQ=255)
                                                                                                  tCCGCTAAACCGAGCAGGGCAGTACAGCATAGACGATCGTCTAGTGCTTTACTGCACACCAGAGCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCACATCCGTATAGCGTCgtcgttgtagttt                                                  >  1:603567/1‑145 (MQ=255)
                                                                                                   ccGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTgg                                                 >  1:1603822/1‑150 (MQ=255)
                                                                                                    cGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTGTGCCTGGAGTTGCGGCGGTTTATACAGAGTGACAGGCGTACCGACCTGAAGGCCCATCCGTATAGCGTCGTCtttggctttggc                                                >  2:863938/1‑150 (MQ=255)
                                                                                                         cacCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATCCCCCTCCGTATAGCGTCGTGTTTGTCTTTGGCACCaa                                           >  2:832650/1‑150 (MQ=255)
                                                                                                           ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAata                                         >  2:39475/1‑150 (MQ=255)
                                                                                                           ccGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAata                                         <  1:1078956/150‑1 (MQ=255)
                                                                                                                cAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAAt                                    >  2:1304604/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATcc                                  >  2:1721334/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATcc                                  >  2:841458/1‑150 (MQ=255)
                                                                                                                  gggCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCAACAATATCAAAcc                                  >  1:1417257/1‑150 (MQ=255)
                                                                                                                     cAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcaca                               <  1:175823/150‑1 (MQ=255)
                                                                                                                     cAGTAAAGCACTAGACGATCGTCTAGGCTGTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCcaca                               >  1:1251166/1‑150 (MQ=39)
                                                                                                                      aGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAg                              <  2:1498321/150‑1 (MQ=255)
                                                                                                                       gTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAAAATCAATCAACAGt                             >  2:1583073/1‑150 (MQ=255)
                                                                                                                            gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                        >  2:1147990/1‑150 (MQ=255)
                                                                                                                            gCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGAGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTc                        <  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)
                                                                                                                                                   |                                                                                                                                                     
CACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAACGGCGATATCAAGCTCCATAGTACTGATAGCATCCGCTACACCGAGCAGGGCAGTACAGCCTAGACGATCGTCTAGTGCTTTACTGCACACCAGATCGTTTGCCAGGAGTTGCGGCGGGTTATACAGAGTGACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGC  >  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: