Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,027 T→C S372S (AGT→AGC)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,0270TC98.6% 198.1 / ‑6.9 70S372S (AGT→AGC) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/38);  total (30/40)
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

GGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTT  >  NC_000913/4626782‑4627273
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
ggagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCg                                                                                                                                                                                                                                                      >  1:438021/4‑250 (MQ=255)
    aTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                                                                                                                                                                                                                                                  >  1:290367/1‑250 (MQ=255)
      gcaGTGACCTCCGCCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGTTGGAGATTCAGGCGTTGGTCGATCACTCGATGATGGCGACTCCGCGCCTCGTGGCAGTGGTGCTGGACCAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATTGCCGATCAGGATGTGTTTGTTAACG‑TGTTCGGCGGCTTGTAGGTAACCGAAACCAGCGCCGACTTAGc                                                                                                                                                                                                                                               <  2:348859/249‑1 (MQ=255)
       aaGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCg                                                                                                                                                                                                                                              >  1:79731/1‑251 (MQ=255)
            cccTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑Act                                                                                                                                                                                                                                        <  2:550322/250‑1 (MQ=255)
             ccTCCGGCAGCTCCGTGATGGTGGTGTGGGCAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGGCTGCAAATGGCCGTTCAGGATGTGTTTGTGACCG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCCGCGCCGACTTAGCGTT‑Actg                                                                                                                                                                                                                                       >  2:313528/1‑251 (MQ=255)
                ccggcagctccggagtgggggtgtGGGAAGGATCGTGTCCGGTGCTGGGGGAGCATCAGGCCCTGGTCGAGCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTTGGGCTGGAACAAAACCGTCTG‑GCAATCCTGAGGGCGGTCTTGCACCGGCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAATCCTGCGCCGACTTAGCGTT‑Tctgct                                                                                                                                                                                                                                     <  2:527580/227‑1 (MQ=255)
                 cggcaGCTCCGTGATGGTGGTGTGGGAAGTAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAATCCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAATCCAGCGCCGACTTATCGGT‑Actgcttg                                                                                                                                                                                                                                   <  1:131991/247‑3 (MQ=255)
                    caGCTCAGTGATGGTGGGGTGGGAAGGAAAGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGGACAAAACAGGCTG‑GCAATACTGCTGGCTGTCTTGCACCGGCACGGTGGTCTGCAACTGGCCGATCAGGATGTGTTTGTGAACG‑TAGTCGTCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTA‑CGTTTACTGCGGGCg                                                                                                                                                                                                                                 >  2:516594/2‑250 (MQ=255)
                       ctccgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCAGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCGGGGTGTGTTTGTGAACGGTGG‑CGGCGGCGGGAAGGTAACCGAAAACCGCTACGACTTAGCGTT‑ACTGCTGGGGATgg                                                                                                                                                                                                                             >  2:454887/5‑251 (MQ=255)
                             gtttcgtgagtGGGATGGAACGCGTTCGTTGCTTGTGGCGAGTCAGGCGCTGGTCGGTCCCGCGATTAGGGCTCATCCGCGCGGCGGGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCGGCTGGCTGGCTTGCACCGTCACGGTGGGCTGCAAATGGCCGATCAGGAGGTGTCTGTGATCG‑TGGTCGGCGGCGTGAAGTTAACCGAAACCAGCGCAGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCa                                                                                                                                                                                                                       <  2:71762/241‑1 (MQ=255)
                             gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCa                                                                                                                                                                                                                       <  1:641681/251‑1 (MQ=255)
                                            aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAAAAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGACGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGGT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                                                                                                                        >  2:523349/1‑251 (MQ=255)
                                             ggaaCGCGTCCGTTGCTGGTGGAGATTCCAGCGCTGCTCCCTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAACCCGTATG‑GCAATCACGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGTTGTGTCTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAGCCAGCGCCGGCTTAGCGTT‑ACTGCTGGCGATGGTTTACCGCCTGCGTGGCAGAcc                                                                                                                                                                                                       >  2:660129/1‑251 (MQ=255)
                                              gaacgctctctttttgggtggAGATTCAGGTGCTGGTCGATCACTGGATGATGGCGAATCCGCGCCGCGTGGCAGGGGGGCTGGAACAAACCCGTCTG‑GCACTCCTGCTGGCTGTCTTGCACCTTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑GGGTCGGCGGCGTGAAGGTAACCGAAACCAGCTCCGACTAAGCGTT‑TCTGCTGGCGATGGGTTCCAGCCTGCGTGACAGACCg                                                                                                                                                                                                      <  2:410471/235‑1 (MQ=255)
                                                 cggttccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTG‑GCAGTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                                                                                                                   <  2:648762/242‑1 (MQ=255)
                                                    gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGACTGCGGGACAGACCGCTGccc                                                                                                                                                                                                >  2:236611/7‑250 (MQ=255)
                                                    gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                                                                                                                 >  1:545516/7‑250 (MQ=255)
                                                       cgttgctgGTGGAGATTCAGGCGCTGTTCGATCACGCGATGATGGCGAATCCGCGCCTCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAg                                                                                                                                                                                             <  1:97439/248‑1 (MQ=255)
                                                          tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                           <  1:50588/250‑1 (MQ=255)
                                                          tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                           <  1:594237/250‑1 (MQ=255)
                                                               gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                                                                                                                     <  1:660129/251‑1 (MQ=255)
                                                                       tCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGttt                                                                                                                                                                              <  1:584484/250‑1 (MQ=255)
                                                                                     aTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                                                                                                                               <  1:313626/251‑1 (MQ=255)
                                                                                        ccTCGATGATGGGGAATCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCa                                                                                                                                                            <  2:438021/250‑1 (MQ=255)
                                                                                            gatgatGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTACGTGACAGGCCGCTGCCACCGGATATGGTGGTGTTTGGtagaggagggctggcggggg                                                                                                                                                        >  2:675795/1‑232 (MQ=255)
                                                                                             atgatgGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGCACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGGAACCGAAACCAGCGCCGACTTAGCGTT‑CCTGCTGGCGATGGTTTCCAGCCTGCGTGGCAGACCGCTGCCACAGGTTCTGGTGGGGTTTTGTGAAGTCGGGcgggacggggg                                                                                                                                                       >  2:31547/1‑241 (MQ=255)
                                                                                              tgatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                                                                                                                        >  1:428578/1‑249 (MQ=255)
                                                                                                 gggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATc                                                                                                                                                   <  1:323335/250‑1 (MQ=255)
                                                                                                      aatccgcgccgcGTGGCAGTGGAGATGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGCACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCAACAGGTTCTGGTGGTGTTTGGGGAGGTCGGGCTGGCGGGGGAGATCCGccc                                                                                                                                              >  2:465289/7‑251 (MQ=255)
                                                                                                         ccgcgccgcGTGGCAGTGGGGCTGGGACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCTCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGt                                                                                                                                           <  1:606953/248‑1 (MQ=255)
                                                                                                                          ggggCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCGCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                                                                                                                          <  1:646331/251‑1 (MQ=255)
                                                                                                                              cTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                      >  1:83922/1‑251 (MQ=255)
                                                                                                                                ggAACAAAACCG‑GTGTGCAATCCGGGTGGCTGTTTTGCTCCGTCACGGTGGTCTGCAAATGGCCGAGCAGGATGTGTTTGTGAACG‑TGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                      <  2:10935/249‑1 (MQ=255)
                                                                                                                                 gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTCGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGACGGGGGGATCCGCCCGGCGCCCAGCgggatggtcagcaatgc                                                                                                                   >  2:446311/1‑236 (MQ=255)
                                                                                                                                 gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                                                                                                                   <  1:651861/251‑4 (MQ=255)
                                                                                                                                     aaaaCCGTCGG‑GCAAACCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcagaa                                                                                                               <  2:428578/251‑8 (MQ=255)
                                                                                                                                                ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACG‑TGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTT‑ACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac                                                                                                     <  2:272108/244‑1 (MQ=255)
                                                                                                                                                   aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑g                                                                                                 <  2:386556/248‑1 (MQ=255)
                                                                                                                                                          cgggctttcttGCACCGTCTCGGTGTTCTGCAAATGGCCGATCTGGATGTGGTTGTTAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTGCCAGCCGGCGTGACATACCGCTGCCATAGGAGCTGGTGGTGTTTGGTGAAGTCTGGCTGGCAGGGGAGATCCGCCGGTTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCg                                                                                          <  1:384155/242‑1 (MQ=255)
                                                                                                                                                               tgtcttGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                    >  1:344473/5‑173 (MQ=255)
                                                                                                                                                               tgtcttGCACCGTCACGGTGGGCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                    <  2:344473/169‑1 (MQ=255)
                                                                                                                                                                gtctttcccgggagtgggggtCTGCAAATGGCCGACCTGGATTTGTTTGTGAACG‑TGGTCGGGGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCGGCGTTAGAGACCGCTGCCTCAGGATCTGGTGGTGTTTGGTGAAGCGGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCggg                                                                                     <  2:505363/232‑1 (MQ=255)
                                                                                                                                                                      cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGTGTGAAGGTAACCGAAACCAGCGCCGTCTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGt                                                                              <  1:273977/251‑1 (MQ=255)
                                                                                                                                                                        ccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTc                                                                            >  1:609796/1‑251 (MQ=255)
                                                                                                                                                                            cacggtgggctgcgagcggcggagcAGGATGTGGTGGTGATCG‑GGGTCGGTGGCGTGAGGGCAGCCGAAACCAGGGCCGACTTAGCGTT‑ATTGCTGGCGTTGGTTTCCAGCCTGCGTGAGAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCATGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGTCGGGCGATTGTTCcgg                                                                         <  2:265148/226‑1 (MQ=255)
                                                                                                                                                                            cACGGTAGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGAAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGCCAGGAACGAATTTCAGAAGCGGCGGAACAA‑GGATTTCGCCGGGCGTTTGTTCcggc                                                                        >  1:632105/1‑251 (MQ=255)
                                                                                                                                                                            aaCGGGGGTCTGAAAATGGCCTACCAGGATGTGTTGGGGAAAG‑TGGTCGGTGGCGTGAAGGTAACCGAATCCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCcgg                                                                         <  2:107367/249‑1 (MQ=255)
                                                                                                                                                                              cGGTGGTCGGCAAATGACCGATCAGGATGTGTTTGGGAACG‑TGGTCGGCGGCGTGAAGGAAACCGACACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGAGCACAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAA‑ACAGGTTTTCGACGGGAGCTTGTTCcggcgg                                                                      >  2:525456/1‑251 (MQ=255)
                                                                                                                                                                                     cTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                  >  1:64893/1‑153 (MQ=255)
                                                                                                                                                                                     cTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCGGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                  <  2:64893/153‑1 (MQ=255)
                                                                                                                                                                                       gaaaggggCCGATCGGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTT‑ACTGCTGGCGATGGTTGCCAGCCGGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGtac                                                             <  2:81565/245‑3 (MQ=255)
                                                                                                                                                                                             ggCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                          <  2:395950/137‑1 (MQ=255)
                                                                                                                                                                                             ggCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                          >  1:395950/1‑137 (MQ=255)
                                                                                                                                                                                              gCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTACAGTACCGaaaa                                                       >  2:491896/1‑250 (MQ=255)
                                                                                                                                                                                                    cAGGATGTGTTGTTGAAAG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGTCTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                 <  1:528094/250‑1 (MQ=255)
                                                                                                                                                                                                          gtgtTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAgg                                          <  1:236611/251‑1 (MQ=255)
                                                                                                                                                                                                               tgtgAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                      <  1:491896/250‑1 (MQ=255)
                                                                                                                                                                                                                gtgAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                        >  2:652604/1‑120 (MQ=255)
                                                                                                                                                                                                                gtgAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                        <  1:652604/120‑1 (MQ=255)
                                                                                                                                                                                                                gtgAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGATGc                                     >  2:232563/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCTGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGATCGAAttgcagaagcgg                                                                                                           <  2:463124/163‑12 (MQ=255)
                                                                                                                                                                                                                                   gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcgg                                                                                                           >  1:463124/1‑163 (MQ=255)
                                                                                                                                                                                                                                     gAAGGGAACCGAAACCAGCCCCGCCTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGTCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAc                <  2:83922/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAct               >  1:149109/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGGTAACCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCt                                                                    >  1:315603/1‑197 (MQ=255)
                                                                                                                                                                                                                                           aaCCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc           <  1:285628/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaCCGAAACCAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTAACGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTCAAAAACTCTCCg          >  1:73171/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt  >  1:576170/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccAGCGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt  <  1:523349/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
GGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACG‑TGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTT‑ACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTT  >  NC_000913/4626782‑4627273

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