Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,439 A→T 7.0% F34I (TTT→ATT)  phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,4390AT 7.0% 168.9 / 3.3 57F34I (TTT→ATT) phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base T (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GCGCAGATAGAAACGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGT  >  NC_000913/4315149‑4315729
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                  
gcgcAGATAGAAACGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggtt                                                                                                                                                                                                                                                                                          <  2:202908/301‑4 (MQ=255)
     gttatatccGGTGCGCGTCGTGGCTCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTCCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggttttaac                                                                                                                                                                                                                                                                                     <  2:224228/293‑9 (MQ=255)
       tAGAAACGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACACACGCGAAACGCGCGGAGGTCAAAcc                                                                                                                                                                                                                                                                                   >  1:64484/1‑300 (MQ=255)
          aaaCGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGCTTGGCGTTAAAACCCACGCGAAACGCATGGTGGTCAAACTccc                                                                                                                                                                                                                                                                                >  2:234063/1‑300 (MQ=255)
           aaCGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGc                                                                                                                                                                                                                                                                               >  2:130781/1‑301 (MQ=255)
                  ggggtggtGGCTTTTCACGTTGTTCGAAAGTTCGGTCGTTTCGGCCCCGGCTTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtttt                                                                                                                                                                                                                                                                        <  1:529534/295‑1 (MQ=255)
                     gtcgtGGCGTTTCACGTTGGTGGAATGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                                                                                                                                                                     <  1:172068/301‑1 (MQ=255)
                     gtcgtGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                                                                                                                                                                     <  2:509091/301‑1 (MQ=255)
                       tgtGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTc                                                                                                                                                                                                                                                                   <  2:557423/300‑1 (MQ=255)
                              tttCACTTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATc                                                                                                                                                                                                                                                            <  2:559838/300‑1 (MQ=255)
                                                    gTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgtgt                                                                                                                                                                                                                                      >  2:18961/1‑301 (MQ=255)
                                                            ggCCCCGGCTTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGc                                                                                                                                                                                                                              <  2:64552/301‑1 (MQ=255)
                                                              ccccGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCCCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGt                                                                                                                                                                                                                            >  1:264942/1‑301 (MQ=255)
                                                                                 ttcttctGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACaa                                                                                                                                                                                                         <  2:265210/301‑1 (MQ=255)
                                                                                  tcttctGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCCCACaac                                                                                                                                                                                                        >  1:523663/1‑300 (MQ=255)
                                                                                     tctGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCaa                                                                                                                                                                                                     >  2:335400/1‑300 (MQ=255)
                                                                                       tGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAgga                                                                                                                                                                                                   >  2:363451/1‑300 (MQ=255)
                                                                                        gCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTAATGCGTGGCCGGGCGAAGCTCACAAGCAGGCa                                                                                                                                                                                                  >  1:193151/1‑301 (MQ=255)
                                                                                        ccccTCGCCAGTATCTTACTGCCGGCTGTCAGACCGGGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCa                                                                                                                                                                                                  <  1:130656/300‑1 (MQ=255)
                                                                                          ccACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAACCTCCGCCTGTTTTAGCTCACAAATCCGCCCGCAAACCGCCTCGGTGTCGTACTGCGTGCCCGGGCCACGCCCCCAAACAGGccag                                                                                                                                                                                                >  1:1299/1‑298 (MQ=255)
                                                                                                 attAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACGCCTGAATTGCGCCGGTCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCACGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                                                                                                                                                                         <  1:231504/299‑1 (MQ=255)
                                                                                                   aaatttACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGAGGAAACTGCAAATGCAGGCCGAGCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                                                                                                                                                                       <  1:609593/297‑1 (MQ=255)
                                                                                                     aCTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                                                                                                                                                                     <  2:319015/301‑1 (MQ=255)
                                                                                                       ttACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                                                                                                                                                                   <  2:193394/301‑1 (MQ=255)
                                                                                                            gCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATg                                                                                                                                                                              >  2:20951/1‑301 (MQ=255)
                                                                                                            gCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATg                                                                                                                                                                              >  2:567313/1‑301 (MQ=255)
                                                                                                                  gTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCaa                                                                                                                                                                        >  2:66611/1‑300 (MQ=255)
                                                                                                                        aCCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAACCTGCAAATGCAGGCCGATCATGCCGACCCCTTCGCCATCCAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAACTCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGACGCCCACAACAAAGCATGGTGATTCCCCtacacactgaacagcctcac                                                                                                                                                                  >  1:373617/1‑282 (MQ=255)
                                                                                                                          cgcgcgCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagca                                                                                                                                                                >  2:442362/1‑301 (MQ=255)
                                                                                                                                 ttGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTc                                                                                                                                                         <  2:524335/300‑1 (MQ=255)
                                                                                                                                 cTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAACCTGCAAATGCCGGCCGATCATCCCGACAACTTCGCCATCAACCAGCGCCACATGGTAGCGCATGTTTGGGTCGCGCAGATTCGCGTTAAAACCCACGCGACACGCGTGGTGCTCAAACTCCGCCTGTTCTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTCCTGCGTGGCCGGGCGAAGCCCACAAGCAGCCATGGTGATTCCCCTTATGACAGATACGCGCCAGCACCgctcc                                                                                                                                                         >  1:119448/1‑297 (MQ=255)
                                                                                                                                     ggCATTACCACCAACGCCTGAATTTCGCCGATCCAGTTGACAGGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc                                                                                                                                                     <  1:66521/301‑1 (MQ=255)
                                                                                                                                                caACTCCTGAATTTCGCCGATCCAGTGGACAGGATGCAGATGAAACGGCAAATGCAGGCCGATCATGCCGACAACGTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAg                                                                                                                                          <  1:441798/301‑1 (MQ=255)
                                                                                                                                                   cTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGAGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCCCAAGcc                                                                                                                                       >  1:191169/1‑300 (MQ=255)
                                                                                                                                                           tttCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCAt                                                                                                                               <  2:135751/301‑1 (MQ=255)
                                                                                                                                                                        ggtGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGAAAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATg                                                                                                                  <  1:18910/299‑1 (MQ=255)
                                                                                                                                                                          ggACTTGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGc                                                                                                                <  1:566599/300‑1 (MQ=255)
                                                                                                                                                                               tgatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGACGTCATTGTTGAGCGGATGGCAATcc                                                                                                           >  1:557586/1‑301 (MQ=255)
                                                                                                                                                                                      ctacacgacgctctgccgatcgtaTCCTGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGGTTGGGTCGCTCAGATTGGCGTTAAAACCAACGCGAAACGCGTGGTGGTCCAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAg                                                                                                    <  1:569644/278‑1 (MQ=255)
                                                                                                                                                                                            aCTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAga                                                                                              >  2:547810/1‑300 (MQ=255)
                                                                                                                                                                                              tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAgcgg                                                                                            >  2:298395/1‑301 (MQ=255)
                                                                                                                                                                                                      cAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAACTCAGCGCGTACACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGAAGCGTGTCGACCGACTGGCGCAAGCTGCCGTAATTGTTGAGCGTATGGCAATCCTGGGGAGTATAGCGGGCGgcgcc                                                                                    >  1:468989/1‑300 (MQ=255)
                                                                                                                                                                                                           cGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCa                                                                               >  2:557499/1‑301 (MQ=255)
                                                                                                                                                                                                            gATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGGGGAgagatcggaagcgcacacgtctg                                                                              >  2:570358/1‑279 (MQ=255)
                                                                                                                                                                                                                 ccaaacatgcgctaccatctggcgctgcttgATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTACACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACACGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCCTCAGCAGCGTGTCGACCGACTGGCGCCGGCTGCCGTCATCGTTGAGCGTATGGCAACCCTGTGGAGAATACCGGGCGGCCCGCCCCCggcggg                                                                         >  1:610303/31‑301 (MQ=255)
                                                                                                                                                                                                                         aCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAAtt                                                                 <  2:1299/301‑1 (MQ=255)
                                                                                                                                                                                                                           ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGACCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCCGGCTGCCCTCATTGTTCAGCGTATTGCACTCCTGTGGAGTATAGCGGGCGGCCCGCGCCCGGCGGGCGTGACTTTc                                                               >  1:162213/1‑301 (MQ=255)
                                                                                                                                                                                                                               ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAATATAGCGGGAGGCGCGGGCCCGGGGGGCGTGAAGTTCCCtt                                                           >  1:570883/1‑300 (MQ=255)
                                                                                                                                                                                                                                cttCAAGCAGGGCCGGTTGGTGGGGCAGGGTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTGAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTgg                                                          <  1:556785/299‑1 (MQ=255)
                                                                                                                                                                                                                                    aaGCAGCGCCAGTTGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGtt                                                      <  2:191412/301‑1 (MQ=255)
                                                                                                                                                                                                                                            ccggATGGTAGCGCGTGTTTGGGTGGCTCAGATTGTCGTTAAAACCCCAGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGTTTTCACGGc                                              <  1:547137/298‑1 (MQ=255)
                                                                                                                                                                                                                                                         cATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGCGGAGTAGAGCTGGCGGCGACCGCCAGGCGGGCGTTAACTTCACTGACGCTTTCCCGCCCAccgccctccac                                 >  1:52771/1‑291 (MQ=255)
                                                                                                                                                                                                                                                                          gATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTCATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCCGGCGGCAGAtcc                >  1:196107/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                cGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCg          <  2:571597/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ccaccACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATCGTTGAGCGTATGGCAATCCTGTGGAGAATCGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTCGCGTTTTCACGGCCACGGTTTTCCAGGCACCGGCGGAGGACCTCCGCCGAAcccc    >  1:194510/4‑300 (MQ=255)
                                                                                                                                                                                                                                                                                        cccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCAGGCGAAAACTGt  >  2:65740/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        accACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTCTTAGCCCACAAATCAGCGCGTAAACCGCGTCGGTGCCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAACCCTGTGGAGTATAGCGCGCGCAGCGCGCCAGGCGGGCGTTAATTTCACAGGCGTTTTCACGGCCACGCTTTTCCAGGCGCGGGCGGAGGATCCCCGGCCAAccccgc  >  1:362492/2‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                  
GCGCAGATAGAAACGGTGCGCGTCGTGGCGCTTCACGTTGGTCGAAAGTTCGGTCATTTCGGCCCCGGCCTGGCGGGCTTCTTCTTCTGCCCACGCCAGTAACTTACTGCCGACGTTCAGACCGCGCGCCTGCGGCATTACCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGCGCCAGGCGGGCGTTAATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGT  >  NC_000913/4315149‑4315729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: