Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,557,473 A→G 14.2% L402P (CTC→CCC)  eutB ← ethanolamine ammonia‑lyase, large subunit, heavy chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,557,4730AG14.2% 177.7 / 4.2 73L402P (CTC→CCC) eutBethanolamine ammonia‑lyase, large subunit, heavy chain
Reads supporting (aligned to +/- strand):  ref base A (20/42);  new base G (9/2);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 5.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACGGTGTTG  >  NC_000913/2557241‑2557709
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
gTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTAACTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCatga                                                                                                                                                                                                                                          <  1:832584/239‑1 (MQ=255)
           gCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATTCTGGAATGCCGTGGGCTGGTAGTTGGGCCTGATGTAATCcccc                                                                                                                                                                                                                               >  1:60037/1‑239 (MQ=255)
                  gcccgCGGACAAATTCTTCAATCTGGTGTTGAACCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGACCGCCCGGTTGGTCAGGCGTCCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCa                                                                                                                                                                                                                        <  1:378727/235‑1 (MQ=255)
                   gCTGGGTTCATTTTCTTGAACTAGTTTTTGATCCATGACATGTTATCTCCGCGTAATCAGAAGAACAGTGACGGATCGACCGCCCGTTGGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCAt                                                                                                                                                                                                                       <  2:140097/239‑1 (MQ=255)
                       cGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTTGGCATGATGTCATCACCCAGCGGCATccc                                                                                                                                                                                                                    >  1:774207/1‑238 (MQ=255)
                       cGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGtca                                                                                                                                                                                                                                     <  1:241548/221‑1 (MQ=255)
                                 cttcAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAAAGCGGGGGTCTGGTAGTTTGGCATGATGTCCTCACCCACCGGCCACCCCACGTTGAAGt                                                                                                                                                                                                         >  1:1075607/1‑239 (MQ=255)
                                 cttcAATCTGTCTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTCGTTTAGCCTGGTGTACTCACCCAGCGGCATccccctactggagt                                                                                                                                                                                                         >  1:863213/1‑229 (MQ=255)
                                   tcAATCTGTTTGTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTg                                                                                                                                                                                                       <  1:118225/239‑1 (MQ=255)
                                                gATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTgg                                                                                                                                                                                         <  1:144955/239‑1 (MQ=255)
                                                      tGATAGGTTAACTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGtt                                                                                                                                                                                                        <  1:290161/219‑1 (MQ=255)
                                                           tGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCATGATGTCCTAACCCACCGGCATACCCATGATGGAGTTGCAGCCCGCGGTGGCGAGCAGGatc                                                                                                                                                                              >  1:1078118/1‑239 (MQ=255)
                                                             ttttctcCGGTTCATCAGAAGAACAGGGGCGGATCCGCCGCCCGTTTGGTCAGGCGCGCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGatcat                                                                                                                                                                            <  2:60037/236‑1 (MQ=255)
                                                                tctcCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGGATCCGGGGAGGGGGGGGGCCGGGAGTTGGGCATTGATGCAACCCCCAGCGGCACCCCCAGGATGTAGTTGCGGCCGGCGggggggggaagggacagcag                                                                                                                                                                         >  2:487376/1‑223 (MQ=255)
                                                                tctcCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGCCGGTATCGTGGAATGCGGGGGTGTGGGTGTTTGGGATGATGTCATACCCCCGCGGCATCCCCCAGATGTAGGTTAAGCCAGCGGGGGGGAGCggggacgccag                                                                                                                                                                         >  2:546303/1‑231 (MQ=255)
                                                                  cccGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGt                                                                                                                                                                       <  1:328034/238‑1 (MQ=255)
                                                                       ggcaccaGAAGCACAGGGGCGGCTCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATCATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCg                                                                                                                                                                  <  2:1098160/234‑1 (MQ=255)
                                                                          cacaaaagaaCCGTGAAGCATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTGCGGTGACGGGCGCAGGTTGAGTAACTGACGCGCAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTGGAGCATGATTTCATCACCCAGCGGGATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:83349/234‑1 (MQ=255)
                                                                          atcaGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:1015421/239‑1 (MQ=255)
                                                                                      gTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTACCTGACGCACAGTGGCGGTATCGTGGAAAGCGGGGGGTTGGGAGTTGGGCACTATGGCATCACCCCACGGCATCCCCCTGGTGTGATTTCACGCCTGGGGGTGCGGCAGGGACAACAGGGTTTTGTTGGGATTCTGggca                                                                                                                                                   >  2:525845/1‑239 (MQ=255)
                                                                                       tGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTg                                                                                                                                                                                          <  1:304742/199‑1 (MQ=255)
                                                                                            gATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCCTGATGTCATCACCCAGCGGCATCCCCCTGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATAAGGGTTTCGTTGAgggttttggcagagtcg                                                                                                                                             >  1:807350/1‑224 (MQ=255)
                                                                                            gATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGGCTGGGAGTTTGGCCTGATGTCATCCCCCCGCGGCATCCCCCCGATGTAGTTGCAGCCTTCGGGGGGGAGCCGGATCATCAGGGTTTTCTTGTGGTTCTGgtcagcgtca                                                                                                                                             >  1:436979/1‑239 (MQ=255)
                                                                                                   cgccccTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                                                                                                      <  1:192174/239‑1 (MQ=255)
                                                                                                   cctccgTTTGGTCTGGCCTCCGTTTGGCAAAATGCCCTGTTTTTCCAGCAGGCGTTCAAACGCCGGTGACGGGCTCAGGTTGAGTAACTGACTCATAGTTGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                                                                                                      <  2:83418/236‑1 (MQ=255)
                                                                                                               ccggctgccttttGCCCTAATGCCCCTGCTTTCCAGCCAGCGTTCAACCTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  2:863213/228‑1 (MQ=255)
                                                                                                               ccggccggctttttCCATAAGGCCCTTTCTTTCCAGCCGGCGTTCAAACCCCGGTGACGGGCGCAGGTTGAGTAACTGACGCCCAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  2:582665/228‑1 (MQ=255)
                                                                                                               cAGGCGACCGTTTGCCATTATGCCCATGCTTGCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  1:17790/238‑1 (MQ=255)
                                                                                                                     tgtgtaaagagacaggccctgggtttGCAGCCAGGGTTCAAACGCCGGGTACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGGTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                                                                                                                    <  2:431835/216‑1 (MQ=255)
                                                                                                                      cctgtttccaatcttccaaTGCTTTCCAGCCAGCGTTCAAACCCCGGTGACGGGCGCAGGTTGAGTAACTGACGCTCAGTGGCGGTACCGGGGAATGCGGTGGTCTGGTAGTTGAGAATGATGTCCACACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGAGCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                                                                                                                    <  2:191828/220‑1 (MQ=255)
                                                                                                                             ccATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:693364/188‑1 (MQ=255)
                                                                                                                              cATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTAGCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcaca                                                                                                                   <  1:394326/231‑1 (MQ=255)
                                                                                                                                   tGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTTAGGTTCTGGTCAGCGTCCGCGTGGTTGGTTTAACAGCAGTCCAAGCCCATAGAGAt                                                                                                       >  1:404065/1‑238 (MQ=255)
                                                                                                                                    gCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCCTGATGTAGTTGCAGCCTGCGGTGGCGAGCCGGCTCATCCGGGTTTTGTTGAGGTTCTGGTCAGCGTCAGCGGGGTTGGATTAAAAGCAGTCACACCCCCTAGAGATg                                                                                                      >  1:1034826/1‑238 (MQ=255)
                                                                                                                                             ttCCGGCCGGGGTTCAAACTCCGGGGGCGGGCGGAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  2:404065/239‑1 (MQ=255)
                                                                                                                                             ttCCAGGCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTTAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  1:357740/239‑1 (MQ=255)
                                                                                                                                               ccagccagCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCCGCGGCATCCCCATGATGTAGGTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCACCGGCAGCGTGGTTGGGTTAACACCAGGTCCAGCCCACAGAGCTggtccttattcaa                                                                                          >  1:431835/1‑227 (MQ=255)
                                                                                                                                                          tCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTg                                                                               <  1:546303/239‑1 (MQ=255)
                                                                                                                                                           catcgtCCGGTTACGGGCGCTGGTTGGGTTACTGAAGCACAGTGTCGGTCTCGTGAAATTCGGTGGTCTTGTAGTTGGGCATGATGTCATCACCCAGCGGCATCCCCATGACGGAGTTGCAGCCTGCGGTGGCGAGCAGGTTCATCAGGTTTTCGTGGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGa                                                                              <  2:92612/234‑1 (MQ=255)
                                                                                                                                                                tCCGGTGGCGGGGTGAGGTTGGGTAACTGACGCACAGTGGCGTGATCGTGGAATGCGGTGGTCTGGGAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGAtcttc                                                                         <  2:436979/239‑1 (MQ=255)
                                                                                                                                                                cccggtggcggggtggggTTGAGTAACTGGCGCACAGTGGCGGTATCGTGGAATGCGGTGTTCTGGTAGTTGGGCATGATGGCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGAtcttc                                                                         <  2:535617/223‑1 (MQ=255)
                                                                                                                                                                 ccGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGGCAGCGTCAGCGTGGTTGGGATAACAGCAGTCACAGCCCATAGAGATGCCGCCCAGCTTGCCCATAAAGGGAacttct                                                                        >  1:535617/1‑239 (MQ=255)
                                                                                                                                                                         ggggGCAGGTTGAGTAACTGACGCACGGTGGCGGTACCGTGGAATGCGGTGGTCTGGTAGTTGAGCTTGATGTCAACCCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGc                                                               <  2:868925/239‑1 (MQ=255)
                                                                                                                                                                          ggCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGGGCATGATGTCATCACCCAGCGGCATCCCCCTGAGGTTGTTGCAGCCCGCGGGGGCGGGCAGGAGCATCAGGGTTTTGTGGAGGTTCTGGTCAACCTCAGCGTGGTTGGGCTAACACCAGTACCACCCCCTAGGAAAGCCGCCCAGCTTTGCCCCAAAATGATCTTCTAAgacccaa                                                               >  2:1010452/1‑233 (MQ=255)
                                                                                                                                                                           gcgcAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGa                                                                              <  1:897382/223‑1 (MQ=255)
                                                                                                                                                                                     aGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCa                                                                                                                                                                                   <  1:58030/112‑1 (MQ=255)
                                                                                                                                                                                     aGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  1:615842/199‑1 (MQ=255)
                                                                                                                                                                                            gccGCACAGGGGCTGTATCGTGGAATGCGGTGGGCTTGTAGTTGAGCATGATGTCACCCCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAAGCTGGCGGTCg                                            <  2:774207/237‑1 (MQ=255)
                                                                                                                                                                                               gCACAGTGGCGGTATCGTGGGATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCCCCCAGCGGCATCCCCCTGATGTAGTTGCAGCCCGCGGTGGCGCGCAGGGtcatcaa                                                                                                                                                                         >  2:58030/1‑111 (MQ=255)
                                                                                                                                                                                               gCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCCGGATCATCCGGGTTTCGTTGAGGTTCTTGTCAGCGGCAGCGGGGTTGGTATAACATCAGTCCAAGCCCCTAGAGAAGCCGCCCAGCTTGCCCATAAAAGGATGTTCTAAGCCAGACACGATAAACTGggggggctgt                                         >  2:223420/1‑233 (MQ=255)
                                                                                                                                                                                                gacaGTGGCGGTATCGTGGAAGTCGGTGGTCTGGTAGTTGAGCATGATGTCATCCCCCAGCGGCATCCACATTATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCAGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGt                                        <  1:512210/238‑1 (MQ=255)
                                                                                                                                                                                                    gTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTTTagag                                    >  1:441185/1‑239 (MQ=255)
                                                                                                                                                                                                     tGGCGGTATCGTGGAATGCGGTGGTCTGGGAGTTGGGCATGATGTCATCCCCCAGCGGCATCCCCCTGATGGAGTTGCAGCCTTCGGTGGCGAGCCGGGTCATCAGGTTTTCGTTGAGGGTTGGGTCAGCGTTCGCGTGGGTGGGTTAACAGCAGGCACAGACCACAGAAGTGCCGGCCAGCTTGCCCCTAAAGGGAGCATCTAACCCAGCACGGCAAAACTTGCgggggtttgaaaga                                   >  2:844771/1‑227 (MQ=255)
                                                                                                                                                                                                      ggCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAgaga                                   <  1:487376/238‑1 (MQ=255)
                                                                                                                                                                                                      ggCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAgaga                                   <  1:1017782/238‑1 (MQ=255)
                                                                                                                                                                                                          gTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCCTGGTGTAGTTGCAGCCTGCGGTGGCGAGGCGGGTCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCCGGGTGGTTGGTATAACACAAGTCACAGCCCACAGAAATTCCGCCCACCTTGCACATGAAGCGAGCTTCCAAACCAGCACGGGAAAACTGGCGGtggggggagagaaaga                               >  2:1085363/1‑223 (MQ=255)
                                                                                                                                                                                                          gTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTTGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCCGCTTGCCCATAAAGTGATCTTCTAATCCAGCCCGGTGAAACTGGGGGTCGGTGTAAAGATaat                               >  1:198060/1‑236 (MQ=255)
                                                                                                                                                                                                           tttCGTGGAATGCGGGGGTTTGGTAGTTGAGCATGTTGTCATCCCCCAGGGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTc                              <  2:441185/236‑1 (MQ=255)
                                                                                                                                                                                                             tCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCgg                           <  1:243320/239‑1 (MQ=255)
                                                                                                                                                                                                              cggtggaTGCGGGGGTTTGGTAGTGGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCgg                           <  2:494260/232‑1 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCa                                                                                                                                                                                   >  2:1053564/1‑86 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCa                                                                                                                                                                                   <  1:1053564/86‑1 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACACCGGcc                         >  1:213680/1‑239 (MQ=255)
                                                                                                                                                                                                                     tggggtggtTTGGTGGTTGAGGATGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATa                    <  1:575989/235‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGc                                                 <  2:879654/209‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGc                                                 >  1:879654/1‑209 (MQ=255)
                                                                                                                                                                                                                       cGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGCCGTTGTAGAGATACTCCGGCCCAATaaa                  >  1:494260/1‑238 (MQ=255)
                                                                                                                                                                                                                              cTGGTAGTTGAGCATGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCa          <  1:844771/239‑1 (MQ=255)
                                                                                                                                                                                                                                ggTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACg        <  1:529111/239‑1 (MQ=255)
                                                                                                                                                                                                                                 gTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCCTAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATAATCCGGCCCAATAAAGCCGGCCACgg       >  1:715081/1‑239 (MQ=255)
                                                                                                                                                                                                                                    gTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCATTAAAGCCGACCACGgggt    >  1:100010/1‑239 (MQ=255)
                                                                                                                                                                                                                                      tGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGCATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGTTAATCTGCCGGTCGTTGTAAAGATACTCCGGCCCAATAAAGCCGACCAGGGGGTTg  >  1:619397/1‑239 (MQ=255)
                                                                                                                                                                                                                                       gAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCTTAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCACATAAAGCCGACCACGGGGTTg  >  1:1102453/1‑238 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
GTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCATGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCAGCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAGCACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACGGTGTTG  >  NC_000913/2557241‑2557709

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