Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,322,786 G→A intergenic (‑364/+550) phnF ← / ← phnD putative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,322,7860GA98.3% 192.7 / NA 58intergenic (‑364/+550)phnF/phnDputative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/31);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCG  >  NC_000913/4322492‑4323078
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
accttGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGatggtg                                                                                                                                                                                                                                                                                                 <  1:184009/296‑1 (MQ=255)
     gACGATGATAAGCACCATCAGGGCGCAGTTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCGTGGAGAGCACGCCGATGGTGTGGAt                                                                                                                                                                                                                                                                                            <  2:240753/300‑1 (MQ=255)
     gACGATCATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCNGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAANGCAGCTTGGAGAGCACGCCGATGGTGTGGAt                                                                                                                                                                                                                                                                                            <  2:626983/300‑1 (MQ=255)
      aCGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATa                                                                                                                                                                                                                                                                                           <  1:154845/300‑1 (MQ=255)
       cGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTNCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCANTCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATaa                                                                                                                                                                                                                                                                                          >  2:804471/1‑300 (MQ=255)
         atgatAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGNGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAAc                                                                                                                                                                                                                                                                                        <  1:647897/300‑1 (MQ=255)
           gatAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAg                                                                                                                                                                                                                                                                                      <  2:512633/300‑1 (MQ=255)
              gcgagatgtgtataagagacagtgttgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTCGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCATTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAgcgc                                                                                                                                                                                                                                                                                   >  2:745665/23‑300 (MQ=255)
               aGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgcc                                                                                                                                                                                                                                                                                  >  1:678682/1‑300 (MQ=255)
                 cACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGAACATGCAGACGACGGTCGCAGAGCGGACGTTGGATTCGAAGCGATAGAGGGCGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTCCGAGCTTGTTGGCACCGGTGTCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGATTCCACCGCTTCGGAACGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAGCAGcgcccg                                                                                                                                                                                                                                                                                >  2:77964/1‑298 (MQ=255)
                   ccATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgccagca                                                                                                                                                                                                                                                                              >  1:392865/1‑300 (MQ=255)
                            cgcAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCgccagcacgccagcac                                                                                                                                                                                                                                                                     <  1:522665/300‑1 (MQ=255)
                                           aCTGGAAACCGCGAATCGCTTCCCACGGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCACCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCTGAGAGCAACTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGggcc                                                                                                                                                                                                                                                      <  2:339029/300‑1 (MQ=255)
                                           aCTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGTAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGggcc                                                                                                                                                                                                                                                      <  1:316260/300‑1 (MQ=255)
                                           aCTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCAAGCACGCCAGCACGCCAGCGAACGggcc                                                                                                                                                                                                                                                      >  2:792398/1‑300 (MQ=255)
                                                            gCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGTGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCAGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTCGGAGGGCACGCAGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGAAGACGACCa                                                                                                                                                                                                                                     >  2:307978/1‑300 (MQ=255)
                                                            gCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTNAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCa                                                                                                                                                                                                                                     >  1:648539/1‑300 (MQ=255)
                                                                 ccACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAAc                                                                                                                                                                                                                                <  1:127482/300‑1 (MQ=255)
                                                                      gggTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGAATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCAACTGTGGCAGCACGCCGTAGAGGATCTCTTAGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCAACGGGCCGGGCTCCATCGCTTCTACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAATAGCGCCAGCATGCCAGCATGACAGCGAACGGGCGGAGGACGGCGGCGACCACGAACAACAt                                                                                                                                                                                                                           >  2:90630/1‑300 (MQ=255)
                                                                         tGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATATAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGTAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGc                                                                                                                                                                                                                        <  2:506371/300‑1 (MQ=255)
                                                                         tGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGc                                                                                                                                                                                                                        <  1:126953/300‑1 (MQ=255)
                                                                          gACGCCGATCCTGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATATAGGGAGTAGTAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCTGAGAGGATCTCTTCGAGCTTGTTGGCACCGGTTGCGCGAATGCCTTCCACCGGGCCGTGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGTAGAGCACGCCGATGGTGTGGATAAGCAGCGCCAGCACGCCAGCATGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCg                                                                                                                                                                                                                       <  1:7072/300‑1 (MQ=255)
                                                                           aCGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCNGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGa                                                                                                                                                                                                                      <  2:762821/300‑1 (MQ=255)
                                                                                gATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCATTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGAcc                                                                                                                                                                                                                 <  2:71162/300‑1 (MQ=255)
                                                                                    ccGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCAttt                                                                                                                                                                                                             <  2:513681/300‑1 (MQ=255)
                                                                                     cGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGNAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTc                                                                                                                                                                                                            >  1:646947/1‑300 (MQ=255)
                                                                                     cGCCTGCGCCGACCATGCCGACGACGGTCGACGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCAGTAGAGGATCTCTTCGAGCTTNTTGGCCCCGGTGGCGCGAATGCCTTCNACCGGGACGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTc                                                                                                                                                                                                            >  2:595457/1‑300 (MQ=255)
                                                                                                                  gCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAggc                                                                                                                                                                               <  1:552598/300‑1 (MQ=255)
                                                                                                                          gACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCACTTCGTTAATGGCGCGGCAGGCGTCCATCAGACGGCGAACg                                                                                                                                                                       >  2:749688/1‑300 (MQ=255)
                                                                                                                                       aatcGATAGAGGGAGTAGGAGATCAGCATTTGCATCACCTTTTGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGTTGGCGCGAATGCCTTCCGCCGGGCCGGGCTCAATCGCTTCCACCGCTGCGGAAAGCAGCTTGGAGAGCACGTCGATGGTGTGGAGAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAAccc                                                                                                                                                          <  1:707687/297‑1 (MQ=255)
                                                                                                                                            aTAGAGGGAGTAGGAGATCAGCAGTGGCATCAACTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGCGGCGCGAATACCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGATGTGGATAAACAGTGCCAGCACGACAGCACGCCAGCGAACGGGCCGAGGCCGACGGCAACCACGAACAGCACGGCGAAGACCATTTCGTTAATAGCACGGCAGGCGTCCATCCGGAGGCAAACGGGATGGCAAACCcaccac                                                                                                                                                     >  2:12666/1‑300 (MQ=255)
                                                                                                                                                 gggAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGAACATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAAACCAACACggaac                                                                                                                                                >  2:664854/1‑297 (MQ=255)
                                                                                                                                                      tAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGt                                                                                                                                           >  1:803475/1‑300 (MQ=255)
                                                                                                                                                          agaTCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTc                                                                                                                                       <  2:732303/300‑1 (MQ=255)
                                                                                                                                                           gccccgcagTGGCATCACCTGCGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAGACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGCCCCCGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCAGCAGGCGTACATAAGGCGGCGAACGGGCTGGTAAAACCCCAACGGAACAAGGTTTTCg                                                                                                                                      >  2:425218/6‑300 (MQ=255)
                                                                                                                                                           gaTCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGTGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCCCCACGGCACCAGGTTTTCg                                                                                                                                      >  2:834070/1‑300 (MQ=255)
                                                                                                                                                            aTCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCgg                                                                                                                                     <  1:529610/300‑1 (MQ=255)
                                                                                                                                                                gcagTGGCATCACCTGTGGCAGCACGCAGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGATTTTCGGCGCt                                                                                                                                 >  2:465219/1‑300 (MQ=255)
                                                                                                                                                                       cATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAgg                                                                                                                          >  1:871171/1‑300 (MQ=255)
                                                                                                                                                                         tCACCTGTGGCAGCACGCCTTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGcc                                                                                                                        <  2:118059/300‑1 (MQ=255)
                                                                                                                                                                         tCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGcc                                                                                                                        <  1:381778/300‑1 (MQ=255)
                                                                                                                                                                              tgtgGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAgg                                                                                                                   <  1:330390/300‑1 (MQ=255)
                                                                                                                                                                                      caCGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGgag                                                                                                           <  1:542614/300‑1 (MQ=255)
                                                                                                                                                                                      caCGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGgag                                                                                                           <  2:792765/300‑1 (MQ=255)
                                                                                                                                                                                       aCGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGgaga                                                                                                          <  2:830734/300‑1 (MQ=255)
                                                                                                                                                                                          ccatAGAGGATCTCTTCGAGCTTGTTGGCACCGGCGGCGAGAATGCCTTCCACAGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGATGGCGACCACGAACAGCATGGCGAAGAACATTTCGTTAACGGCGAGGAAGGCCTCCATCAGGCGGCGACCGGACTGGTAAACCCACCACAGCACCAGGTTTTAGGCGCTCATCAGGCAAAAAGGGAGGGAGAGaa                                                                                                       >  2:525007/4‑300 (MQ=255)
                                                                                                                                                                                           cGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAac                                                                                                      <  1:439771/300‑1 (MQ=255)
                                                                                                                                                                                            gTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGGTGGAGAGAacc                                                                                                     >  2:89004/1‑300 (MQ=255)
                                                                                                                                                                                            gTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAacc                                                                                                     >  1:731361/1‑300 (MQ=255)
                                                                                                                                                                                                agGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCg                                                                                                 <  2:499906/300‑1 (MQ=255)
                                                                                                                                                                                                    tctcTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAg                                                                                             <  1:810787/300‑1 (MQ=255)
                                                                                                                                                                                                      tctTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTTATCAGGACCAAGAGGATGGAGGGAACCACCGCCAgcg                                                                                           >  2:468691/1‑300 (MQ=255)
                                                                                                                                                                                                                 gttgGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAg                                                                                <  1:875414/300‑1 (MQ=255)
                                                                                                                                                                                                                   tgGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGTGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTGCGTTAATGGCGCGGCAGGCGTCCATCAGTCGGCGAACGGGCTGGTAAGCCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCTAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGAc                                                                              <  1:876430/300‑1 (MQ=255)
                                                                                                                                                                                                                   tgGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGAc                                                                              <  1:467817/300‑1 (MQ=255)
                                                                                                                                                                                                                                           ccACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCAt                                                      >  1:502954/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                           tCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCNTCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAAt                      >  2:770015/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                     gcTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGAGGgaa            >  2:316910/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                               cACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCACAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGAAGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTaggcg  >  2:26780/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    
ACCGTGACGATGATAAGCACCATCAGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCG  >  NC_000913/4322492‑4323078

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