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→CGTradADNA 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

AGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGA  >  NC_000913/4626826‑4627307
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAAAAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGACGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGGTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                                                                                                                                                          >  2:523349/1‑251 (MQ=255)
 ggaaCGCGTCCGTTGCTGGTGGAGATTCCAGCGCTGCTCCCTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAACCCGTATGGCAATCACGCTGGCTGTCTTGCACCGTCACGGAGGTCTGCAAATGGCCGATCAGGTTGTGTCTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAGCCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTACCGCCTGCGTGGCAGAcc                                                                                                                                                                                                                                         >  2:660129/1‑251 (MQ=255)
  gaacgctctctttttgggtggAGATTCAGGTGCTGGTCGATCACTGGATGATGGCGAATCCGCGCCGCGTGGCAGGGGGGCTGGAACAAACCCGTCTGGCACTCCTGCTGGCTGTCTTGCACCTTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGGGGTCGGCGGCGTGAAGGTAACCGAAACCAGCTCCGACTAAGCGTTTCTGCTGGCGATGGGTTCCAGCCTGCGTGACAGACCg                                                                                                                                                                                                                                        <  2:410471/235‑1 (MQ=255)
     cggttccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTGGCAGTCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTg                                                                                                                                                                                                                                     <  2:648762/242‑1 (MQ=255)
        gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGACTGCGGGACAGACCGCTGccc                                                                                                                                                                                                                                  >  2:236611/7‑250 (MQ=255)
        gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                                                                                                                                                   >  1:545516/7‑250 (MQ=255)
           cgttgctgGTGGAGATTCAGGCGCTGTTCGATCACGCGATGATGGCGAATCCGCGCCTCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAg                                                                                                                                                                                                                               <  1:97439/248‑1 (MQ=255)
              tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                                                             <  1:594237/250‑1 (MQ=255)
              tgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGa                                                                                                                                                                                                                             <  1:50588/250‑1 (MQ=255)
                   gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGTTGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                                                                                                                                                       <  1:660129/251‑1 (MQ=255)
                           tCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGttt                                                                                                                                                                                                                <  1:584484/250‑1 (MQ=255)
                                         aTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                                                                                                                                                                 <  1:313626/251‑1 (MQ=255)
                                            ccTCGATGATGGGGAATCCGCGCCGCGTGGCAGTGGGGCGGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCa                                                                                                                                                                                              <  2:438021/250‑1 (MQ=255)
                                                gatgatGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGACCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTACGTGACAGGCCGCTGCCACCGGATATGGTGGTGTTTGGtagaggagggctggcggggg                                                                                                                                                                                          >  2:675795/1‑232 (MQ=255)
                                                 atgatgGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGCACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGGAACCGAAACCAGCGCCGACTTAGCGTTCCTGCTGGCGATGGTTTCCAGCCTGCGTGGCAGACCGCTGCCACAGGTTCTGGTGGGGTTTTGTGAAGTCGGGcgggacggggg                                                                                                                                                                                         >  2:31547/1‑241 (MQ=255)
                                                  tgatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                                                                                                                                                          >  1:428578/1‑249 (MQ=255)
                                                     gggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATc                                                                                                                                                                                     <  1:323335/250‑1 (MQ=255)
                                                          aatccgcgccgcGTGGCAGTGGAGATGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGCACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCAACAGGTTCTGGTGGTGTTTGGGGAGGTCGGGCTGGCGGGGGAGATCCGccc                                                                                                                                                                                >  2:465289/7‑251 (MQ=255)
                                                             ccgcgccgcGTGGCAGTGGGGCTGGGACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCTCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGt                                                                                                                                                                             <  1:606953/248‑1 (MQ=255)
                                                                              ggggCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCGCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAAc                                                                                                                                                            <  1:646331/251‑1 (MQ=255)
                                                                                  cTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                                                        >  1:83922/1‑251 (MQ=255)
                                                                                    ggAACAAAACCGGTGTGCAATCCGGGTGGCTGTTTTGCTCCGTCACGGTGGTCTGCAAATGGCCGAGCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAt                                                                                                                                                        <  2:10935/249‑1 (MQ=255)
                                                                                     gAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTCGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGACGGGGGGATCCGCCCGGCGCCCAGCgggatggtcagcaatgc                                                                                                                                                     >  2:446311/1‑236 (MQ=255)
                                                                                     gAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttc                                                                                                                                                     <  1:651861/251‑4 (MQ=255)
                                                                                         aaaaCCGTCGGGCAAACCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcagaa                                                                                                                                                 <  2:428578/251‑8 (MQ=255)
                                                                                                    ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac                                                                                                                                       <  2:272108/244‑1 (MQ=255)
                                                                                                       aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACg                                                                                                                                   <  2:386556/248‑1 (MQ=255)
                                                                                                              cgggctttcttGCACCGTCTCGGTGTTCTGCAAATGGCCGATCTGGATGTGGTTGTTAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGTGACATACCGCTGCCATAGGAGCTGGTGGTGTTTGGTGAAGTCTGGCTGGCAGGGGAGATCCGCCGGTTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCg                                                                                                                            <  1:384155/242‑1 (MQ=255)
                                                                                                                   tgtcttGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                                                      >  1:344473/5‑173 (MQ=255)
                                                                                                                   tgtcttGCACCGTCACGGTGGGCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCg                                                                                                                                                                                                      <  2:344473/169‑1 (MQ=255)
                                                                                                                    gtctttcccgggagtgggggtCTGCAAATGGCCGACCTGGATTTGTTTGTGAACGTGGTCGGGGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTTAGAGACCGCTGCCTCAGGATCTGGTGGTGTTTGGTGAAGCGGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACACGGTTTTCGCCggg                                                                                                                       <  2:505363/232‑1 (MQ=255)
                                                                                                                          cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGTGTGAAGGTAACCGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGt                                                                                                                <  1:273977/251‑1 (MQ=255)
                                                                                                                            ccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTc                                                                                                              >  1:609796/1‑251 (MQ=255)
                                                                                                                                cacggtgggctgcgagcggcggagcAGGATGTGGTGGTGATCGGGGTCGGTGGCGTGAGGGCAGCCGAAACCAGGGCCGACTTAGCGTTATTGCTGGCGTTGGTTTCCAGCCTGCGTGAGAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCATGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGTCGGGCGATTGTTCcgg                                                                                                           <  2:265148/226‑1 (MQ=255)
                                                                                                                                cACGGTAGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGCCAGGAACGAATTTCAGAAGCGGCGGAACAAGGATTTCGCCGGGCGTTTGTTCcggc                                                                                                          >  1:632105/1‑251 (MQ=255)
                                                                                                                                aaCGGGGGTCTGAAAATGGCCTACCAGGATGTGTTGGGGAAAGTGGTCGGTGGCGTGAAGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCcgg                                                                                                           <  2:107367/249‑1 (MQ=255)
                                                                                                                                  cGGTGGTCGGCAAATGACCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGAAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGAGCACAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAAACAGGTTTTCGACGGGAGCTTGTTCcggcgg                                                                                                        >  2:525456/1‑251 (MQ=255)
                                                                                                                                         cTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                                                    >  1:64893/1‑153 (MQ=255)
                                                                                                                                         cTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCGGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCggg                                                                                                                                                                                                    <  2:64893/153‑1 (MQ=255)
                                                                                                                                           gaaaggggCCGATCGGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGtac                                                                                               <  2:81565/245‑3 (MQ=255)
                                                                                                                                                 ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                                                            <  2:395950/137‑1 (MQ=255)
                                                                                                                                                 ggCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTg                                                                                                                                                                                                            >  1:395950/1‑137 (MQ=255)
                                                                                                                                                  gCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTACAGTACCGaaaa                                                                                         >  2:491896/1‑250 (MQ=255)
                                                                                                                                                        cAGGATGTGTTGTTGAAAGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                                                   <  1:528094/250‑1 (MQ=255)
                                                                                                                                                              gtgtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAgg                                                                            <  1:236611/251‑1 (MQ=255)
                                                                                                                                                                   tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                                                        <  1:491896/250‑1 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                                                          <  1:652604/120‑1 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGaa                                                                                                                                                                                                          >  2:652604/1‑120 (MQ=255)
                                                                                                                                                                    gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGATGc                                                                       >  2:232563/1‑250 (MQ=255)
                                                                                                                                                                                      gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCTGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGATCGAAttgcagaagcgg                                                                                                                                             <  2:463124/163‑12 (MQ=255)
                                                                                                                                                                                      gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcgg                                                                                                                                             >  1:463124/1‑163 (MQ=255)
                                                                                                                                                                                        gAAGGGAACCGAAACCAGCCCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGTCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAc                                                  <  2:83922/251‑1 (MQ=255)
                                                                                                                                                                                         aaGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAct                                                 >  1:149109/1‑251 (MQ=255)
                                                                                                                                                                                          aGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCt                                                                                                      >  1:315603/1‑197 (MQ=255)
                                                                                                                                                                                              aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc                                             <  1:285628/250‑1 (MQ=255)
                                                                                                                                                                                              aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTAACGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTCAAAAACTCTCCg                                            >  1:73171/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                    >  1:576170/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                    <  1:523349/251‑1 (MQ=255)
                                                                                                                                                                                                       cagcgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTa                                   <  1:446095/247‑1 (MQ=255)
                                                                                                                                                                                                        agcgCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGGTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGAAGCATTTTGGCGTTAAAAAGctcctcaccccgctt                                    >  2:163077/4‑237 (MQ=255)
                                                                                                                                                                                                              gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGGGtt                            >  1:166709/1‑251 (MQ=255)
                                                                                                                                                                                                              gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAAGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAGAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAAGCGCTTAGCGTGtt                            >  1:521929/1‑251 (MQ=255)
                                                                                                                                                                                                               agttAGCGTTACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGTGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCGCAGCGGTCAGGGACGAAGTGCAGAAGCGTCGAAACACGGTTTTCGCCGGGCGATGGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTGTGGCGTTAAAAAACTCTCCGACGCGGGTAGTGTGTTc                           <  2:79731/249‑1 (MQ=255)
                                                                                                                                                                                                                    ggtgtACTGCTGGCGGGGGTTTCCAGCCTGCGTGACAGATCGCTGCCACAGTATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCTCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTGTGCGTGTTcgacg                       <  2:290367/246‑1 (MQ=255)
                                                                                                                                                                                                                       ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACtt                   >  1:303629/1‑251 (MQ=255)
                                                                                                                                                                                                                            gctgggtttGGTTTCCAGCCTTCGTGAAAGATCTCGGCCACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACACGGGTTTCGCCTGGCGATTGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGTCGCGCTTAGCGTGTTCGACGATTTATAAt              <  2:539755/243‑1 (MQ=255)
                                                                                                                                                                                                                            gctgGCGATGGTGTCCAGCCGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATaa               <  1:609692/250‑1 (MQ=255)
                                                                                                                                                                                                                             ctgGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATg             >  1:636584/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCGATGGTTTCAAGCCTGCTGGACAGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTAGAATgagg          <  2:609796/251‑2 (MQ=255)
                                                                                                                                                                                                                                         tCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgga  <  1:525456/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
AGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGA  >  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: