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,7860GA96.7% 175.2 / NA 61intergenic (‑364/+550)phnF/phnDputative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/28);  minor base C (1/1);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGAC  >  NC_000913/4322541‑4323020
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
aaCCGCGAACTGCTCCCCACAGGGAGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGatggtg                                                                                                                                                                                                                                          <  1:63148/251‑1 (MQ=255)
     cgatccggctccctctggtgggcGCCGTTCCCGTCTGCGCCGACTAGGCCGACGACGGTCGCCGAGCGGCCGTTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCTGCAGTGGCATCACCTGTGGCGGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCCCCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGAt                                                                                                                                                                                                                                     <  2:397085/229‑1 (MQ=255)
        aTCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATaa                                                                                                                                                                                                                                   >  1:197077/1‑250 (MQ=255)
        aTCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACATGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCCCCGGGCCGGGCTCAATCGCTTCCCCCGCTTCGGAAAGCAGCATGGGGAGCACGCCGATGGGGTGGaaaa                                                                                                                                                                                                                                   >  2:496667/1‑247 (MQ=255)
         ccGCTTCCCTCAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAAc                                                                                                                                                                                                                                 <  1:277163/250‑1 (MQ=255)
          cGCTACCCACAGGGTAACGCCGGTCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGGTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGATCACGCCGATGGTGTGGATAAAc                                                                                                                                                                                                                                 <  1:47347/250‑1 (MQ=255)
            cTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGCTTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGCTGGTGTGGATCAACAgc                                                                                                                                                                                                                              >  2:339843/1‑251 (MQ=255)
               cccACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgcc                                                                                                                                                                                                                           >  1:228119/1‑251 (MQ=255)
               ccaccgggggagcgcgtaccacGCCTGCGCCGACCAGGCCGACGACGGGCGCCGCGGGGACGGTGGATACGAAGCGATAGAGGGAGTAGGAGATCGGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGAGCGCTTCGAGCTTGTGGGCCCCGGTGGCGCGAATGCCTCCCACCGGGCCGGGCTCAATCGCTTCCACCGCTACGGAAAGCCGCTTGGAGAGCACGCCGATGGTTTGGATAACCAGcgcc                                                                                                                                                                                                                           <  2:39964/230‑1 (MQ=255)
               ccaaaggggtgcggcggtgcccGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGCCGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgcc                                                                                                                                                                                                                           <  2:176424/232‑1 (MQ=255)
                    aGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGCGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGCGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGCGTGGATAAACAGcgccagcac                                                                                                                                                                                                                      >  2:111768/1‑251 (MQ=255)
                        ggccgccaactccGCCTGCGCCGACCCTGCCGACGACTGTCGCCGAGCGTACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCAGGAGTGGCATCACCGGTGGCAGCACGCCGTAGAGGAGCTCGTCGATCTTGGTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGAGCTCAATCGCTGCCCCCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGcgccagcacgcca                                                                                                                                                                                                                  <  2:104412/240‑1 (MQ=255)
                               gATCCCGCCTGCGCCGAACATGCCGACGACGGTCGCAGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGAAGCACGCAGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGACTGCCTTACACCGGGCAGGGCTCAATCGCTTAAACCGCTTCGGCAAGCAGCTTGGCGAGCACGCCGATGGTGTGGATAAACAGCGCCagaccgcccgcacgcc                                                                                                                                                                                                           >  1:496988/1‑251 (MQ=255)
                                 tCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCagcacgccagcacgcc                                                                                                                                                                                                           >  1:428575/1‑249 (MQ=255)
                                        tGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACgg                                                                                                                                                                                                  >  1:369592/1‑251 (MQ=255)
                                               aCCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGggccgagg                                                                                                                                                                                           >  1:485389/1‑251 (MQ=255)
                                                    gCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCAGAGGCCGAc                                                                                                                                                                                      >  1:232182/1‑251 (MQ=255)
                                                             ggTCGCCGAGCGGGCGTTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCAGCAGGGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCAc                                                                                                                                                                             <  2:68200/251‑1 (MQ=255)
                                                                 cccGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCTCCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGAGAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGTCGACGGCGACCACGAAc                                                                                                                                                                         <  1:173350/250‑1 (MQ=255)
                                                                   cGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAAc                                                                                                                                                                                                    >  2:373607/1‑222 (MQ=255)
                                                                   cGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAAc                                                                                                                                                                                                    <  1:373607/222‑1 (MQ=255)
                                                                    gAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGc                                                                                                                                                                      <  1:371015/251‑1 (MQ=255)
                                                                       gggACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGAGCAGCAGTGGCATCACCTGTGGCAGCACGCCGGAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATg                                                                                                                                                                   <  2:239528/250‑1 (MQ=255)
                                                                        ggACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCCCGAACAGCATgg                                                                                                                                                                  >  2:93113/1‑251 (MQ=255)
                                                                        ggACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGGGGATAAACAGCGCCAGAAAGCCAGCAACCCAGCGACCGGGCCGGGGCCGACGGGCACAACGAACAACAttg                                                                                                                                                                  >  2:346945/1‑249 (MQ=255)
                                                                                    cGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCCCGCCAGCGAACGGGCCGCGGCCGACGGCGACCACGCACAGCATCGCGAAGACCAttt                                                                                                                                                      >  2:157090/1‑251 (MQ=255)
                                                                                    cGAAGCGAGAGAGGGAGTAGGAGATCAGCAGTGGCATCCCCTGGGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCAttt                                                                                                                                                      <  2:441288/251‑1 (MQ=255)
                                                                                    cGAACCGATAGAGGGAGTGGGAGCTCAGTAGGGGCAACACTTGTGGCAGCCTGCCTTAGATGATTGCTCCGAGGTTGTTGGCACCGGTGGCGAGAATGCTTTCCGCCGGGCAGGGCTCAATCGCTTCACCCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTTAGGATAAACAGCGCAAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCAtt                                                                                                                                                       <  2:114136/250‑1 (MQ=255)
                                                                                                 gtaGTAGGAGATCAGCAGTGGCATCCCCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCCCCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAACCGCTTCCACCGTTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGcgcg                                                                                                                                         <  2:289027/249‑1 (MQ=255)
                                                                                                        gagaTCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCg                                                                                                                                  <  1:198616/251‑1 (MQ=255)
                                                                                                          gatcagacgtgtcatgtccgggggCAGCACGCCGGAGAGGATCTCGTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCGCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGGGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGTCCGAGGCCGACGGGGACCTCGAACAGCATGGCGATGACCATTACGTTAATGGCGCGGGAGGCGTc                                                                                                                                <  2:59654/229‑1 (MQ=255)
                                                                                                                   tGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACCGCGACCCCGCACAGCATGCCGACGACCATTTCGTTAAGGGggcggccgcgcggcctacggcc                                                                                                                       >  2:177115/1‑235 (MQ=255)
                                                                                                                    ggCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAggcgg                                                                                                                      >  1:152393/1‑251 (MQ=255)
                                                                                                                     gcagcccttgtgGCAGCCCGGCGTATAGGCTCGCTCCGAGCTTGTTGGCCCCGGTGGCGCGAATGCCTTCTTCCGGGCCGGGTGCAATCGCTTCCCCCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGAATGTGTGGATAAACAGCGCCAGCACGGCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCAGGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCTTCCATCAggcggc                                                                                                                     <  2:314895/243‑1 (MQ=255)
                                                                                                                      cATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGgcggcg                                                                                                                    <  1:97079/251‑1 (MQ=255)
                                                                                                                                gtcagtaagacGTATAGGATCTTGTCGAGCTCGTTGGCTACGGGGGCGCGAATGCCGTCCACCGGGCCGGGTTCAATCGCTCCCCCCGATTCGGAAACCAGCGTGGAGAGCACGTCGTGGGGGTGGATAAACAGCGCCAGCAAGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTgg                                                                                                          <  2:333868/241‑1 (MQ=255)
                                                                                                                                ggCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAAAAGCGCCAGCACGCCCGCACGCCAGCGAACGGGCCGGGGCCGAAGGAGACCACGAACAGCATGGCGAAGACCAATTCGTTAAAGGCGCGGCCGGCGTCCCTCAggcggagggcgggccgg                                                                                                          >  2:321753/1‑239 (MQ=255)
                                                                                                                                     caCGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTGCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAAc                                                                                                     <  1:338909/251‑1 (MQ=255)
                                                                                                                                        gCCGTAGAGGATCTCGTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACcca                                                                                                  <  1:463395/251‑1 (MQ=255)
                                                                                                                                                  aTCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCCGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCCCAACGGCAcc                                                                                      >  2:89072/1‑251 (MQ=255)
                                                                                                                                                      ctTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCa                                                                                     <  1:97740/248‑1 (MQ=255)
                                                                                                                                                          gAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCCGGCGTCCATCAGGCGGCGAACCGGCTGGGAAACCCACCACGGCACCCGGGTTTc                                                                              >  2:106436/1‑251 (MQ=255)
                                                                                                                                                           aGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGGTTTCg                                                                             >  1:137639/1‑251 (MQ=255)
                                                                                                                                                                gttgGCTCCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCGCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCt                                                                        <  2:446455/251‑1 (MQ=255)
                                                                                                                                                                gttgGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGgcgcc                                                                        >  1:359466/1‑250 (MQ=255)
                                                                                                                                                                   ggcgccGGTGGCGCGAATGCCTTCCACCGGGCCGGGCGCAATCGCTTCCACCGCTTCGGAAAGCGGTTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCtcat                                                                    <  2:135733/247‑1 (MQ=255)
                                                                                                                                                                     cACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCtcatca                                                                  <  1:96674/251‑1 (MQ=255)
                                                                                                                                                                          gTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCATGcc                                                              >  1:193225/1‑250 (MQ=255)
                                                                                                                                                                                 gAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCCCGCCCGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCCTCCGGCGGCGAACGGGCTGGTAAACCCACCCATGGACCAGGTTTTCGGGGCTCATCCGGCCAAAggggg                                                      >  2:283003/1‑250 (MQ=255)
                                                                                                                                                                                 gAATGCCTTCCACCGGGCCGGACTCAATCGCTTCCACCGCTTCGGAACGCAGCTTGGAGCGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACAATTTCGGTCATGGCGCGGCAGGCGCACATCAGGCGGCCACCGGGCTGGTAAGCCCCCCACACGCCCAGGTTTTCCGCGCtcctccagctcaaggggg                                                       >  2:259295/1‑235 (MQ=255)
                                                                                                                                                                                    tGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGGACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCCTCAGGCCAAAGGGGATgg                                                   >  1:89692/1‑251 (MQ=255)
                                                                                                                                                                                     gCCTTCCACCGGGCCGGGGTCAAACGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATgg                                                   <  1:34864/250‑1 (MQ=255)
                                                                                                                                                                                      ccTTCCCCCGGGCCGGGCTCACACTCTTCCTCCGCTTCGGAAAGCGGCTTGGATAGCACGCCTCGGGTGGGGAGAAACAGCGCCAGCCCGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGGTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCCCGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGgag                                                 <  2:238023/251‑1 (MQ=255)
                                                                                                                                                                                       cTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGACGACCATTTCGTTAATGGCGCGGCCGGCGTCCATCCGGCGGCGAACGGGCTGGTAAACCCACCCCCGCACCCGGGTTTTGGCGCCCACCCGGCCAAAGGGGATGgagg                                                >  2:62010/1‑250 (MQ=255)
                                                                                                                                                                                             ccgggccgggCGCAATCGCTGCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAacca                                          <  2:443152/251‑1 (MQ=255)
                                                                                                                                                                                               gggccgggcTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAaccacc                                        >  1:413436/1‑251 (MQ=255)
                                                                                                                                                                                                   cgggcTCACTCGCTGCCCCCCGTTCGGAAAGCGGCTGGGAGAGCACGCCGAGGGGGGGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCAGGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCCCCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCa                                    <  2:152393/251‑1 (MQ=255)
                                                                                                                                                                                                    gggcTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCAAGCCAGCACGCAAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCAACAGGCGGGGGAAGGGCTGGAAAACCCACCCCTGCACCAGGTTTTCGGCGTCAATCCGGCCCAAGGGGCAGGTGAGAACCACCGcccg                                   >  2:116002/1‑249 (MQ=255)
                                                                                                                                                                                                          aaTCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGGAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCggg                             >  1:168564/1‑250 (MQ=255)
                                                                                                                                                                                                                   cACCGCTGCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGAc                    <  1:157090/251‑1 (MQ=255)
                                                                                                                                                                                                                        cTTCGGAAAGCAGCTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCCTCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCg                >  1:451674/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cTTGGAGAGCACGCCGATGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGAAGGAGTCCATCAAGCGGCGAACGGGCTGGGAAACCCAACACGGCACCAGGTTATCGGAGATCATAAGGCCAAAGGGGATGGAGAGAACCACCGCCAACGCGGAGACCCAGACGGCGATTTGCAGCGTGAc  >  1:365594/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
AACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGCGACCACGAACAGCATGGCGAAGACCATTTCGTTAATGGCGCGGCAGGCGTCCATCAGGCGGCGAACGGGCTGGTAAACCCACCACGGCACCAGGTTTTCGGCGCTCATCAGGCCAAAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGAC  >  NC_000913/4322541‑4323020

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