Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,554,517 A→C 5.7% intergenic (+156/‑59) 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,5170AC 5.7% 309.7 / 10.3 88intergenic (+156/‑59)uxuB/uxuRD‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (44/39);  new base C (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

ATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCT  >  NC_000913/4554369‑4554662
                                                                                                                                                    |                                                                                                                                                 
aTTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATaa                                                                                                                                                   <  1:739989/149‑1 (MQ=255)
 ttACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCAGGATAAg                                                                                                                                                  <  2:441028/149‑1 (MQ=255)
   aCCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGc                                                                                                                                                <  1:669070/149‑1 (MQ=255)
     ccTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGt                                                                                                                                              >  2:421576/1‑149 (MQ=255)
      cTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGtt                                                                                                                                             <  1:421576/149‑1 (MQ=255)
       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                                                                                                                                 <  2:497520/149‑1 (MQ=255)
                  tGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                 >  3:167016/1‑149 (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                                                                                                                             <  2:344646/140‑2 (MQ=255)
                               tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                             >  1:344646/1‑139 (MQ=255)
                               tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                              >  4:362566/1‑139 (MQ=255)
                               tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                              <  3:362566/139‑1 (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)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGc                                                                                                                                                >  4:101454/1‑97 (MQ=255)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGc                                                                                                                                                <  3:101454/97‑1 (MQ=255)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAt                                                                                                                   >  3:595434/1‑126 (MQ=255)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAt                                                                                                                   <  4:595434/126‑1 (MQ=255)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaac                                                                                            >  2:104497/1‑149 (MQ=255)
                                                       tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaac                                                                                            >  2:592178/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                                                                                   >  3:503691/1‑149 (MQ=255)
                                                                ccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaaa                                                                                   >  2:375433/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:22539/1‑149 (MQ=255)
                                                                            gCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                                                                       >  2:208749/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                                                                                                               >  4:2199/1‑93 (MQ=255)
                                                                                            aaTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                               <  3:2199/93‑1 (MQ=255)
                                                                                             aTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCa                                                      <  2:251666/149‑1 (MQ=255)
                                                                                                      tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCgg                                             >  4:337796/1‑149 (MQ=255)
                                                                                                         gagaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGc                                          <  3:276058/149‑1 (MQ=255)
                                                                                                           gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                        >  2:737828/1‑149 (MQ=255)
                                                                                                           gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                        >  3:402316/1‑149 (MQ=255)
                                                                                                           gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                        <  3:244266/149‑1 (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                         <  1:718759/149‑1 (MQ=255)
                                                                                                                          tAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGa                         <  4:167016/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)
                                                                                                                                                    |                                                                                                                                                 
ATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCT  >  NC_000913/4554369‑4554662

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