Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,069 C→T R386R (CGC→CGT)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,0690CT95.7% 207.7 / NA 70R386R (CGC→CGT) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/39);  minor base G (1/2);  total (29/41)
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

AGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGA  >  NC_000913/4626826‑4627307
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAAAAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGACGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGGTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                                                                                                                                                          >  2:523349/1‑251 (MQ=255)
 ggaaCGCGTCCGTTGCTGGTGGAGATTCCAGCGCTGCTCCCTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAACCCGTATG‑GCAATCACGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGTTGTGTCTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAGCCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTACCGCCTGCGTGGCAGAcc                                                                                                                                                                                                                                         >  2:660129/1‑251 (MQ=255)
  gaacgctctctttttgggtggAGATTCAGGTGCTGGTCGATCACTGGATGATGGCGAATCCGCGCCGCGTGGCAGGGGGGCTGGAACAAACCCGTCTG‑GCACTCCTGCTGGCTGTCTTGCACCTTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGGGGTCGGCGGCGTGAAGGTAACCGAAACCAGCTCCGACTAAGCGTTTCTGCTGGCGATGGGTTCCAGCCTGCGTGACAGACCg                                                                                                                                                                                                                                        <  2:410471/235‑1 (MQ=255)
     cggttccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTG‑GCAGTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                                                                                                                                                     <  2:648762/242‑1 (MQ=255)
        gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGACTGCGGGACAGACCGCTGccc                                                                                                                                                                                                                                  >  2:236611/7‑250 (MQ=255)
        gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                                                                                                                                                   >  1:545516/7‑250 (MQ=255)
           cgttgctgGTGGAGATTCAGGCGCTGTTCGATCACGCGATGATGGCGAATCCGCGCCTCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAg                                                                                                                                                                                                                               <  1:97439/248‑1 (MQ=255)
              tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                                                             <  1:594237/250‑1 (MQ=255)
              tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                                                             <  1:50588/250‑1 (MQ=255)
                   gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                                                                                                                                                       <  1:660129/251‑1 (MQ=255)
                           tCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGttt                                                                                                                                                                                                                <  1:584484/250‑1 (MQ=255)
                                         aTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                                                                                                                                                                 <  1:313626/251‑1 (MQ=255)
                                            ccTCGATGATGGGGAATCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCa                                                                                                                                                                                              <  2:438021/250‑1 (MQ=255)
                                                gatgatGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTACGTGACAGGCCGCTGCCACCGGATATGGTGGTGTTTGGtagaggagggctggcggggg                                                                                                                                                                                          >  2:675795/1‑232 (MQ=255)
                                                 atgatgGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGCACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGGAACCGAAACCAGCGCCGACTTAGCGTTCCTGCTGGCGATGGTTTCCAGCCTGCGTGGCAGACCGCTGCCACAGGTTCTGGTGGGGTTTTGTGAAGTCGGGcgggacggggg                                                                                                                                                                                         >  2:31547/1‑241 (MQ=255)
                                                  tgatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                                                                                                                                                          >  1:428578/1‑249 (MQ=255)
                                                     gggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATc                                                                                                                                                                                     <  1:323335/250‑1 (MQ=255)
                                                          aatccgcgccgcGTGGCAGTGGAGATGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGCACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCAACAGGTTCTGGTGGTGTTTGGGGAGGTCGGGCTGGCGGGGGAGATCCGccc                                                                                                                                                                                >  2:465289/7‑251 (MQ=255)
                                                             ccgcgccgcGTGGCAGTGGGGCTGGGACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCTCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGt                                                                                                                                                                             <  1:606953/248‑1 (MQ=255)
                                                                              ggggCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCGCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                                                                                                                                                            <  1:646331/251‑1 (MQ=255)
                                                                                  cTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                                                        >  1:83922/1‑251 (MQ=255)
                                                                                    ggAACAAAACCG‑GTGTGCAATCCGGGTGGCTGTTTTGCTCCGTCACGGTGGTCTGCAAATGGCCGAGCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                                                        <  2:10935/249‑1 (MQ=255)
                                                                                     gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTCGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGACGGGGGGATCCGCCCGGCGCCCAGCgggatggtcagcaatgc                                                                                                                                                     >  2:446311/1‑236 (MQ=255)
                                                                                     gAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                                                                                                                                                     <  1:651861/251‑4 (MQ=255)
                                                                                         aaaaCCGTCGG‑GCAAACCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcagaa                                                                                                                                                 <  2:428578/251‑8 (MQ=255)
                                                                                                    ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac                                                                                                                                       <  2:272108/244‑1 (MQ=255)
                                                                                                       aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑g                                                                                                                                   <  2:386556/248‑1 (MQ=255)
                                                                                                              cgggctttcttGCACCGTCTCGGTGTTCTGCAAATGGCCGATCTGGATGTGGTTGTTAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGTGACATACCGCTGCCATAGGAGCTGGTGGTGTTTGGTGAAGTCTGGCTGGCAGGGGAGATCCGCCGGTTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCg                                                                                                                            <  1:384155/242‑1 (MQ=255)
                                                                                                                   tgtcttGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                                                      >  1:344473/5‑173 (MQ=255)
                                                                                                                   tgtcttGCACCGTCACGGTGGGCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                                                      <  2:344473/169‑1 (MQ=255)
                                                                                                                    gtctttcccgggagtgggggtCTGCAAATGGCCGACCTGGATTTGTTTGTGAACGTGGTCGGGGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTTAGAGACCGCTGCCTCAGGATCTGGTGGTGTTTGGTGAAGCGGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCggg                                                                                                                       <  2:505363/232‑1 (MQ=255)
                                                                                                                          cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGTGTGAAGGTAACCGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGt                                                                                                                <  1:273977/251‑1 (MQ=255)
                                                                                                                            ccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTc                                                                                                              >  1:609796/1‑251 (MQ=255)
                                                                                                                                cacggtgggctgcgagcggcggagcAGGATGTGGTGGTGATCGGGGTCGGTGGCGTGAGGGCAGCCGAAACCAGGGCCGACTTAGCGTTATTGCTGGCGTTGGTTTCCAGCCTGCGTGAGAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCATGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGTCGGGCGATTGTTCcgg                                                                                                           <  2:265148/226‑1 (MQ=255)
                                                                                                                                cACGGTAGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGCCAGGAACGAATTTCAGAAGCGGCGGAACAA‑GGATTTCGCCGGGCGTTTGTTCcggc                                                                                                          >  1:632105/1‑251 (MQ=255)
                                                                                                                                aaCGGGGGTCTGAAAATGGCCTACCAGGATGTGTTGGGGAAAGTGGTCGGTGGCGTGAAGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCcgg                                                                                                           <  2:107367/249‑1 (MQ=255)
                                                                                                                                  cGGTGGTCGGCAAATGACCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGAAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGAGCACAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAA‑ACAGGTTTTCGACGGGAGCTTGTTCcggcgg                                                                                                        >  2:525456/1‑251 (MQ=255)
                                                                                                                                         cTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                                                    >  1:64893/1‑153 (MQ=255)
                                                                                                                                         cTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCGGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                                                    <  2:64893/153‑1 (MQ=255)
                                                                                                                                           gaaaggggCCGATCGGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGtac                                                                                               <  2:81565/245‑3 (MQ=255)
                                                                                                                                                 ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                                                            <  2:395950/137‑1 (MQ=255)
                                                                                                                                                 ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                                                            >  1:395950/1‑137 (MQ=255)
                                                                                                                                                  gCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTACAGTACCGaaaa                                                                                         >  2:491896/1‑250 (MQ=255)
                                                                                                                                                        cAGGATGTGTTGTTGAAAGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                                                   <  1:528094/250‑1 (MQ=255)
                                                                                                                                                              gtgtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAgg                                                                            <  1:236611/251‑1 (MQ=255)
                                                                                                                                                                   tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                                                        <  1:491896/250‑1 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                                                          <  1:652604/120‑1 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                                                          >  2:652604/1‑120 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGATGc                                                                       >  2:232563/1‑250 (MQ=255)
                                                                                                                                                                                      gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCTGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGATCGAAttgcagaagcgg                                                                                                                                             <  2:463124/163‑12 (MQ=255)
                                                                                                                                                                                      gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcgg                                                                                                                                             >  1:463124/1‑163 (MQ=255)
                                                                                                                                                                                        gAAGGGAACCGAAACCAGCCCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGTCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAc                                                  <  2:83922/251‑1 (MQ=255)
                                                                                                                                                                                         aaGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAct                                                 >  1:149109/1‑251 (MQ=255)
                                                                                                                                                                                          aGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCt                                                                                                      >  1:315603/1‑197 (MQ=255)
                                                                                                                                                                                              aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc                                             <  1:285628/250‑1 (MQ=255)
                                                                                                                                                                                              aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTAACGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTCAAAAACTCTCCg                                            >  1:73171/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                    >  1:576170/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                    <  1:523349/251‑1 (MQ=255)
                                                                                                                                                                                                       cagcgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTa                                   <  1:446095/247‑1 (MQ=255)
                                                                                                                                                                                                        agcgCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGGTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGAAGCATTTTGGCGTTAAAAAGctcctcaccccgctt                                    >  2:163077/4‑237 (MQ=255)
                                                                                                                                                                                                              gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGGGtt                            >  1:166709/1‑251 (MQ=255)
                                                                                                                                                                                                              gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAA‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAGAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAAGCGCTTAGCGTGtt                            >  1:521929/1‑251 (MQ=255)
                                                                                                                                                                                                               agttAGCGTTACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGTGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCGCAGCGGTCAGGGACGAAGTGCAGAAGCGTCGAAACAC‑GGTTTTCGCCGGGCGATGGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTGTGGCGTTAAAAAACTCTCCGACGCGGGTAGTGTGTTc                           <  2:79731/249‑1 (MQ=255)
                                                                                                                                                                                                                    ggtgtACTGCTGGCGGGGGTTTCCAGCCTGCGTGACAGATCGCTGCCACAGTATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCTCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTGTGCGTGTTcgacg                       <  2:290367/246‑1 (MQ=255)
                                                                                                                                                                                                                       ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACtt                   >  1:303629/1‑251 (MQ=255)
                                                                                                                                                                                                                            gctgggtttGGTTTCCAGCCTTCGTGAAAGATCTCGGCCACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACAC‑GGGTTTCGCCTGGCGATTGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGTCGCGCTTAGCGTGTTCGACGATTTATAAt              <  2:539755/243‑1 (MQ=255)
                                                                                                                                                                                                                            gctgGCGATGGTGTCCAGCCGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATaa               <  1:609692/250‑1 (MQ=255)
                                                                                                                                                                                                                             ctgGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATg             >  1:636584/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCGATGGTTTCAAGCCTGCTGGACAGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTAGAATgagg          <  2:609796/251‑2 (MQ=255)
                                                                                                                                                                                                                                         tCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgga  <  1:525456/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
AGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTG‑GCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGA  >  NC_000913/4626826‑4627307

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