Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 495,622 C→T T168I (ACC→ATC)  htpG → chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913495,6220CT95.1% 163.6 / 0.1 63T168I (ACC→ATC) htpGchaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (28/30);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CGACGAAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTT  >  NC_000913/495380‑495865
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cgacgaAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAaaga                                                                                                                                                                                                                                                 >  1:255734/1‑249 (MQ=255)
   agcaGTGATTGCCCATCTGGGGTCTAGCGCTAAATCCGGTACACAATCATTCCTCGAATCCCTGTGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCATTTTGGTGTTGGTTTCTGCGCTGCGTTTAGCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATCCACCGGTGCCGACATCATCAAAGAAGATc                                                                                                                                                                                                                                           <  2:116347/248‑1 (MQ=255)
    gaAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAaagaa                                                                                                                                                                                                                                                <  2:227886/246‑1 (MQ=255)
    gaAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAaagaa                                                                                                                                                                                                                                                >  1:227886/1‑246 (MQ=255)
                cATCTGGGGACTATCGCTAAATCCGGTACCAATTCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAACACCCAGACCATGGCGTCTTCTGGGAATCGGCTGGCGAAGGAGAATACACCGTTGCCGACATCATCAAAGCAGGTCGTGGTACTGAAAt                                                                                                                                                                                                                             >  1:270331/1‑251 (MQ=255)
                   cggggtcccaTCGCTAAATCCGGTACCAAACCAGTCCGCGAATCCCTGGGTGCTGACCATGCGAACGCCAGCCATCTGATCGGCCAGTTGTGTGGTGTGTTCTACCCTGCGTTTATCGTGGCCGACAAAGTGACCGGGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAACACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACGGaaagca                                                                                                                                                                                                                           <  2:270331/241‑4 (MQ=255)
                   cTGGGGACAATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACACAGTGACCGTGCGTACTCGAGCGGCAGGCGAAAAACCAGCAAATGGCGTCTTCTGGGAATCGGCTGGAGAAGGGGAATACACCGTTGCCGACATCATCAAAGAAGATcaggggaatgaaataac                                                                                                                                                                                                                           >  2:474909/1‑235 (MQ=255)
                      ggggcTATCGCTAAATCCGTTACCAAATCCTTCCTCGAATCCCTGGGTTCTGCCCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCCGGCGTAAAACCAGACAATGGCTTCTTCTGGTCGTCTGCTGTCGCCGGAGAATCCCCCGTTTCCGACCTCATCCAAGAAATGCGTGGTACTGCAATCACCCt                                                                                                                                                                                                                       >  2:272336/5‑251 (MQ=255)
                      gggACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGCCCAGGCGAAAGACAGACAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGAAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGCATCGGCTGGCGCAAGTGAATACACCGTTGCCCACATCATCACAGACGGTGTGGGTACTGAAATCAcccc                                                                                                                                                                                                                       >  2:41618/1‑250 (MQ=255)
                                  aaaTCCGGTCCCAAATCATTCCGCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTTATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGGTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCTCCCTGCATCTGCGTGa                                                                                                                                                                                                           <  2:228670/251‑1 (MQ=255)
                                            ccAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGAc                                                                                                                                                                                                 <  1:601480/251‑1 (MQ=255)
                                             cAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGGGAATACACCGTTGCCGACATCCTCAAAGAAGATAGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACg                                                                                                                                                                                                >  1:639876/1‑251 (MQ=255)
                                              aaaTCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACg                                                                                                                                                                                                >  1:302055/1‑250 (MQ=255)
                                              aaaTCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACg                                                                                                                                                                                                >  1:631651/1‑250 (MQ=255)
                                                   atgccTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTcc                                                                                                                                                                                          <  1:352151/248‑1 (MQ=255)
                                                           aaTCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGAc                                                                                                                                                                                  <  1:150009/251‑1 (MQ=255)
                                                              cccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCTGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTg                                                                                                                                                                                <  2:631651/250‑1 (MQ=255)
                                                              cccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTg                                                                                                                                                                                <  1:41618/250‑1 (MQ=255)
                                                                                    aGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCAATGACTGGCGCGTGCGTTCCATCATCAGCa                                                                                                                                                         <  1:145073/251‑1 (MQ=255)
                                                                                      aCAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCaaa                                                                                                                                                       <  1:164146/251‑1 (MQ=255)
                                                                                       cagccagcTGATCGGTTAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAAt                                                                                                                                                      <  1:370418/251‑1 (MQ=255)
                                                                                       cagccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAAt                                                                                                                                                      >  1:347272/1‑251 (MQ=255)
                                                                                        agccagcTGATCGGTCAGTTTGGTGTTGGTTACTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCTGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATa                                                                                                                                                     <  2:161772/251‑1 (MQ=255)
                                                                                            agcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTcc                                                                                                                                                 <  1:38497/251‑1 (MQ=255)
                                                                                              cTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCg                                                                                                                                                >  1:618541/1‑250 (MQ=255)
                                                                                                 aTCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGTCGAAAACCCAGAAAATGGCGTCTTCTGGAAATCGGCTGGCGAAGGTGAATACAACGGTGCCGACATCATCAAAGAACATCGTGGTACTGAAATCACCCAGCCTCTGCGTGAAGGCGGAGACGAGTTCCTCGATTACTGGCGGGTTCGTTCCAACATAAGCAAATACTcctgccc                                                                                                                                            >  2:163220/1‑246 (MQ=255)
                                                                                                    ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCatat                                                                                                                                         <  2:255734/251‑1 (MQ=255)
                                                                                                     gTCAGTTTGGTGTTGGTTTCGCCTCTGCGGTTAGCGTGGCCGACAATGTTACCGTGAGTAGTCGCGCGGTAGGAGAAAATCCAGAAAATGGAGTTTGCGGGGAATCGGCTGGCGAATGTGAATACACCGTTGCCGACATCATCAAAGAAGAGCGTGGTACTGAAATCATCCTGTATCTTCGTGAAGGCTAAGACGAGTTCCTCGATGACTGGCGCGTGCGGGCCATCAGAAGCATATACGCCGattatatc                                                                                                                                        <  2:402305/251‑8 (MQ=255)
                                                                                                          tttGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGGTGCCGGCATTATCAAAGAAGATCGGGGTAATGAAATCACCCTGCATCTGCGTGAAGGCGAAAACGAGTTCCTCGATGACCGGCGCGTGCGTTCCATCCTCAGGCAATAATCCGACCataacccgcc                                                                                                                                   >  2:21011/1‑244 (MQ=255)
                                                                                                            tGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGc                                                                                                                                 >  1:205046/1‑251 (MQ=255)
                                                                                                                    gTTTCTACTCTGCGTTTATCGTTGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTagag                                                                                                                         <  2:250499/251‑1 (MQ=255)
                                                                                                                          aCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACGCCGACCATATCGCGCTGCCGGTAGAGATc                                                                                                                      <  1:137607/248‑1 (MQ=255)
                                                                                                                               gCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAAc                                                                                                               <  1:347579/250‑1 (MQ=255)
                                                                                                                                  tttATCGTGGCCGACAAAGTGACCGTTCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGGGAATACACCGTTGCCGACATCATCAAAAAAGATCGGGGTAATGACAACACCCCGCATCAGCGTGAAGGCGAAAACGAGTTACGCGATGACATGCGCGGGCGTTCCCACATCAGAAAATAATACGCCCATATCACGCTGCCggcggtgagcgaaaaacacg                                                                                                            >  2:7041/1‑232 (MQ=255)
                                                                                                                                   tgaTCGTGGACGACAAAGTGACCGGGCGTACCCGCGCGGCAGGCGAAAAACCAGAAAATGGTGGCCTCTGAGAATCGGATGGCGATTGTTAATACACCTTGTCCGACATCAGCAAAGAAGATCGTGGTACTGAAATCCCCCCTCATCTGCGTGCAGGCGTAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACGCCGACCATATCGCGCTGCCGGTTGAGATCGtaaatcttga                                                                                                           <  2:6624/248‑7 (MQ=255)
                                                                                                                                    tATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGcc                                                                                                                                <  2:291039/228‑1 (MQ=255)
                                                                                                                                    tATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGcc                                                                                                                                >  1:291039/1‑228 (MQ=255)
                                                                                                                                     aTCGTGGCCGACAAAGTGACCGTGCGTCCTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGCGATCGGAAAACGCGa                                                                                                           >  1:163358/1‑248 (MQ=255)
                                                                                                                                       cGTGGCCGACAAAGTGACCGTGCGTACTCGTGCGGCAGGCGAAAAACCAGAAAATGGAGTATTGTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTACTCGATGACTGGCGCGTGCGTTACATCATCAGCAAATACTCCGACCATAGCGCGCTGCCGGTAGAGATCGAAAAACGCGAagag                                                                                                       <  1:102769/250‑1 (MQ=255)
                                                                                                                                         gcgctgaaaaaGTGCCCGTGTGTGCTCGCGCGGCAGGCTAAAAAACAGAAAAGGGCGTCTGCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGAACGTGGTACTGAAATCTCCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACGGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATAGCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGaaa                                                                                                    <  2:400545/243‑1 (MQ=255)
                                                                                                                                           gCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGAGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGaa                                                                                                     >  1:606931/1‑248 (MQ=255)
                                                                                                                                                              gTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGCGAATAAACCGTTGCCGACATCATCAAAGAAGATCGTGGTAATGAAATCAACCTGCATCTGCGTGGAGGCGTAGACGAGTTCCTCGGTGACTGGCGAGTGCGTTCCATCATCGGCAAATACTCCGACCATTTCGCGCCGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAACCCGTTaccacc                                                                               >  2:334159/1‑246 (MQ=255)
                                                                                                                                                                acgcgcgcgGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCt                                                                              <  1:71642/247‑1 (MQ=255)
                                                                                                                                                                aCTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTg                                                                             <  2:347272/251‑1 (MQ=255)
                                                                                                                                                                     cgcgGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTggg                                                                           >  1:604477/1‑248 (MQ=255)
                                                                                                                                                                     cgcgGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTggg                                                                           >  1:495263/1‑248 (MQ=255)
                                                                                                                                                                     cgcgGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGgaga                                                                        <  1:341030/251‑1 (MQ=255)
                                                                                                                                                                          cAGGCGAAAAACAAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGaaagtc                                                                   >  1:210660/1‑248 (MQ=255)
                                                                                                                                                                                  aaaCCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAAACGGCGAAACCGTTATCTCCTGGGGGAAAATCAACaaaac                                                           >  2:195413/1‑249 (MQ=255)
                                                                                                                                                                                        gAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGgc                                                     <  1:487128/251‑1 (MQ=255)
                                                                                                                                                                                         aaaaTGGCGTCTTCTGGGAATCGGCTGGCCAAGGTGAATAAACCGTTGCCGACATCATAAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGCGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGgcg                                                    >  1:241706/1‑251 (MQ=255)
                                                                                                                                                                                           aattgCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGc                                                                                                                                 <  2:395941/168‑1 (MQ=255)
                                                                                                                                                                                           aaTGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGc                                                                                                                                 >  1:395941/1‑172 (MQ=255)
                                                                                                                                                                                             tgggtttatctgggatgcggatgctgAAGGTGAATCCTCCGTTGCCGCCATCATCAAATAAGATCGTGGTACTTAAATCCCCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATCCTCCGACCATATCGCGCTTCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAACTCACCAAAGCGCGGTCGCtg                                                 <  2:654768/225‑1 (MQ=255)
                                                                                                                                                                                                       tGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCAGCATCAGCAAATACTCCGACCACATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTcc                                                                               <  1:460837/210‑1 (MQ=255)
                                                                                                                                                                                                        gcgAATCGGCTGGCGAAGGTAAATACACCGTTGCCGTCATCATCAAAGAAGATCGCGGGACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGCCGAGTTCCTCGATGACTGGCGCGTTCGTTCCATCCACAACAAATACTCCGACCCTTTCGCGCTGAAGGTAGAACTCGAAAAACGCGTAGAGAAAGACGGCGAAAACGTTATCGCCTGGGAGACGATCAACAAAGCGCAGACGCTGTGGACTCGTAAc                                     >  2:488802/3‑251 (MQ=255)
                                                                                                                                                                                                                            gAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGAGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAACGACGGCGAAACCGTTATCTCCTGGGAGAAACTCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATg                  >  2:241749/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cacCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAag               >  1:162411/1‑248 (MQ=255)
                                                                                                                                                                                                                                    cGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAGTACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTAc          >  1:204266/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACa         >  2:141187/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAag      <  1:195413/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cGACATCATCAAAGAAGATCTTGGTACTGAAATCACCCTGCATCTGCGGGAAGTCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGt   <  1:47871/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gACATCAGCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCGCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGATGAGTACAAAGAGtt  <  1:450185/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
CGACGAAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTT  >  NC_000913/495380‑495865

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