Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 533,849 G→A 7.0% intergenic (+23/‑67) allR → / → gcl glyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913533,8490GA 7.0% 613.1 / 26.0 196intergenic (+23/‑67)allR/gclglyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base G (93/89);  new base A (7/7);  total (100/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-02

Only 100 of 198 total aligned reads displayed.

CCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGA  >  NC_000913/533702‑533996
                                                                                                                                                   |                                                                                                                                                   
  gCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAg                                                                                                                                                  >  4:27694/1‑149 (MQ=255)
  gCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGATGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAg                                                                                                                                                  >  3:911750/1‑149 (MQ=255)
     atctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCgg                                                                                                                                               <  2:243091/149‑1 (MQ=255)
      tctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCggg                                                                                                                                              >  1:132277/1‑149 (MQ=255)
       ctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                             <  3:101738/149‑2 (MQ=255)
         ccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttt                                                                                                                                           <  2:137615/149‑1 (MQ=255)
            tctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                        <  2:108309/149‑1 (MQ=255)
             ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                       <  2:1093981/149‑1 (MQ=255)
             ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                       <  4:619984/149‑1 (MQ=255)
             ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                       <  4:1453855/149‑1 (MQ=255)
             ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                       <  1:997742/149‑1 (MQ=255)
                         atcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGa                                                                                                                           >  4:1027763/1‑149 (MQ=255)
                         atcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGa                                                                                                                           >  3:655063/1‑149 (MQ=255)
                          tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCATATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                          >  3:1380999/1‑149 (MQ=255)
                          tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                          >  2:1163896/1‑149 (MQ=255)
                              gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                      <  3:845486/149‑1 (MQ=255)
                              gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                      <  1:203729/149‑1 (MQ=255)
                              gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                      <  1:1261863/149‑1 (MQ=255)
                                 tGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGa                                                                                                                   >  1:910315/1‑149 (MQ=255)
                                      gaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGAGTGAAAGTACATCCATAATGTCTGTAGCATCCCGCTCTTCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAt                                                                                                              >  1:1059368/1‑149 (MQ=255)
                                      gaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAt                                                                                                              <  3:1100516/149‑1 (MQ=255)
                                        agaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAttt                                                                                                            >  2:926918/1‑149 (MQ=255)
                                         gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAATtttt                                                                                                           >  1:829707/1‑149 (MQ=255)
                                         gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                           <  3:813663/149‑1 (MQ=255)
                                          aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTgg                                                                                                                    >  2:731643/1‑139 (MQ=255)
                                          aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                          >  2:185429/1‑149 (MQ=255)
                                           tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTcc                                                                                                         <  1:265909/149‑1 (MQ=255)
                                             gTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaa                                                                                                       >  1:572620/1‑149 (MQ=255)
                                              ttttGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAATTTTTCCaac                                                                                                      <  2:829707/149‑2 (MQ=255)
                                                 tGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaa                                                                                                   <  4:827552/149‑1 (MQ=255)
                                                    cagtcagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATag                                                                                                <  3:458519/149‑1 (MQ=255)
                                                         agGGTGAGCTGGTCGGAGACACCGCCCGCGATATCAGCACGGCGTGGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAAGAGAGGTa                                                                                           >  4:1364315/1‑149 (MQ=255)
                                                                 cTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaa                                                                                   >  1:1088209/1‑149 (MQ=255)
                                                                 cTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaa                                                                                   <  4:281194/149‑1 (MQ=255)
                                                                     tCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaa                                                                                    <  4:162445/144‑1 (MQ=255)
                                                                     tCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATgg                                                                               >  3:984549/1‑149 (MQ=255)
                                                                        gagaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaa                                                                            >  4:440248/1‑149 (MQ=255)
                                                                         agaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaaa                                                                           >  4:471480/1‑149 (MQ=255)
                                                                         agaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaaa                                                                           >  2:1336339/1‑149 (MQ=255)
                                                                         agaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaaa                                                                           >  1:809132/1‑149 (MQ=255)
                                                                          gaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaaaa                                                                          >  4:219012/1‑149 (MQ=255)
                                                                             acCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATg                                                                       >  4:1214666/1‑149 (MQ=255)
                                                                             acCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATg                                                                       >  2:485060/1‑149 (MQ=255)
                                                                              ccgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATga                                                                      <  1:1357903/149‑1 (MQ=255)
                                                                              ccgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGAGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATga                                                                      <  2:1235756/149‑1 (MQ=255)
                                                                              ccgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATACCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATga                                                                      <  2:132277/149‑1 (MQ=255)
                                                                                gcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATgaga                                                                    <  1:251629/149‑1 (MQ=255)
                                                                                gcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATgaga                                                                    <  3:1043352/149‑1 (MQ=255)
                                                                                   cgcgATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGcc                                                                 <  1:539970/149‑1 (MQ=255)
                                                                                      gATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGtt                                                              <  3:471480/149‑1 (MQ=255)
                                                                                      gATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGtt                                                              <  1:603662/149‑1 (MQ=255)
                                                                                          tCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAATATTTTCCAATAAATAGAGGTAg                                                                                          >  2:326424/1‑117 (MQ=255)
                                                                                          tCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAg                                                                                          <  1:326424/117‑1 (MQ=255)
                                                                                          tCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGAcg                                                          >  3:589280/1‑149 (MQ=255)
                                                                                           cAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGAcgc                                                         <  1:1332280/149‑1 (MQ=255)
                                                                                                   cGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgt                                                 <  3:1437991/149‑1 (MQ=255)
                                                                                                   cGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgt                                                 <  2:130544/149‑1 (MQ=255)
                                                                                                   cGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgt                                                 <  3:386265/149‑1 (MQ=255)
                                                                                                    gTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAAt                                                                        <  2:427559/125‑1 (MQ=255)
                                                                                                    gTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAAt                                                                        >  1:427559/1‑125 (MQ=255)
                                                                                                       gggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAACAATGAGAGCCGTTGACGCGGCAatgtatgt                                             >  3:1450075/1‑149 (MQ=255)
                                                                                                         gACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                                  <  1:1127702/110‑1 (MQ=255)
                                                                                                         gACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGc                                                      <  4:1413519/138‑1 (MQ=255)
                                                                                                            tGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTgg                                        <  4:1133743/149‑1 (MQ=255)
                                                                                                              aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGcc                                                                 <  2:59/122‑1 (MQ=255)
                                                                                                              aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGcc                                                                 >  1:59/1‑122 (MQ=255)
                                                                                                              aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGgag                                      >  3:229584/1‑149 (MQ=255)
                                                                                                              aaaGCACATCCATAATGTCTGTCGCAACCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTGTACAGTTGTAAAAATTTTCCAATAAATAGAGGTAGGAAAAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGgag                                      >  4:1060329/1‑149 (MQ=255)
                                                                                                                    caTCCATAATGTCTGTCGCATCCCGCTCTGCGGGGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGagaaagaa                                <  1:695045/149‑1 (MQ=255)
                                                                                                                    caTCCATAATGTCTGGCGCATCCCGCGCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGGTGACGCGGCAAGGTAGGTGCTGGaggaagaa                                <  2:254863/149‑1 (MQ=255)
                                                                                                                     aTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAaagaag                               >  3:19932/1‑149 (MQ=255)
                                                                                                                     aTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAaagaag                               >  2:1406979/1‑149 (MQ=255)
                                                                                                                     aTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAaagaag                               >  1:533196/1‑149 (MQ=255)
                                                                                                                      tCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                           >  4:260060/1‑72 (MQ=255)
                                                                                                                      tCCATAATGCCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                           <  3:260060/72‑1 (MQ=255)
                                                                                                                        cATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATagag                                                                                              <  3:1275607/83‑1 (MQ=255)
                                                                                                                         aTAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTAt                           <  3:881774/149‑1 (MQ=255)
                                                                                                                          tAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATc                          >  2:809210/1‑149 (MQ=255)
                                                                                                                           aaTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGaa                                                                                       <  1:814195/87‑1 (MQ=255)
                                                                                                                           aaTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGaa                                                                                       >  2:814195/1‑87 (MQ=255)
                                                                                                                           aaTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCa                         <  1:1442971/149‑1 (MQ=255)
                                                                                                                            aTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCAc                        >  1:948864/1‑149 (MQ=255)
                                                                                                                            aTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCAc                        <  4:1137717/149‑1 (MQ=255)
                                                                                                                               tctgtcGCATCCCGCTCTGCGGAGCGTGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTAc                     <  3:95150/149‑1 (MQ=255)
                                                                                                                                 tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCg                   >  2:232367/1‑149 (MQ=255)
                                                                                                                                    aaaaacccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGagaaagaa                                >  4:886446/6‑133 (MQ=255)
                                                                                                                                    aaaaacccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCt                >  4:701273/6‑149 (MQ=255)
                                                                                                                                      aaacccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTc              >  3:758017/4‑149 (MQ=255)
                                                                                                                                       aTCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCg             <  4:1312832/149‑1 (MQ=255)
                                                                                                                                       aTCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCg             <  2:856907/149‑1 (MQ=255)
                                                                                                                                        tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCgg            >  2:639927/1‑149 (MQ=255)
                                                                                                                                        tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCgg            >  4:3601/1‑149 (MQ=255)
                                                                                                                                        tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAG‑‑AGAAGGTATCACTACCGCCTTCGgtg          <  2:561053/149‑1 (MQ=255)
                                                                                                                                         cccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGACTTCGgt           <  4:758017/149‑1 (MQ=255)
                                                                                                                                          ccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGgtg          >  3:719281/1‑149 (MQ=255)
                                                                                                                                          ccGATCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGgtg          <  4:984549/149‑1 (MQ=255)
                                                                                                                                             ctctGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTc       >  4:532299/1‑149 (MQ=255)
                                                                                                                                                cgCGGAGCGGGTTTTTTTGTCAAAATTTGAAAGTTGGAAAAATTTTCCAGTAAATAGAGGCAGGCATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGCAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCgg    <  2:1059368/148‑1 (MQ=255)
                                                                                                                                                 gCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCggg   >  1:1035550/1‑149 (MQ=255)
                                                                                                                                                  cGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGa  <  2:539119/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGA  >  NC_000913/533702‑533996

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