Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 533,846 T→C 11.0% intergenic (+20/‑70) 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,8460TC11.0% 305.6 / 21.9 100intergenic (+20/‑70)allR/gclglyoxylate‑inducible transcriptional repressor of all and gcl operons/glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base T (42/47);  new base C (4/7);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTT  >  NC_000913/533699‑533990
                                                                                                                                                   |                                                                                                                                                
ttgCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTg                                                                                                                                                 >  1:1443349/1‑149 (MQ=255)
ttgCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTg                                                                                                                                                 <  4:244879/149‑1 (MQ=255)
 tgCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGc                                                                                                                                                <  4:546979/149‑1 (MQ=255)
      cTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGc                                                                                                                                           <  3:754323/149‑1 (MQ=255)
        atctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCgg                                                                                                                                         <  2:1399596/149‑1 (MQ=255)
         tctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggg                                                                                                                                        <  4:767805/149‑1 (MQ=255)
          ctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       <  2:1056496/149‑2 (MQ=255)
          ctccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       <  1:286224/149‑2 (MQ=255)
           tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                      >  1:741220/1‑149 (MQ=255)
           tccatctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtt                                                                                                                                      >  4:281187/1‑149 (MQ=255)
              atctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       <  4:367911/145‑2 (MQ=255)
              atctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       >  3:367911/1‑144 (MQ=255)
               tctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                  >  3:22171/1‑149 (MQ=255)
                ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                 <  3:281187/149‑1 (MQ=255)
                ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                 <  4:22171/149‑1 (MQ=255)
                ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                 <  2:1443349/149‑1 (MQ=255)
                ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                 <  2:741220/149‑1 (MQ=255)
                ctccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttttt                                                                                                                                 <  2:1423542/149‑1 (MQ=255)
                 tccGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTg                                                                                                                                <  2:1195599/149‑1 (MQ=255)
                      gCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       <  2:577597/137‑2 (MQ=255)
                      gCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCggga                                                                                                                                       >  1:577597/1‑136 (MQ=255)
                        cTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCAGATCCATAATGTCTGTCGCATCCCGCTCTGCGGa                                                                                                                                             >  3:1011785/1‑129 (MQ=255)
                        cTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCAGATCCATAATGTCTGTCGCATCCCGCTCTGCGGa                                                                                                                                             <  4:1011785/129‑1 (MQ=255)
                             tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCggg                                                                                                                                        <  4:973172/129‑1 (MQ=255)
                             tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaa                                                                                                                    >  4:1441726/1‑149 (MQ=255)
                             tcaAGACTGACAGAAGATCGTTTTGTCAGTCAGGATGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCggg                                                                                                                                        >  3:973172/1‑129 (MQ=255)
                              caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                  >  1:1166980/1‑134 (MQ=255)
                              caAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttttt                                                                                                                                  <  2:1166980/134‑1 (MQ=255)
                                 gACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                <  3:24458/149‑1 (MQ=255)
                                  aCTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTg                                                                                                               <  3:598323/149‑1 (MQ=255)
                                    tGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                   <  4:1014808/127‑1 (MQ=255)
                                    tGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGttttt                                                                                                                                   >  3:1014808/1‑127 (MQ=255)
                                        agaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttt                                                                                                                                    <  1:152813/122‑1 (MQ=255)
                                        agaagaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGtttt                                                                                                                                    >  2:152813/1‑122 (MQ=255)
                                            gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATACTGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTg                                                                                                               >  1:421805/1‑139 (MQ=255)
                                            gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTg                                                                                                               <  2:421805/139‑1 (MQ=255)
                                            gaTCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAAtttt                                                                                                     >  3:60736/1‑149 (MQ=255)
                                                    tGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGCCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaa                                                                                             <  3:1441726/149‑1 (MQ=255)
                                                           cagGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGt                                                                                      <  2:700430/149‑1 (MQ=255)
                                                            agGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTa                                                                                     >  4:1130995/1‑149 (MQ=255)
                                                             gggTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATagag                                                                                        >  4:1176505/1‑145 (MQ=255)
                                                             gggTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATagag                                                                                        <  3:1176505/145‑1 (MQ=255)
                                                             gggTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAg                                                                                    <  3:1130995/149‑1 (MQ=255)
                                                               gTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCgggat                                                                                                                                      >  4:331744/1‑95 (MQ=255)
                                                               gTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCCGCAGAGCgggat                                                                                                                                      <  3:331744/97‑3 (MQ=255)
                                                                  aGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                >  1:1378337/1‑116 (MQ=255)
                                                                  aGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                <  2:1378337/116‑1 (MQ=255)
                                                                     tGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            <  4:60736/149‑1 (MQ=255)
                                                                          agagACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAtt                                                                                                                        >  1:856339/1‑100 (MQ=255)
                                                                          agagACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAAtt                                                                                                                        <  2:856339/100‑1 (MQ=255)
                                                                          agagACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAGTTTTTCCaac                                                                                                >  3:998964/1‑123 (MQ=255)
                                                                          agagACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGACGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAATTTTTCCaac                                                                                                <  4:998964/124‑2 (MQ=255)
                                                                            agaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATaga                                                                                         <  2:698970/129‑1 (MQ=255)
                                                                            agaCACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATaga                                                                                         >  1:698970/1‑129 (MQ=255)
                                                                              acacCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTgg                                                                                                              <  4:604754/106‑1 (MQ=255)
                                                                              acacCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTgg                                                                                                              >  3:604754/1‑106 (MQ=255)
                                                                                 ccgcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATga                                                                >  4:580829/1‑149 (MQ=255)
                                                                                   gcccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATgaga                                                              <  3:580829/149‑1 (MQ=255)
                                                                                     ccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAg                                                                                    >  2:158416/1‑125 (MQ=255)
                                                                                     ccgcGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAg                                                                                    <  1:158416/125‑1 (MQ=255)
                                                                                             tCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTa                                                                                     <  2:576995/116‑1 (MQ=255)
                                                                                             tCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTa                                                                                     >  1:576995/1‑116 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaaa                                                                                                                   <  2:1204797/83‑1 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGaaa                                                                                                                   >  1:1204797/1‑83 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAAt                                                                           <  3:558107/123‑1 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAAt                                                                           >  4:558107/1‑123 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATgaga                                                              >  4:466862/1‑136 (MQ=255)
                                                                                                gCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATgaga                                                              <  3:466862/136‑1 (MQ=255)
                                                                                                  aCGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            >  2:299739/1‑120 (MQ=255)
                                                                                                  aCGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            <  1:299739/120‑1 (MQ=255)
                                                                                                           ggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATg                                                                          >  3:684326/1‑113 (MQ=255)
                                                                                                           ggACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATg                                                                          <  4:684326/113‑1 (MQ=255)
                                                                                                             aCTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                <  1:932272/73‑1 (MQ=255)
                                                                                                             aCTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGtt                                                                                                                >  2:932272/1‑73 (MQ=255)
                                                                                                             aCTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCt                                    >  4:1410470/1‑149 (MQ=255)
                                                                                                                 aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTGCAATAAATag                                                                                          >  4:870960/1‑91 (MQ=255)
                                                                                                                 aaaGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATag                                                                                          <  3:870960/91‑1 (MQ=255)
                                                                                                                    gCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            >  3:409239/1‑102 (MQ=255)
                                                                                                                    gCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            <  4:409239/102‑1 (MQ=255)
                                                                                                                      acaTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGa                                                      <  1:207902/122‑1 (MQ=255)
                                                                                                                      acaTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGa                                                      >  2:207902/1‑122 (MQ=255)
                                                                                                                      acaTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                         >  3:541224/1‑135 (MQ=255)
                                                                                                                      acaTCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                         <  4:541224/135‑1 (MQ=255)
                                                                                                                         tCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAgg                        >  3:548572/1‑149 (MQ=255)
                                                                                                                            aTAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatg                                            >  3:1453929/1‑126 (MQ=255)
                                                                                                                            aTAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatg                                            <  4:1453929/126‑1 (MQ=255)
                                                                                                                                 gtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                         <  2:1397802/124‑1 (MQ=255)
                                                                                                                                 gtctgtcGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAatgtat                                         >  1:1397802/1‑124 (MQ=255)
                                                                                                                                         cATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTc        >  4:684915/1‑149 (MQ=255)
                                                                                                                                          aTCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCg       <  3:684915/149‑1 (MQ=255)
                                                                                                                                           tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                    >  1:1254785/1‑55 (MQ=255)
                                                                                                                                           tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTc                                                                                                    <  2:1254785/55‑1 (MQ=255)
                                                                                                                                           tCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCgg      >  2:1203435/1‑149 (MQ=255)
                                                                                                                                             ccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaa                                                                                             <  2:615226/60‑1 (MQ=255)
                                                                                                                                             ccGCTCCGCAGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCaataaa                                                                                             >  1:615226/1‑60 (MQ=255)
                                                                                                                                              cGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGagaaagaa                          >  1:181316/1‑126 (MQ=255)
                                                                                                                                              cGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGagaaagaa                          <  2:181316/126‑1 (MQ=255)
                                                                                                                                               gCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            <  3:694663/75‑1 (MQ=255)
                                                                                                                                               gCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATaaaa                                                                            >  4:694663/1‑75 (MQ=255)
                                                                                                                                               gCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGtt  >  1:365866/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGCCGCTATCTCCATCTCCGGGCCTTCATCAAGACTGACAGAAGATCGTTTTGTCAGTCAGGGTGAGCTGGTCAGAGACACCGCCCGCGATATCAGCACGGCGTTGGGACTGAAAGCACATCCATAATGTCTGTCGCATCCCGCTCTGCGGAGCGGGTTTTTTTGACAAAATTTGAAAGTTGGAAAAATTTTCCAATAAATAGAGGTAGGAATAAAATGGCAAAAATGAGAGCCGTTGACGCGGCAATGTATGTGCTGGAGAAAGAAGGTATCACTACCGCCTTCGGTGTT  >  NC_000913/533699‑533990

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: