Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,749,465 T→A 5.0% R78R (CGT→CGA sgbU → putative L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,4650TA 5.0% 2346.0 / 71.3 635R78R (CGT→CGAsgbUputative L‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (305/298);  new base A (17/15);  total (322/313)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 692 total aligned reads displayed.

GCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCA  >  NC_000913/3749317‑3749613
                                                                                                                                                    |                                                                                                                                                    
       cataacctgccagctggatggtgcgaatgccgagatcgcgcgccaggcgaatcgctttactcatgatttcccgcgcccgttcgcgcactgcctcgtcacggctaccgaagggatcccgtCGATGTGCCTGTCCGCACATCGTCGCTTTc                                                                                                                                               >  1:631769‑M1/120‑149 (MQ=255)
        gATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCtt                                                                                                                                                 <  4:19520/146‑1 (MQ=255)
          ttttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCt                                                                                                                                            >  1:927298/1‑149 (MQ=255)
          ttttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCt                                                                                                                                            >  2:1364849/1‑149 (MQ=255)
          ttttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCt                                                                                                                                            >  1:316810/1‑149 (MQ=255)
          ttttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCt                                                                                                                                            >  3:844677/1‑149 (MQ=255)
            ttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCttt                                                                                                                                          <  2:750805/149‑1 (MQ=255)
            ttGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCttt                                                                                                                                          <  3:320908/149‑1 (MQ=255)
                        cGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGAc                                                                                                                              <  1:829376/149‑1 (MQ=255)
                           tGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg                                                                                                                           <  3:462221/149‑1 (MQ=255)
                              aTGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCa                                                                                                                        <  4:313463/149‑1 (MQ=255)
                               tGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTgg                                                                                                                                        <  4:235564/132‑1 (MQ=255)
                                gAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGt                                                                                                                      >  4:1273665/1‑149 (MQ=255)
                                gAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGt                                                                                                                      <  3:1013053/149‑1 (MQ=255)
                                  aaCCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTcc                                                                                                                                              >  1:159414/1‑123 (MQ=255)
                                  aaCCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAgg                                                                                                                          >  3:359813/1‑143 (MQ=255)
                                  aaCCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAgg                                                                                                                          <  4:359813/143‑1 (MQ=255)
                                    ccGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTa                                                                                                                                      >  1:1427206/1‑129 (MQ=255)
                                          aaCGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACgg                                                                                                            <  3:1227173/149‑1 (MQ=255)
                                                tctcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACggg                                                                                                           <  3:278974/144‑1 (MQ=255)
                                                tctcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACggg                                                                                                           >  4:278974/1‑144 (MQ=255)
                                                tctcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCgg                                                                                                      >  1:1291939/1‑149 (MQ=255)
                                                 ctcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGa                                                                                                                            <  4:975654/126‑1 (MQ=255)
                                                  tcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTccc                                                                                                                                             <  4:923771/108‑1 (MQ=255)
                                                  tcACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGcgcg                                                                                                       <  4:538140/146‑1 (MQ=255)
                                                       tCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCa                                                                                               >  1:368892/1‑149 (MQ=255)
                                                          tGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGa                                                                                            <  3:739744/149‑1 (MQ=255)
                                                          tGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGa                                                                                            <  4:1374783/149‑1 (MQ=255)
                                                           gACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGc                                                                                                                         >  4:74283/1‑119 (MQ=255)
                                                           gACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAg                                                                                           <  3:542422/149‑1 (MQ=255)
                                                                gAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCTGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAg                                                                                      <  3:974111/149‑1 (MQ=255)
                                                                gAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAg                                                                                      >  1:786703/1‑149 (MQ=255)
                                                                   tcatgatttcccgcgcccgttcgcgcactgcctcgtcacggctaccaaagggaaagcgtCGATGTGCCTGTCCGCACAtcgtc                                                                                                                                                     >  3:418991‑M1/60‑83 (MQ=255)
                                                                    gccgcACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGACGCTTTCCCTTTGGTAGACGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGAt                                                                                  <  4:561162/149‑1 (MQ=255)
                                                                      tgatttcccgcgcccgttcgcgcactgcctcgtcacggctaccaaagggaaagcgtCGATGTGCCTGTCCGCACATCGTCGc                                                                                                                                                   >  3:1400466‑M1/57‑82 (MQ=255)
                                                                        cacaAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGa                                                                                   <  4:607986/144‑1 (MQ=255)
                                                                        cacaAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGa                                                                                   >  3:607986/1‑144 (MQ=255)
                                                                         acaAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGcc                                                                             <  2:316810/149‑1 (MQ=255)
                                                                         acaAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGcc                                                                             <  2:368892/149‑1 (MQ=255)
                                                                           aaaGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGa                                                                                                                            <  3:1185804/100‑1 (MQ=255)
                                                                           aaaGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTg                                                                           >  3:564952/1‑149 (MQ=255)
                                                                           aaaGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTg                                                                           <  2:786703/149‑1 (MQ=255)
                                                                            aaGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTgg                                                                          <  4:417643/149‑1 (MQ=255)
                                                                               gACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGaa                                                                                                               <  3:658730/109‑1 (MQ=255)
                                                                               gACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGaa                                                                                                               >  4:658730/1‑109 (MQ=255)
                                                                                aCTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcg                                                                      <  2:510618/149‑1 (MQ=255)
                                                                                  ttCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTa                                                                                                                                      <  4:112883/83‑1 (MQ=255)
                                                                                  ttCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCgg                                                                                                      <  3:190111/115‑1 (MQ=255)
                                                                                  ttCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCAt                                                                                              >  3:225782/1‑123 (MQ=255)
                                                                                  ttCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTg                                                                           >  1:1005866/1‑142 (MQ=255)
                                                                                        ggTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCg                                                              <  1:583273/149‑1 (MQ=255)
                                                                                          ttGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGa                                                                                            <  4:1457628/117‑1 (MQ=255)
                                                                                              cgcgATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTc                                                        >  1:1111790/1‑149 (MQ=255)
                                                                                                 gatgatCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCa                                                     <  2:89590/149‑1 (MQ=255)
                                                                                                  atgatCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTg                                                                                                                                         <  4:1075763/64‑1 (MQ=255)
                                                                                                  atgatCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCAc                                                    <  3:1360322/149‑1 (MQ=255)
                                                                                                     atCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGt                                                                                                                                       <  1:1062277/63‑1 (MQ=255)
                                                                                                     atCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCAt                                                 >  2:412371/1‑149 (MQ=255)
                                                                                                      tCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCAc                                                    <  2:1404466/145‑1 (MQ=255)
                                                                                                      tCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATc                                                >  2:472031/1‑149 (MQ=255)
                                                                                                         aaaCAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAAc                                                                                                              >  1:1374066/1‑84 (MQ=255)
                                                                                                         aaaCAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAg                                             >  3:1284287/1‑149 (MQ=255)
                                                                                                          aaCAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGaaa                                                                                                  <  1:1243310/95‑1 (MQ=255)
                                                                                                            cAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTg                                          <  1:412371/149‑1 (MQ=255)
                                                                                                               gtgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCa                                       <  2:1205136/149‑1 (MQ=255)
                                                                                                                tgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCa                                       <  1:230139/148‑1 (MQ=255)
                                                                                                                tgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAg                                      <  2:1283638/149‑1 (MQ=255)
                                                                                                                tgtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAg                                      >  1:1325543/1‑149 (MQ=255)
                                                                                                                 gtTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGAt                                                                  >  4:1035347/1‑120 (MQ=255)
                                                                                                                  ttGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCg                                                              <  2:215141/123‑1 (MQ=255)
                                                                                                                    gggATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAg                                                                                           <  2:949260/92‑1 (MQ=255)
                                                                                                                     ggATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTgg                                         <  2:848702/141‑1 (MQ=255)
                                                                                                                      gATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATg                                >  2:100849/1‑149 (MQ=255)
                                                                                                                       aTCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcgc                                                                     >  4:859151/1‑111 (MQ=255)
                                                                                                                        tCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGt                                    >  1:727068/1‑143 (MQ=255)
                                                                                                                          ccGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGt                            >  2:148196/1‑149 (MQ=255)
                                                                                                                           cGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTc                           >  1:959785/1‑149 (MQ=255)
                                                                                                                            gacgatgtgcggacaggcacatcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCa                                                                                               <  4:1351368‑M2/56‑1 (MQ=255)
                                                                                                                            gacgatgtgcggacaggcacatcgACGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGa                                                                                   >  1:889292‑M2/25‑92 (MQ=255)
                                                                                                                            gacgatgtgcggacaggcacatcgACGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCt                          >  3:1454531‑M2/25‑149 (MQ=255)
                                                                                                                            gTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAg                                                                                      <  1:474249/89‑1 (MQ=255)
                                                                                                                            gTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCAt                                                 <  3:791934/126‑1 (MQ=255)
                                                                                                                            gTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGt                            <  2:1311108/147‑1 (MQ=255)
                                                                                                                               GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTa                                                                                                                                      >  1:513294‑M1/1‑38 (MQ=255)
                                                                                                                                 tgtgcggacaggcacatcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTAc                     >  3:1092949‑M2/20‑149 (MQ=255)
                                                                                                                                   TGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTAt                                 <  4:1197808‑M1/135‑1 (MQ=255)
                                                                                                                                   TGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCACACCATCCAGCTGGCAGGTTATGACGTCTATTAc                     >  2:645695‑M1/1‑147 (MQ=255)
                                                                                                                                     CCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGa                                                                                   >  2:256717‑M1/1‑83 (MQ=255)
                                                                                                                                     CCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTg                                                                           >  4:79826‑M1/1‑91 (MQ=255)
                                                                                                                                        gTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGAGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTAt                        >  3:1062357/1‑139 (MQ=255)
                                                                                                                                            gcacatcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCacga          >  4:449076‑M2/9‑149 (MQ=255)
                                                                                                                                              acatcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTa                      <  2:1126464‑M2/129‑1 (MQ=255)
                                                                                                                                              acatcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCacgacg        >  2:440301‑M2/7‑149 (MQ=255)
                                                                                                                                                ctcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATg                                >  1:1146904‑M2/5‑123 (MQ=255)
                                                                                                                                                atcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTc                           >  2:1075888‑M2/5‑128 (MQ=255)
                                                                                                                                                 tcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAg     <  3:890125‑M2/146‑1 (MQ=255)
                                                                                                                                                 tcgTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAg     <  2:558940‑M2/146‑1 (MQ=255)
                                                                                                                                                   gTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCAt                                                                                              >  1:446139‑M2/2‑58 (MQ=255)
                                                                                                                                                    tcgCTTTCCCTTTGGTAGCTGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGaa                  <  3:999523/133‑1 (MQ=255)
                                                                                                                                                    tcgCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCa  >  1:258545/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCA  >  NC_000913/3749317‑3749613

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: