Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,554,523 T→G 10.8% intergenic (+162/‑53) 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,5230TG10.8% 232.2 / 17.9 74intergenic (+162/‑53)uxuB/uxuRD‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base G (5/3);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAAC  >  NC_000913/4554375‑4554669
                                                                                                                                                    |                                                                                                                                                  
cTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGtt                                                                                                                                                    <  4:687114/149‑1 (MQ=255)
  aCAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                                  <  4:178594/149‑1 (MQ=255)
        ttACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                            >  3:236676/1‑149 (MQ=255)
          aCTGGTCAATCAAACTGATATTTGGTTGCCCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGCATGCGGCGTGAAcgccttat                                                                                                                                          <  2:319829/149‑5 (MQ=255)
            tGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGACAAGGCGTTCACGCCGCATcc                                                                                                                                        <  1:551132/149‑1 (MQ=255)
              gTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGTCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                      >  2:418343/1‑149 (MQ=255)
                      aaCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAtt                                                                                                                              >  1:317797/1‑149 (MQ=255)
                         tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                          >  1:376453/1‑130 (MQ=255)
                          gATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAt                                                                                                                          >  1:372427/1‑149 (MQ=255)
                            tatTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTa                                                                                                                        >  4:332814/1‑149 (MQ=255)
                             atTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTaa                                                                                                                       <  2:616685/149‑1 (MQ=255)
                             atTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTaa                                                                                                                       <  2:355263/149‑1 (MQ=255)
                              tttGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                               <  3:47605/124‑1 (MQ=255)
                              tttGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                                      >  2:104645/1‑149 (MQ=255)
                              tttGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAAc                                                                                                                      <  4:501300/149‑1 (MQ=255)
                              tttGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAAGGGGATAGTGTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                               >  4:47605/1‑124 (MQ=255)
                                  gTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCa                                                                                                                  >  4:304945/1‑149 (MQ=255)
                                  gTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCa                                                                                                                  >  4:211191/1‑149 (MQ=255)
                                        cAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                                                                            >  3:99715/1‑149 (MQ=255)
                                         aGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACg                                                                                                           >  4:483759/1‑149 (MQ=255)
                                             ttCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGaca                                                                                                       <  1:138561/149‑1 (MQ=255)
                                                 tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                          >  3:90739/1‑106 (MQ=255)
                                                 tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                          <  4:90739/110‑5 (MQ=255)
                                                 tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                     >  1:402205/1‑115 (MQ=38)
                                                 tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                     <  2:402205/115‑1 (MQ=38)
                                                  ttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                 <  3:75938/102‑1 (MQ=255)
                                                  ttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                 >  4:75938/1‑102 (MQ=255)
                                                     tttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCa                                                                                               <  4:236676/149‑1 (MQ=255)
                                                          ccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaaa                                                                                          >  1:402575/1‑149 (MQ=255)
                                                            aCCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATc                                                                                        >  3:183236/1‑149 (MQ=255)
                                                             ccTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                                                                       <  1:282554/149‑1 (MQ=255)
                                                             ccTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                                                                       >  1:504628/1‑149 (MQ=255)
                                                              cTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTg                                                                                      >  4:109189/1‑149 (MQ=255)
                                                               tGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                                                                                     <  3:483759/149‑1 (MQ=255)
                                                               tGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                                                                                     <  4:99715/149‑1 (MQ=255)
                                                               tGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                                                                                     <  3:604901/149‑1 (MQ=255)
                                                               tGTACGTGCCAACTTCCAGCGCTAATGGTATAGTTTGAGAGTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCAGGAAATCTGc                                                                                     <  4:567267/149‑1 (MQ=255)
                                                                    gTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACct                                                                                >  4:553240/1‑149 (MQ=255)
                                                                        cAACTTCCAGTGCTAATGGTATAGTTTGAGATTGACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTgc                                                                            >  3:437582/1‑149 (MQ=255)
                                                                        cAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTgc                                                                            <  2:372427/149‑1 (MQ=255)
                                                                         aaCTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTgcg                                                                           >  3:259184/1‑149 (MQ=255)
                                                                            ttCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGa                                                                                                                           >  1:579759/1‑98 (MQ=255)
                                                                            ttCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGa                                                                                                                           <  2:579759/98‑1 (MQ=255)
                                                                                 gTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGAcc                                                                   >  2:397376/1‑149 (MQ=255)
                                                                                    cTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTa                                                                >  1:425981/1‑149 (MQ=255)
                                                                                      aaTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTAcc                                                              >  3:685697/1‑149 (MQ=255)
                                                                                       aTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGaca                                                                                                       >  4:636142/1‑107 (MQ=255)
                                                                                       aTGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGaca                                                                                                       <  3:636142/107‑1 (MQ=255)
                                                                                         ggTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTa                                                                >  2:337587/1‑144 (MQ=255)
                                                                                         ggTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCGTa                                                                <  1:337587/144‑1 (MQ=255)
                                                                                           tataGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGaa                                                         >  1:540316/1‑149 (MQ=255)
                                                                                                tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACct                                                                                >  2:515691/1‑121 (MQ=255)
                                                                                                tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACct                                                                                <  1:515691/121‑1 (MQ=255)
                                                                                                tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCgg                                                    >  1:367629/1‑149 (MQ=255)
                                                                                                 ttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCggg                                                   <  1:418370/149‑1 (MQ=255)
                                                                                                     gaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCga                                               <  2:317797/149‑1 (MQ=255)
                                                                                                       ttAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatg                                             >  3:441860/1‑149 (MQ=255)
                                                                                                         aaCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatgat                                           <  2:402575/149‑1 (MQ=255)
                                                                                                         aaCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatgat                                           <  3:59386/149‑1 (MQ=255)
                                                                                                                    tAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGa                                <  3:304945/149‑1 (MQ=255)
                                                                                                                    tAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGa                                >  4:29124/1‑149 (MQ=255)
                                                                                                                     aaaaTTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGAt                               <  4:685697/149‑1 (MQ=255)
                                                                                                                       aaTTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCa                             <  3:553240/149‑1 (MQ=255)
                                                                                                                           gCCCGTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATaaa                         <  4:437582/149‑1 (MQ=255)
                                                                                                                            cccGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAg                        >  3:160194/1‑149 (MQ=255)
                                                                                                                               gTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACg                     >  1:271444/1‑149 (MQ=255)
                                                                                                                                ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGc                    <  3:29124/149‑1 (MQ=255)
                                                                                                                                 tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACTCAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGTGCGATGATCCGCGATCTGATCATAAAGACGcc                   <  4:441860/149‑1 (MQ=255)
                                                                                                                                         ggATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCAGGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATc           <  4:259184/149‑1 (MQ=255)
                                                                                                                                          gatgcggCGTGAACGCCTTATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATcc          >  4:415440/6‑149 (MQ=255)
                                                                                                                                          gatgcggCGTGAACGCCTTATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATcc          >  2:334876/6‑149 (MQ=255)
                                                                                                                                           aTAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCt         >  2:312549/1‑149 (MQ=255)
                                                                                                                                           aTAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCt         <  3:211191/149‑1 (MQ=255)
                                                                                                                                              aGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGc      <  1:418343/149‑1 (MQ=255)
                                                                                                                                                  gTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAAc  >  4:223036/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCTGGCGAAC  >  NC_000913/4554375‑4554669

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: