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,7860GA89.8% 125.2 / 1.1 49intergenic (‑364/+550)phnF/phnDputative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (20/24);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGT  >  NC_000913/4322539‑4323017
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
gaacccgctactgctgtCCCACGGGTTGCCGCAGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCgatgg                                                                                                                                                                                                                                      <  1:111991/235‑1 (MQ=255)
         aaTCGCTTCCCACAGGGCGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATa                                                                                                                                                                                                                              >  1:300745/1‑250 (MQ=255)
             cctacctcggtggtcgcctggtgcccGCCTGCGCCGACCCTGCCGCCATCGTCCTCCGAGCTGACGTGGGATTCGATGTGAGAGAGGGAGTAGTAGATCAGCAGTGGCATCCCCTTGGGCATCACGCCGTAGAGGATCCCTTCGAGCTTGTTGGAGACTGTGGCGCGAATGCCTTCCACCGGGACGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGTTTAAACAg                                                                                                                                                                                                                         <  2:297879/228‑1 (MQ=255)
                  ccccatggTGACGCCGAACCCGCCTGCGCCGTCCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgcca                                                                                                                                                                                                                    <  1:170178/245‑1 (MQ=255)
                  ccACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgcc                                                                                                                                                                                                                     >  1:377694/1‑250 (MQ=255)
                       gggTGACGCCGATCACGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgccagcacg                                                                                                                                                                                                               >  1:124573/1‑251 (MQ=255)
                                 gATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCagcacgccagcacgcc                                                                                                                                                                                                     <  1:130690/251‑1 (MQ=255)
                                 gATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCagcacgccagcacgcc                                                                                                                                                                                                     <  1:33162/251‑1 (MQ=255)
                                  aTCCCGGCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAgcactccagcacgcca                                                                                                                                                                                                    <  1:16702/251‑1 (MQ=255)
                                   cccccccTGCGCCGACCATGCCGACGACGGTGGCCGAGCGGACGTGGGATGCGAAGCGAGAGAGAGAGTAGGAGAGCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTGCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGcacgccagcacgccag                                                                                                                                                                                                   <  2:66615/246‑1 (MQ=255)
                                    ccagCCTGCGCCGACCATGCCGACGACGGTAGCAGAGCGGACGTTGGATTCGAAGCGATGGAGGGAGTAGGCGATCAGCAGTGGCATCACCTGTGGCAGAACGACGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTACCACCGGGCCGGGCTCAATCGCTTCCACCGCATCGGAAAGCAGCTTGGAGAGCACGCCGGTGGGGTGGATAACCAGCGCCACCacgccggcacgccagc                                                                                                                                                                                                  >  2:169424/4‑251 (MQ=255)
                                        ccTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTGCGAGCTTGTGGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAAc                                                                                                                                                                                              <  1:52444/251‑1 (MQ=255)
                                                    aTGCCTACGACGGTCGCCGAGTGGACGTTGGATTCGATGCGAAAGAGGGATTAGGATATCAGCAGTGGCATCACCTGGGGCAGCACGCCGTAGAGGATCTCTTCGGGCTTGTTTGCACCGGTTGCTCGAATGCCTTCCACCGGGCCGGGCTCAAGCGCTGCCACCGCTGCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGggccgaggccg                                                                                                                                                                                  <  1:385043/251‑1 (MQ=255)
                                                             acgGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGa                                                                                                                                                                           <  1:281715/249‑1 (MQ=255)
                                                                   gCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCAGCACCTGTGTCAGCACGCCGGAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCCCCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCGGATGGTGTGGAAAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGCGCCGAGGCCGACGGCGACCACGAAc                                                                                                                                                                   <  1:308204/251‑1 (MQ=255)
                                                                          ggactgtGGAGTCGATGCTATAGAGGGAGTAGGAGATCAGCAGTGGCATCTCCTGGGGCAGCACGCTGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCGTCCACCGCTGCGGAAAGCAGCTTGGAGATCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGTGCCGAGGCCGACGGCGACCACGAACAGCATgg                                                                                                                                                            <  2:179583/245‑1 (MQ=255)
                                                                             cGTTGGATTCGAAGCGATAGAGGGCGTAGGTGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGCATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGCAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCCCGCCAGCACGCCCGCTAACGGGCCGAGGCCGACGGCGACCACGaccacaatggcgc                                                                                                                                                         >  2:74933/1‑239 (MQ=255)
                                                                              gTTGGATTCGAAGCGATAGAGGGAGTCGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGAACGCCGATGGTGGGGATAAACAGCGCCAGCACGCCCGCCCGCCAGCGAACCGACCGAGGCCGACGACGACCACGCACAGCCTGGCGa                                                                                                                                                         >  1:337768/1‑250 (MQ=255)
                                                                                        aaGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCg                                                                                                                                              >  1:360428/1‑251 (MQ=255)
                                                                                          gCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGtt                                                                                                                                            <  1:371131/251‑1 (MQ=255)
                                                                                            gATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTaa                                                                                                                                          <  1:248961/251‑1 (MQ=255)
                                                                                                   ggAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACTATTTCGTTAATGGcgcg                                                                                                                                   >  1:284840/1‑251 (MQ=255)
                                                                                                      gTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCa                                                                                                                                <  1:248217/251‑1 (MQ=255)
                                                                                                      gTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCa                                                                                                                                <  1:185040/251‑1 (MQ=255)
                                                                                                      gTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCATTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTAGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCCGCGAACGGGCCGAGGACGGCGGCGACCACGAAAAGCATGGCGGAGAACCTTTCGTTAATGGCGCGGc                                                                                                                                 >  2:303200/1‑250 (MQ=255)
                                                                                                        aGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGCCGGCGACCACGAACAGCATg                                                                                                                                                             >  1:220806/1‑220 (MQ=255)
                                                                                                         ggAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGc                                                                                                                             >  1:5445/1‑251 (MQ=255)
                                                                                                          gagtacagcagGGGCCTCCCCTTTGGCAGATCGCCTTAGAGGATCTCTGCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGTTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCTGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCTCGGCAGGCg                                                                                                                            <  2:373228/246‑1 (MQ=255)
                                                                                                          gagaTCAGCAGTGGCATCACCTGTGGCAGCACGCCTGAGAGGATCTCTTCGAGCTTGTTGGCTCCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAATACCATTTCGTTAATGGCGCGGCAGGCg                                                                                                                            <  1:333015/251‑1 (MQ=255)
                                                                                                                    gTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAggg                                                                                                                  >  1:110753/1‑250 (MQ=255)
                                                                                                                    gTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAggc                                                                                                                  <  1:108015/251‑1 (MQ=255)
                                                                                                                     tGGCATCACCTGTGGGAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGACCATCAggcg                                                                                                                 >  1:366108/1‑251 (MQ=255)
                                                                                                                             ccTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGaa                                                                                                            <  1:314531/248‑1 (MQ=255)
                                                                                                                               tgtgGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCAACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCAACAGGCGGCGACCGGGc                                                                                                       >  1:165566/1‑251 (MQ=255)
                                                                                                                                  ggCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGCCAGCGAACGGGccgg                                                                                                    >  1:291098/1‑248 (MQ=255)
                                                                                                                                     agcaAGCCGTAGAGGATCTCTTCGATCTTGTGGTCCCCGGTGGCGCGAATGCCGTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTaa                                                                                                 <  1:242215/251‑1 (MQ=255)
                                                                                                                                             gTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACccaccc                                                                                         >  1:15399/1‑250 (MQ=255)
                                                                                                                                               agagGATCTCTTCGATCTTGTTGGCTCCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCCCCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCcaccac                                                                                        <  1:13436/250‑1 (MQ=255)
                                                                                                                                                        ctTCGAGCTTATTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCa                                                                                 <  1:279628/248‑1 (MQ=255)
                                                                                                                                                           cGAGCTTGTTGGCACCGGGGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCTTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGtttt                                                                           <  1:134591/251‑1 (MQ=255)
                                                                                                                                                                  gttgGCACCGGTGGCTCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTGGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCt                                                                    <  1:352926/251‑1 (MQ=255)
                                                                                                                                                                       cACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAACCAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCCCCAGGTTTTCGGCGCtcatc                                                                >  1:159931/1‑250 (MQ=255)
                                                                                                                                                                          cgtaggggcgattgcctcccgtcgggccgggcTCAATCTCTTCCCCCGCTTCGGAAAGCCGGTTGGAGAGTACGGCGAGGGTGTGGATAAACAGCGCCAGCACGCCTGCCCGCCAGCGAACGGGCTGAGGCCGACGGCGACCACGAACAGCATGGTGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCTGCGCTCATCAGGc                                                            <  2:381869/229‑1 (MQ=255)
                                                                                                                                                                          cGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGAAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGc                                                            >  1:337746/1‑251 (MQ=255)
                                                                                                                                                                          cGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCAACAGGCGGCGAACGGGCTGGTAAAACCACCAAGGCACCAGGTTTTCGGCGCTCATCAGGc                                                            >  1:273006/1‑251 (MQ=255)
                                                                                                                                                                                   gAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGGTTTCGGCGCTCATCAGGCCAAAggggg                                                   >  1:186871/1‑250 (MQ=255)
                                                                                                                                                                                             cACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAAACCACCACGGCACCAGGTTTTCGGCGCTCCTCAGGCCAAAGGGGATGGAGAGAac                                         >  1:290246/1‑251 (MQ=255)
                                                                                                                                                                                                 ggtccggtctcaactgcctcacccGCTTCGGAAAGCAGCGTGGAGGGAACGCCGATGGGGTGGATAAACAGCGCCTGCACGCCAGCACGTCAGCGAACGGGCCGAGGCCGACGGCTACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAaccacc                                     <  2:336350/229‑1 (MQ=255)
                                                                                                                                                                                                                                  ggcgcgtgggtcgccggcgcGGTGGTTGGATAACCTGCCCCGGCGCGCCAGCGCGCCAGCGAACGGGCGGGGGCGGGTGGCGACCAAGAAACGGATGGCGAAGACCATCTCGTTAATGGCGCGGAGGGCGTCCATCAGGCGGCGAACGTGCGGGTAAACCCCCCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGAGGGAGCGAACCTCCGCCAGCGCGGTGCCCCAGACGGTGATTTGCAGCg   <  2:83632/232‑1 (MQ=255)
                                                                                                                                                                                                                                    cagcTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCCTCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAAAGAACCACCGCCAGCGCGGTGCCCCCGCCGGCGATTTGCAGCGt  >  1:288363/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       
GAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGT  >  NC_000913/4322539‑4323017

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