Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,554,524 C→A 5.7% intergenic (+163/‑52) uxuB → / → uxuR D‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,554,5240CA 5.7% 295.3 / 7.5 89intergenic (+163/‑52)uxuB/uxuRD‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base A (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

TACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAACG  >  NC_000913/4554376‑4554670
                                                                                                                                                    |                                                                                                                                                  
tACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                    >  3:369788/1‑149 (MQ=255)
tACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                    >  3:779091/1‑149 (MQ=255)
    gACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                >  4:244266/1‑149 (MQ=255)
         aCTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                           <  4:752674/149‑5 (MQ=255)
           tGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                         >  3:167016/1‑149 (MQ=255)
           tGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                         <  2:497520/149‑1 (MQ=255)
             gTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                       <  4:369788/149‑1 (MQ=255)
              tcaatcaaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                      >  3:382200/1‑149 (MQ=255)
                aatcaaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                                    >  4:258520/1‑149 (MQ=255)
                     aaCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAtt                                                                                                                               >  4:632658/1‑149 (MQ=255)
                      aCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAAttt                                                                                                                              >  3:509897/1‑149 (MQ=255)
                      aCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                              <  3:760360/149‑1 (MQ=255)
                        tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGa                                                                                                                            <  3:258520/149‑1 (MQ=255)
                        tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  1:344646/1‑139 (MQ=255)
                        tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     <  2:344646/140‑2 (MQ=255)
                        tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                      <  3:362566/139‑1 (MQ=255)
                        tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                      >  4:362566/1‑139 (MQ=255)
                            atTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTaa                                                                                                                        >  2:737482/1‑149 (MQ=255)
                             tttGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                                       >  4:276058/1‑149 (MQ=255)
                                     aCCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTa                                                                                                               >  3:603753/1‑149 (MQ=255)
                                       cAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                                                                             >  2:643139/1‑149 (MQ=255)
                                        aGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACg                                                                                                            <  4:622093/149‑1 (MQ=255)
                                        aGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACg                                                                                                            >  1:62470/1‑149 (MQ=255)
                                               gTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaa                                                                                                     <  3:451984/149‑1 (MQ=255)
                                               gTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaa                                                                                                     >  2:132058/1‑149 (MQ=255)
                                                tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAt                                                                                                                           <  4:595434/126‑1 (MQ=255)
                                                tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAt                                                                                                                           >  3:595434/1‑126 (MQ=255)
                                                tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaac                                                                                                    >  2:592178/1‑149 (MQ=255)
                                                tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaac                                                                                                    >  2:104497/1‑149 (MQ=255)
                                                     ttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCAt                                                                                               <  2:62470/149‑1 (MQ=255)
                                                        cccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                              >  4:331618/1‑115 (MQ=255)
                                                        cccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                              <  3:331618/115‑1 (MQ=255)
                                                         ccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaaa                                                                                           >  2:375433/1‑149 (MQ=255)
                                                         ccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaaa                                                                                           >  3:503691/1‑149 (MQ=255)
                                                            ccTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                                                                        <  4:503691/149‑1 (MQ=255)
                                                              tGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                                                                                      <  1:132058/149‑1 (MQ=255)
                                                              tGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                                                                                      <  4:382200/149‑1 (MQ=255)
                                                               gTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                >  4:72860/1‑90 (MQ=255)
                                                               gTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                                <  3:72860/90‑1 (MQ=255)
                                                               gTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGcc                                                                                     >  4:291350/1‑149 (MQ=255)
                                                                   gTGCCAACTTCCTGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACct                                                                                 <  2:612753/149‑1 (MQ=255)
                                                                    tGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctc                                                                                <  1:768944/149‑1 (MQ=255)
                                                                     gCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                                                                               >  2:208749/1‑149 (MQ=255)
                                                                     gCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                                                                               >  2:22539/1‑149 (MQ=255)
                                                                      ccAACTTCCAGTTCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTg                                                                              <  1:643139/149‑1 (MQ=255)
                                                                        aaCTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTgcg                                                                            >  1:638330/1‑149 (MQ=255)
                                                                            tCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCaaa                                                                        <  1:592178/149‑1 (MQ=255)
                                                                            tCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCaaa                                                                        <  1:737482/149‑1 (MQ=255)
                                                                               aGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGAc                                                                     >  2:302684/1‑149 (MQ=255)
                                                                               aGTGCTAAGGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACTGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGAc                                                                     >  4:447235/1‑149 (MQ=255)
                                                                                     aaTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                                       <  3:2199/93‑1 (MQ=255)
                                                                                     aaTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                                       >  4:2199/1‑93 (MQ=255)
                                                                                      aTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCa                                                              <  2:251666/149‑1 (MQ=255)
                                                                                               tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCgg                                                     >  4:337796/1‑149 (MQ=255)
                                                                                                  gagaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGc                                                  <  3:276058/149‑1 (MQ=255)
                                                                                                    gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                                <  3:244266/149‑1 (MQ=255)
                                                                                                    gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                                >  3:402316/1‑149 (MQ=255)
                                                                                                    gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                                >  2:737828/1‑149 (MQ=255)
                                                                                                      ttagggggagCCGCATAATTGGCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatg                                              >  2:532851/6‑149 (MQ=38)
                                                                                                       tAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatga                                             <  3:447235/149‑1 (MQ=255)
                                                                                                         aCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                                                                             >  2:778435/1‑83 (MQ=255)
                                                                                                         aCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                                                                             <  1:778435/83‑1 (MQ=255)
                                                                                                         aCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATc                                           <  2:638330/149‑1 (MQ=255)
                                                                                                           gggggCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTa                                                                 >  4:185179/1‑125 (MQ=255)
                                                                                                           gggggCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTa                                                                 <  3:185179/125‑1 (MQ=255)
                                                                                                           gggggCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCcg                                         >  4:43921/1‑149 (MQ=255)
                                                                                                           gggggCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCcg                                         >  4:586271/1‑149 (MQ=255)
                                                                                                               gCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGAt                                     <  1:737828/149‑1 (MQ=255)
                                                                                                                   tAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGa                                 <  4:167016/149‑1 (MQ=255)
                                                                                                                   tAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGa                                 <  1:718759/149‑1 (MQ=255)
                                                                                                                    aaaaTTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGAt                                <  3:632658/149‑1 (MQ=255)
                                                                                                                       aTTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCAt                             >  4:100362/1‑149 (MQ=255)
                                                                                                                             cGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGAc                       >  4:237888/1‑149 (MQ=255)
                                                                                                                              gTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACg                      >  4:770003/1‑149 (MQ=255)
                                                                                                                                tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCCAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGcc                    <  4:509897/149‑1 (MQ=255)
                                                                                                                                tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATAAACCGCACCTAACCGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGcc                    <  4:779091/149‑1 (MQ=255)
                                                                                                                                   aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTa                 <  4:402316/149‑1 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACTAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACa               >  3:366181/1‑149 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatgat                                            >  2:144482/1‑120 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatgat                                            <  1:144482/120‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCt          <  3:337796/149‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCt          >  3:249574/1‑149 (MQ=255)
                                                                                                                                              ggCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                                                                        <  2:310297/67‑1 (MQ=255)
                                                                                                                                              ggCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                                                                        >  1:310297/1‑67 (MQ=255)
                                                                                                                                              ggCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTa                 <  3:739101/138‑1 (MQ=255)
                                                                                                                                              ggCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTa                 >  4:739101/1‑138 (MQ=255)
                                                                                                                                                 gTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCt                                   >  2:468571/1‑117 (MQ=255)
                                                                                                                                                 gTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCt                                   <  1:468571/117‑1 (MQ=255)
                                                                                                                                                 gTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAAc   >  4:101750/1‑149 (MQ=255)
                                                                                                                                                 gTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGACCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAAc   <  4:603753/149‑1 (MQ=255)
                                                                                                                                                  ttCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAACg  <  1:532851/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAACG  >  NC_000913/4554376‑4554670

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: