Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 533,849 G→A 5.4% 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 5.4% 305.8 / 5.4 93intergenic (+23/‑67)allR/gclglyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base A (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

ATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAG  >  NC_000913/533707‑533997
                                                                                                                                              |                                                                                                                                                    
atctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCgg                                                                                                                                                <  1:346273/149‑1 (MQ=255)
atctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCgg                                                                                                                                                <  3:260372/149‑1 (MQ=255)
   tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                             >  1:1087/1‑149 (MQ=255)
   tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                             >  2:492040/1‑149 (MQ=255)
   tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                             >  3:543023/1‑149 (MQ=255)
   tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                             >  2:306764/1‑149 (MQ=255)
   tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                             >  3:521862/1‑149 (MQ=255)
    ccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttt                                                                                                                                            >  4:147662/1‑149 (MQ=255)
    ccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttt                                                                                                                                            >  4:331406/1‑149 (MQ=255)
      atctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                          >  4:462555/1‑149 (MQ=255)
      atctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                          >  1:136514/1‑149 (MQ=255)
      atctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                          >  1:151291/1‑149 (MQ=255)
       tctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                         <  1:468193/149‑1 (MQ=255)
       tctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                         >  3:258042/1‑149 (MQ=255)
        ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                        <  2:1087/149‑1 (MQ=255)
          ccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                          <  3:50878/145‑1 (MQ=255)
          ccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                          >  4:50878/1‑145 (MQ=255)
           cGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGAc                                                                                                                                     <  1:301664/149‑1 (MQ=255)
               ccTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACaaaa                                                                                                                                 >  2:477636/1‑149 (MQ=255)
                cTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAt                                                                                                                                >  1:461163/1‑149 (MQ=255)
                cTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAt                                                                                                                                >  2:12565/1‑149 (MQ=255)
                 ttCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAtt                                                                                                                               >  1:284743/1‑149 (MQ=255)
                 ttCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAtt                                                                                                                               <  2:551475/149‑1 (MQ=255)
                  tcatcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAttt                                                                                                                              <  2:170439/149‑1 (MQ=255)
                    atcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGa                                                                                                                            >  2:489829/1‑149 (MQ=255)
                     tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                           >  1:597386/1‑149 (MQ=255)
                         gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                       <  2:42/149‑1 (MQ=255)
                         gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                       <  4:194616/149‑1 (MQ=255)
                         gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                       <  2:429965/149‑1 (MQ=255)
                         gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                       <  4:258042/149‑1 (MQ=255)
                         gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                       <  3:379202/149‑1 (MQ=255)
                                agaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATATCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGaaaaa                                                                                                                >  2:586697/1‑149 (MQ=255)
                                 gaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAt                                                                                                               <  1:492040/149‑1 (MQ=255)
                                  aagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtt                                                                                                              >  1:19350/1‑149 (MQ=255)
                                    gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                            <  3:175684/149‑1 (MQ=255)
                                    gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                            <  1:65895/149‑1 (MQ=255)
                                    gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                            >  4:390920/1‑149 (MQ=255)
                                     aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                           <  3:147662/149‑1 (MQ=255)
                                     aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                           <  3:488348/149‑1 (MQ=255)
                                     aTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                           >  4:434110/1‑149 (MQ=255)
                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                        >  1:210050/1‑119 (MQ=255)
                                      tCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                        <  2:210050/119‑1 (MQ=255)
                                        gTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaa                                                                                                        <  4:16075/149‑1 (MQ=255)
                                        gTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaa                                                                                                        >  3:259784/1‑149 (MQ=255)
                                          tttGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaata                                                                                                      >  2:85534/1‑149 (MQ=255)
                                             gtcagtcaGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaat                                                                                                   >  3:339672/1‑149 (MQ=255)
                                             gtcagtcaGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaat                                                                                                   >  2:147174/1‑149 (MQ=255)
                                             gtcagtcaGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaat                                                                                                   <  3:171157/149‑1 (MQ=255)
                                                  tcagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGTACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAgg                                                                                              >  3:562482/1‑149 (MQ=255)
                                                  tcagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAgg                                                                                              >  3:13446/1‑149 (MQ=255)
                                                  tcagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAgg                                                                                              >  3:209101/1‑149 (MQ=255)
                                                     gggTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAg                                                                                           <  2:597386/149‑1 (MQ=255)
                                                        tGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGaa                                                                                        <  2:202241/149‑1 (MQ=255)
                                                             tGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                                   <  3:20282/149‑1 (MQ=255)
                                                              ggTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGATTAAAAt                                                                                  <  2:422712/149‑1 (MQ=255)
                                                                tCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATgg                                                                                >  1:335540/1‑149 (MQ=255)
                                                                      acacCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCaaaaa                                                                          >  2:407107/1‑149 (MQ=255)
                                                                         ccgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATga                                                                       <  4:454890/149‑1 (MQ=255)
                                                                              cgcgATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGcc                                                                  <  2:397668/149‑1 (MQ=255)
                                                                                       aGCACGGCGTTGGGACTGAAAGCACTTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGACAGTTGGAAAAAttt                                                                                                             <  2:286066/97‑1 (MQ=255)
                                                                                       aGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAttt                                                                                                             >  1:286066/1‑97 (MQ=255)
                                                                                        gCACGGCGTTGGGACTGAAAGCATATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCgg                                                        >  3:420859/1‑149 (MQ=255)
                                                                                        gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCgg                                                        <  2:19350/149‑1 (MQ=255)
                                                                                           cGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCaa                                                     <  2:151291/149‑1 (MQ=255)
                                                                                                ttGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                                <  1:147174/149‑1 (MQ=255)
                                                                                                 tGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtatg                                               <  1:12565/149‑1 (MQ=255)
                                                                                                  gggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATGTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtatgt                                              >  4:349469/1‑149 (MQ=255)
                                                                                                    gACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAACGAGAGCCGTTGACGCGGCAATGTATGTGc                                            <  3:349469/149‑1 (MQ=255)
                                                                                                      cTGAAAGCACATCTATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTg                                          <  3:434110/149‑1 (MQ=255)
                                                                                                      cTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTg                                          <  3:507350/149‑1 (MQ=255)
                                                                                                        gAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGga                                        <  3:571615/149‑1 (MQ=255)
                                                                                                         aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGgag                                       >  1:161066/1‑149 (MQ=255)
                                                                                                         aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGgag                                       >  1:320710/1‑149 (MQ=255)
                                                                                                                 tCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAgg                               >  2:279923/1‑149 (MQ=255)
                                                                                                                    ataataTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTAt                            <  1:586697/149‑1 (MQ=255)
                                                                                                                        tgtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                                <  3:537729/125‑1 (MQ=255)
                                                                                                                        tgtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                                >  4:537729/1‑125 (MQ=255)
                                                                                                                         gtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTa                       >  4:50646/1‑149 (MQ=255)
                                                                                                                            tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTAt                            >  3:150964/1‑141 (MQ=255)
                                                                                                                            tgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTAt                            <  4:150964/141‑1 (MQ=255)
                                                                                                                              tcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGcc                  >  2:298378/1‑149 (MQ=255)
                                                                                                                              aaaaaacccGCTCCGCAGAGCGGGTTTCTTTGACAAAATTTGAAAGTTGTAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGcc                  >  2:177075/7‑149 (MQ=255)
                                                                                                                                  aacccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCg              >  1:210133/3‑149 (MQ=255)
                                                                                                                                  aTCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCg              <  1:279923/149‑1 (MQ=255)
                                                                                                                                   tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTACGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCgg             <  4:420859/149‑1 (MQ=255)
                                                                                                                                     ccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTc               <  4:165737/145‑1 (MQ=255)
                                                                                                                                     ccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTc               >  3:165737/1‑145 (MQ=255)
                                                                                                                                         tctGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTcc       <  3:453573/149‑1 (MQ=255)
                                                                                                                                          ctGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCg      <  1:321446/149‑1 (MQ=255)
                                                                                                                                           tGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCgg     <  4:209101/149‑1 (MQ=255)
                                                                                                                                             cGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGa   <  4:543023/149‑1 (MQ=255)
                                                                                                                                             cGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGa   <  3:202558/149‑1 (MQ=255)
                                                                                                                                             cGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGa   <  3:542034/149‑1 (MQ=255)
                                                                                                                                              ggAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAACAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAg  <  3:495824/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTTCCGGGAG  >  NC_000913/533707‑533997

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: