Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,749 T→A 9.7% T100S (ACA→TCA)  allE ← S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,7490TA 9.7% 169.8 / 9.2 62T100S (ACA→TCA) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base A (2/4);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-02

TGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAA  >  NC_000913/543456‑544038
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                  
tgtTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCt                                                                                                                                                                                                                                                                                            <  1:241507/301‑1 (MQ=255)
actggagttcagacgtgtgctcttccgatctGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCt                                                                                                                                                                                                                                                                                            <  2:1178/271‑1 (MQ=255)
  tttCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCtta                                                                                                                                                                                                                                                                                          <  2:506220/301‑2 (MQ=255)
    tCGATATAACCGTGGCTTGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCttaag                                                                                                                                                                                                                                                                                        <  2:442183/301‑4 (MQ=255)
         ctaACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTg                                                                                                                                                                                                                                                                                   <  1:298385/300‑1 (MQ=255)
           aaCCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTgaa                                                                                                                                                                                                                                                                                 >  2:298789/1‑299 (MQ=255)
             ccGTGGCTGGTACCTGGTGCAAAAGAGAGGTTATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAg                                                                                                                                                                                                                                                                               <  2:429078/301‑1 (MQ=255)
                              tGCAAACGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCAAAAGAta                                                                                                                                                                                                                                                              >  2:517073/1‑300 (MQ=255)
                                cAAAAGAGAGGATATTCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATcac                                                                                                                                                                                                                                                            <  2:372287/301‑1 (MQ=255)
                                aaaaaGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATcac                                                                                                                                                                                                                                                            <  2:390184/300‑1 (MQ=255)
                                               gCATGTTCATATCAAAACCTAACTCGTGGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTGCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTGCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTGAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTc                                                                                                                                                                                                                                             <  1:378512/301‑1 (MQ=255)
                                                cATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCa                                                                                                                                                                                                                                            >  2:602153/1‑301 (MQ=255)
                                                 ttGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATACCATCTTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCaa                                                                                                                                                                                                                                           <  2:491838/300‑1 (MQ=255)
                                                   tttCATATCAAAACCTAACTCTTTTGGCAGAAAAGCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTTGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATg                                                                                                                                                                                                                                         <  1:271912/300‑1 (MQ=255)
                                                   gTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATg                                                                                                                                                                                                                                         >  2:420691/1‑301 (MQ=255)
                                                                           gggAAGTAAAGCCAGAAGAATAACATCGTCAATGCCTTTATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTgttgg                                                                                                                                                                                                                 <  1:161463/301‑1 (MQ=255)
                                                                               aGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCGTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGGCATTAAGGCGCCTGGCGGGCAATAACGATAGCCACCTTCCCTTAAGACAAAAGTTTAGCCTTCGCAGTTGGCAGAGACATTTACAGAAAACACATACAGAAACGATTCAAAGCACTCGCCACCCACGcccacttagttac                                                                                                                                                                                                             >  1:64998/1‑291 (MQ=255)
                                                                                   aaTCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTAGCCACCGAAACCCTGTTGGTTGCCAcc                                                                                                                                                                                                         >  1:478654/1‑301 (MQ=255)
                                                                                         acagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTg                                                                                                                                                                                                   <  1:241981/300‑1 (MQ=255)
                                                                                             aaTAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCTCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGATTGCCACCGTTTTGATga                                                                                                                                                                                               >  2:121684/1‑300 (MQ=255)
                                                                                                      gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGAcc                                                                                                                                                                                      >  1:308542/1‑301 (MQ=255)
                                                                                                                 cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATcaaa                                                                                                                                                                           >  2:25009/1‑300 (MQ=255)
                                                                                                                                 ctgTTGGCTGGCATTGCCAGATGCGTGCCACGGTGCCTTGCCTTCTACCGGGAGATAGCGGCGTTTATATAGAAAGATTTGGCTGNNTTCGGCCTGGGNNNTAACAACCGTCATTNNNNGGCCTGGCGGGCAAGAAAGATAGCCACCNNCGCTTAAGGCAANTGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGc                                                                                                                                                           <  1:131561/299‑1 (MQ=255)
                                                                                                                                                    aaaaCCAGCCTCGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTg                                                                                                                                        <  1:601392/300‑1 (MQ=255)
                                                                                                                                                     aaaCCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACACACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAACTTTGGCGTGGAGAGGagggt                                                                                                                                       >  1:351728/1‑297 (MQ=255)
                                                                                                                                                       aCCAGCCATGGTTCATAGCCTTCGACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTTACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTcg                                                                                                                                     <  1:245687/301‑1 (MQ=255)
                                                                                                                                                         cAGCCACGGTGCATAGCCTTTTACCGTAACATAGCGGCGCGTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTcgcg                                                                                                                                   <  1:420128/301‑1 (MQ=255)
                                                                                                                                                         cAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTcgcg                                                                                                                                   >  2:270930/1‑301 (MQ=255)
                                                                                                                                                           gCCACGGTTCATAGCCTTCTACCGGAACATTGCGGGGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTc                                                                                                                                 <  1:270632/301‑1 (MQ=255)
                                                                                                                                                                    cATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAAAGGCACCCAGCTTTGGCGTGGAGCGGATTGTCGCGTCAAAAtttta                                                                                                                        >  2:26222/1‑300 (MQ=255)
                                                                                                                                                                      tAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGCGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAAAGGCACCCAGCTTTGGCCTGGACAGGCATATCACCTCACAATTTTCaa                                                                                                                      >  1:631181/1‑301 (MQ=255)
                                                                                                                                                                       aGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTAGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCAACGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATACTCAACAAAAGAGGCACCCAGATTTGGCGTGGAGAGGAATGTCGCGTACCAACTTTCaaa                                                                                                                     >  1:389380/1‑301 (MQ=255)
                                                                                                                                                                                ccGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTTGCGTCACAATTTTCAAAGCCCGGaat                                                                                                            <  2:308946/301‑1 (MQ=255)
                                                                                                                                                                                   taaCATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGaataat                                                                                                         <  2:352191/300‑1 (MQ=255)
                                                                                                                                                                                       aTAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTAACAATTTTCAAACCCCGAAATAATAttt                                                                                                     >  1:414488/1‑301 (MQ=255)
                                                                                                                                                                                                     ataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAACCCGGAATAATATTTTTAAACAGACCATc                                                                                       >  2:213948/1‑301 (MQ=255)
                                                                                                                                                                                                      taAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTCACAAACGTCATTAAGGAGCCTGGCGGGAAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACCGGAACGTTTCAATGACTTCGCCACCGAAGCCACATTGGGTGCCAACGTTTTGATGCAGTGTGACCAGATAAACAAAAAAATAAGAACCAAGCTTGGACGTGGAGAGGAGTGTCGAGACAAACATATCCAAGCCCGGAAAAAAATAATTTACCAAAcacacc                                                                                      >  1:369226/1‑296 (MQ=255)
                                                                                                                                                                                                         aaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCggg                                                                                   >  2:323022/1‑301 (MQ=255)
                                                                                                                                                                                                                    cagacgtgtgctcttccgatctaaaCGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGc                                                                        <  2:616528/279‑1 (MQ=255)
                                                                                                                                                                                                                       cTTCGGCCTGGGCGTTAACAAACGTCATTAATGAGCCTGGCGGGCAATAAAGATAGCCACCTTCTCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTACAGAGATCACATACAGGAACGTTTCAATACCTTCGCCACCGAAGACCTGTTATTTGCCACCGTTTTGATGCAGTGTGACCAACTAATCAACAAAAGATGCACCCAGCTTTGGCGTGGAGATGATTGTCGCGTCACAATTTTCAACTCCCGACATAATATTTTTTACCATACCATCATAGGTTAACAGTGCGaa                                                                     >  2:429829/1‑301 (MQ=255)
                                                                                                                                                                                                                              ttctgatctcaCAAACGTGATTAAGTAGCCTGTCGTGCAATATAGATTTGCACCGTCGCGTAAGGAAAATGTTTTACATTTGGCTTTGGCGGTGATATGTCCAGAGAGTACATACAGGAACGTTTCAAGTCATTCGCTACCGAAGCCCTGTGGGTTTCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCTGGTTTAACAGTGCGAAATTACCg                                                              <  1:149926/291‑1 (MQ=255)
                                                                                                                                                                                                                              cTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGAACATCCGGGGTTAACAGTGCGAAATTACCg                                                              >  2:510834/1‑301 (MQ=255)
                                                                                                                                                                                                                                gggCGTTTATAAAAGTAATTAAGGGGCCTGTCTGGCAATAAAGATAGTCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAGGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCgtg                                                            <  2:369710/301‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTAACAAACGTTATTAAGGTGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTGTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATGATATTTTTTACCAGACCATCCGGGGTTAACAGGGCGAAATTACCGTGttt                                                         <  1:516401/301‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGACAACGTTTTGATGCAATGTGACCAGATAATAAACAAAAGAGGCACCCAGCTTTAGCGTGGAGAGGATTGTCGCATCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCAACCGGGATTAACAGTACGAAATTACCgaaatcgaaa                                                     >  2:150077/1‑292 (MQ=255)
                                                                                                                                                                                                                                          aaaCGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACACAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGcagatcgtaagagcgtcgagtag                                                  >  1:615670/1‑279 (MQ=255)
                                                                                                                                                                                                                                           aaCGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCAGGAATAATCTTTTTTACCAGACCATCAGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATcc                                                 >  1:407299/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      ggggCCGGGCGGGAAAAAAGTGTAGCCACCTTCGCTGAATGCAATGTGTTTGCCTTCGGCTTTGGTAGTGATATTTCCCGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAAGTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGcc                                      <  1:322618/298‑1 (MQ=255)
                                                                                                                                                                                                                                                           gaaggtggctatctttattgCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGGGGAGAGGATTGTCGCGGCACAATTTTCAAAGCCAGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTTTTTAACAATACCACGGTTACCCAGTaa                                 >  1:534234/20‑301 (MQ=255)
                                                                                                                                                                                                                                                            tGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTaaa                                <  2:407862/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCACAGCCCGGAATAATATTTTTAACCAGACCATCCGGGGTTAACAGTGCGAAATTAACGTGTTTAACAATCGCACGGTTAGCAAGTAAAta                              >  2:210726/1‑300 (MQ=255)
                                                                                                                                                                                                                                                              gCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATc                              <  2:479255/301‑1 (MQ=255)
                                                                                                                                                                                                                                                               ggggCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCt                             <  2:65066/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaTAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTcgcg                        >  2:314997/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                        aaGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCATTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGTAATAATATTTTTTACCAGACCATCCGTGGTTAACAGTACGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTaa                    >  2:143917/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             aGCCATCTTCGCTTTAGGCAAATGTGTTGCTGTCGGTTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTGCGCCACCGACGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGAACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGt               <  1:213705/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ccACCTTCGCGTTAGGCAAATTGTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGTGATCACATACAGGAACGTTGCAATGCCTTCGCCACCGAAGCCCTGTTGGGTGCCACCGGTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGa             <  1:510162/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ttgtgttgaGGCAAATGTTTTGCCGTCGGCTTTTGCAGTGATAGTTCCCGGGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTa        <  1:26171/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        tttAAGGCAAATGTTTTTCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTa    <  2:389894/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ccgaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTa    >  2:58024/4‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         cgaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTAAGCCAGTAAATCTTCGCGGTAACCGGTAACGTTATTaaa   >  1:222935/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          aaaGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTaaa  <  1:24958/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                  
TGTTTCGATATAACCGTGGCTGGCACCTGGTGCAAAAGAGAGGATATGCATGTTCATATCAAAACCTAACTCTTTGGGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAA  >  NC_000913/543456‑544038

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