Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,318,215 G→A L184L (CTG→TTG)  phnK ← carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,318,2150GA98.3% 186.2 / ‑5.4 59L184L (CTG→TTG) phnKcarbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (34/24);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCC  >  NC_000913/4317927‑4318510
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                       
aCGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCa                                                                                                                                                                                                                                                                                              >  1:423771/1‑300 (MQ=255)
    tGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCAc                                                                                                                                                                                                                                                                                             <  2:71860/297‑1 (MQ=255)
      aCGCTTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGaca                                                                                                                                                                                                                                                                                        <  1:116237/300‑1 (MQ=255)
         cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGCGCACGCGGTCGGTTCACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCATGCACCGACACAt                                                                                                                                                                                                                                                                                     >  2:118361/1‑300 (MQ=255)
          gTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGNGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATc                                                                                                                                                                                                                                                                                    >  1:133218/1‑300 (MQ=255)
                   gACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCccgcc                                                                                                                                                                                                                                                                           >  1:431220/1‑300 (MQ=255)
                     ccTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCGATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCAACAGCAAACGGTCCGCCGGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCATGCACCGACACATCCAGCccgc                                                                                                                                                                                                                                                                            >  2:274563/1‑297 (MQ=255)
                      cctgcgagatgtgtataagagacaggATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCcggt                                                                                                                                                                                                                                                                        >  1:686255/25‑300 (MQ=255)
                                  cATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACAANCGCGAGGTTCAGCTCCACCAACAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATTCAGACCGCCGGTAGGTTCATCCAt                                                                                                                                                                                                                                                            >  2:554731/1‑300 (MQ=255)
                                               gATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAATGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTc                                                                                                                                                                                                                                               >  1:186896/1‑300 (MQ=255)
                                               gATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATaaa                                                                                                                                                                                                                                                         <  1:196177/290‑1 (MQ=255)
                                               gATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGAGCAGCAGGTCGAGCAAGCGGGACTGTACCGACACATCCAGCCCGCCGGTAGGTTCATCCATaaa                                                                                                                                                                                                                                                         >  2:196177/1‑290 (MQ=255)
                                                  aaGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGa                                                                                                                                                                                                                                            <  1:354480/300‑1 (MQ=255)
                                                                gcATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATTGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAACAGGtt                                                                                                                                                                                                                              >  1:639662/1‑300 (MQ=255)
                                                                  aTCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGTTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACCATCACCACCGCGAGGTTCAGCTACACCAACAGGCCGCGCAGCAGGTCGAGCAAGCGGGCATGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCCTAAACACCAGCTTCGGATGCGTCAAAAGGTTGc                                                                                                                                                                                                                            >  2:170132/1‑300 (MQ=255)
                                                                  aTCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGc                                                                                                                                                                                                                            >  1:696153/1‑300 (MQ=255)
                                                                                   ttCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgc                                                                                                                                                                                                           <  1:6136/300‑1 (MQ=255)
                                                                                                       cagcagCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaa                                                                                                                                                                                       >  1:666863/1‑300 (MQ=255)
                                                                                                       cagcagCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAACAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGaaaa                                                                                                                                                                                       >  1:355377/1‑300 (MQ=255)
                                                                                                               gtgtATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATCAACACCAGCTTCGGATGCGTCAACAGGTTGTGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAggt                                                                                                                                                                                    >  2:386698/1‑295 (MQ=255)
                                                                                                                 gtATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGCCGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGc                                                                                                                                                                             >  2:608133/1‑300 (MQ=255)
                                                                                                                 gtATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCAGGAAAAGGTGGTAGGc                                                                                                                                                                             >  2:286137/1‑300 (MQ=255)
                                                                                                                 gtATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTAGGATGCGTCAACAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGc                                                                                                                                                                             >  2:680653/1‑300 (MQ=255)
                                                                                                                                   tcgtcgAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTg                                                                                                                                                           <  1:42859/300‑1 (MQ=255)
                                                                                                                                                    ggttaaatcACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAACCGGTCCACCAACAAGCGGGCGACGCCTAAATCATGGGTGACAATCAACACCGCGAGGTTCAGCTCCCCCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGTACCGACACATCCCTCCCGACGGTCTGTTCCTCCCTAAACAACAGCTTCGGATGCGTCACAAGGTTGCGGTCAATCTGCAAACGCTGTTGCATACCGCCGGAAAAGGTGGTAGGCAGGTCGTCGATCAGGTTGGTGGGAAtctaaacctc                                                                                                                                          >  2:135713/9‑293 (MQ=255)
                                                                                                                                                    ggtTAACCCACTCTCCACCAATTGCACCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAACTCCACCACCAGGCCGCGCAGCAGGTTGAGCAGGCGGGCCTACACCGACACACCCAGcaagcaggtc                                                                                                                                                                                                                                                                       >  2:135877/1‑166 (MQ=255)
                                                                                                                                                        aaCCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                                                                                                                      <  1:286989/300‑1 (MQ=255)
                                                                                                                                                                ctcCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACtt                                                                                                                              >  1:103016/1‑300 (MQ=255)
                                                                                                                                                                 tcCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTc                                                                                                                             <  1:399370/300‑1 (MQ=255)
                                                                                                                                                                       acTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGTAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGTGCAATCTACAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCAGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcgc                                                                                                                        <  1:18352/299‑1 (MQ=255)
                                                                                                                                                                       acTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTgcgc                                                                                                                        <  1:18496/299‑1 (MQ=255)
                                                                                                                                                                            ccccTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGc                                                                                                                  >  1:30179/1‑300 (MQ=255)
                                                                                                                                                                                          ccAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGt                                                                                                    >  1:227025/1‑300 (MQ=255)
                                                                                                                                                                                          ccAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCCCCCCCGGTGGCGGGAATCTCAAAATATTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGt                                                                                                    >  1:227083/1‑300 (MQ=255)
                                                                                                                                                                                          ccAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCCATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGcagt                                                                                                    >  1:227097/1‑297 (MQ=255)
                                                                                                                                                                                                 aCGGTCCACCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGAGGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGGTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACg                                                                                             <  2:326629/300‑1 (MQ=255)
                                                                                                                                                                                                           agcagGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTc                                                                                   <  2:686255/300‑1 (MQ=255)
                                                                                                                                                                                                                    gcgACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAACAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGc                                                                          <  2:587682/300‑1 (MQ=255)
                                                                                                                                                                                                                    gcgACGCCTAAATCAGGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGc                                                                          <  1:554536/300‑1 (MQ=255)
                                                                                                                                                                                                                        cGCCTAAATCATGGGTNACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATANCGCCGTAATGACGTGCCCNGGTCGCCATCAGCCGCt                                                                      >  1:586602/1‑300 (MQ=255)
                                                                                                                                                                                                                        cGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCt                                                                      >  2:584172/1‑300 (MQ=255)
                                                                                                                                                                                                                          ccTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTc                                                                     >  1:117210/1‑299 (MQ=255)
                                                                                                                                                                                                                             aaaTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGcc                                                                  >  1:223569/1‑299 (MQ=255)
                                                                                                                                                                                                                               aTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGat                                                               >  2:9730/1‑300 (MQ=255)
                                                                                                                                                                                                                                       tGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGTTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTgccg                                                        <  2:296186/299‑1 (MQ=255)
                                                                                                                                                                                                                                         aCAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCTACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCt                                                     >  1:177778/1‑300 (MQ=255)
                                                                                                                                                                                                                                         aCAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCt                                                     >  1:376893/1‑300 (MQ=255)
                                                                                                                                                                                                                                                ccaccGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGaca                                              >  1:355075/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        aGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTggc                                        >  2:246614/1‑298 (MQ=255)
                                                                                                                                                                                                                                                         ggTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGACGCTCGCCGATATTGCCGCCTGCCGACACCTggcggc                                     <  1:97777/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccaccaGGCAGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCAGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGACTGCCGACACCTGGCGGCGCAGGCCGTc                           >  2:13678/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                               gcgcAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGc                 <  1:246614/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gcgcAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCt                <  2:423771/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gcgcAGCAGGTCGAGCAAGCGGGCCTGCACCGACACACCCAGCCCGCCGGTCGGTGCATCCAGAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGc                 <  2:30179/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      agGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGcac         <  1:30935/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      agGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCAGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGcac         <  1:446288/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ggTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGTTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACTAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGTGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGcacca       <  2:668564/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            aGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGTTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGccc  <  1:123101/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            aGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGccc  <  1:123132/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                       
ACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTGCCGACACCTGGCGGCGCAGGCCGTCGAGTGGATGCTGATGCACCACGCCC  >  NC_000913/4317927‑4318510

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