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

GTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACGGTGTTG  >  NC_000913/2557241‑2557709
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
gTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTAACTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCa‑tga                                                                                                                                                                                                                                          <  1:832584/239‑1 (MQ=255)
           gCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATTCTGGAATGCCGTGGGCTGGTAGTTGGGCC‑TGATGTAATCcccc                                                                                                                                                                                                                               >  1:60037/1‑239 (MQ=255)
                  gcccgCGGACAAATTCTTCAATCTGGTGTTGAACCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGACCGCCCGGTTGGTCAGGCGTCCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCa                                                                                                                                                                                                                        <  1:378727/235‑1 (MQ=255)
                   gCTGGGTTCATTTTCTTGAACTAGTTTTTGATCCATGACATGTTATCTCCGCGTAATCAGAAGAACAGTGACGGATCGACCGCCCGTTGGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCAt                                                                                                                                                                                                                       <  2:140097/239‑1 (MQ=255)
                       cGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTTGGCA‑TGATGTCATCACCCAGCGGCATccc                                                                                                                                                                                                                    >  1:774207/1‑238 (MQ=255)
                       cGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGtca                                                                                                                                                                                                                                     <  1:241548/221‑1 (MQ=255)
                                 cttcAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAAAGCGGGGGTCTGGTAGTTTGGCA‑TGATGTCCTCACCCACCGGCCACCCCACGTTGAAGt                                                                                                                                                                                                         >  1:1075607/1‑239 (MQ=255)
                                 cttcAATCTGTCTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTCGTTTAGCC‑TGGTGTACTCACCCAGCGGCATccccctactggagt                                                                                                                                                                                                         >  1:863213/1‑229 (MQ=255)
                                   tcAATCTGTTTGTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTg                                                                                                                                                                                                       <  1:118225/239‑1 (MQ=255)
                                                gATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTgg                                                                                                                                                                                         <  1:144955/239‑1 (MQ=255)
                                                      tGATAGGTTAACTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGtt                                                                                                                                                                                                        <  1:290161/219‑1 (MQ=255)
                                                           tGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCA‑TGATGTCCTAACCCACCGGCATACCCATGATGGAGTTGCA‑GCCCGCGGTGGCGAGCAGGatc                                                                                                                                                                              >  1:1078118/1‑239 (MQ=255)
                                                             ttttctcCGGTTCATCAGAAGAACAGGGGCGGATCCGCCGCCCGTTTGGTCAGGCGCGCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGatcat                                                                                                                                                                            <  2:60037/236‑1 (MQ=255)
                                                                tctcCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGGATCCGGGGAGGGGGGGGGCCGGGAGTTGGGCATTGATG‑CAACCCCCAGCGGCACCCCCAGGATGTAGTTGCG‑GCCGGCGggggggggaagggacagcag                                                                                                                                                                         >  2:487376/1‑223 (MQ=255)
                                                                tctcCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGCCGGTATCGTGGAATGCGGGGGTGTGGGTGTTTGGGA‑TGATGTCATACCCCCGCGGCATCCCCCAGATGTAGGTTAA‑GCCAGCGGGGGGGAGCggggacgccag                                                                                                                                                                         >  2:546303/1‑231 (MQ=255)
                                                                  cccGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGt                                                                                                                                                                       <  1:328034/238‑1 (MQ=255)
                                                                       ggcaccaGAAGCACAGGGGCGGCTCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATCATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCg                                                                                                                                                                  <  2:1098160/234‑1 (MQ=255)
                                                                          cacaaaagaaCCGTGAAGCATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTGCGGTGACGGGCGCAGGTTGAGTAACTGACGCGCAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTGGAGCA‑TGATTTCATCACCCAGCGGGATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:83349/234‑1 (MQ=255)
                                                                          atcaGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:1015421/239‑1 (MQ=255)
                                                                                      gTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTACCTGACGCACAGTGGCGGTATCGTGGAAAGCGGGGGGTTGGGAGTTGGGCA‑CTATGGCATCACCCCACGGCATCCCCCTGGTGTGATTTCACGCCTGGGGGTGCG‑GCAGGGACAACAGGGTTTTGTTGGGATTCTGggca                                                                                                                                                   >  2:525845/1‑239 (MQ=255)
                                                                                       tGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTg                                                                                                                                                                                          <  1:304742/199‑1 (MQ=255)
                                                                                            gATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCC‑TGATGTCATCACCCAGCGGCATCCCCCTGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATAAGGGTTTCGTTGAgggttttggcagagtcg                                                                                                                                             >  1:807350/1‑224 (MQ=255)
                                                                                            gATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGGCTGGGAGTTTGGCC‑TGATGTCATCCCCCCGCGGCATCCCCCCGATGTAGTTGCA‑GCCTTCGGGGGGGAGCCGGATCATCAGGGTTTTCTTGTGGTTCTGgtcagcgtca                                                                                                                                             >  1:436979/1‑239 (MQ=255)
                                                                                                   cgccccTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                                                                                                      <  1:192174/239‑1 (MQ=255)
                                                                                                   cctccgTTTGGTCTGGCCTCCGTTTGGCAAAATGCCCTGTTTTTCCAGCAGGCGTTCAAACGCCGGTGACGGGCTCAGGTTGAGTAACTGACTCATAGTTGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGtggt                                                                                                                                      <  2:83418/236‑1 (MQ=255)
                                                                                                               ccggctgccttttGCCCTAATGCCCCTGCTTTCCAGCCAGCGTTCAACCTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  2:863213/228‑1 (MQ=255)
                                                                                                               ccggccggctttttCCATAAGGCCCTTTCTTTCCAGCCGGCGTTCAAACCCCGGTGACGGGCGCAGGTTGAGTAACTGACGCCCAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  2:582665/228‑1 (MQ=255)
                                                                                                               cAGGCGACCGTTTGCCATTATGCCCATGCTTGCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAAcag                                                                                                                           <  1:17790/238‑1 (MQ=255)
                                                                                                                     tgtgtaaagagacaggccctgggtttGCAGCCAGGGTTCAAACGCCGGGTACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGGTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                                                                                                                    <  2:431835/216‑1 (MQ=255)
                                                                                                                      cctgtttccaatcttccaaTGCTTTCCAGCCAGCGTTCAAACCCCGGTGACGGGCGCAGGTTGAGTAACTGACGCTCAGTGGCGGTACCGGGGAATGCGGTGGTCTGGTAGTTGAGAA‑TGATGTCCACACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGAGCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcac                                                                                                                    <  2:191828/220‑1 (MQ=255)
                                                                                                                             ccATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTg                                                                                                                                                               <  1:693364/188‑1 (MQ=255)
                                                                                                                              cATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTAGCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTcaca                                                                                                                   <  1:394326/231‑1 (MQ=255)
                                                                                                                                   tGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGGGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTTAGGTTCTGGTCAGCGTCCGCGTGGTTGGTTTAACAGCAGTCCAAGCCCATAGAGAt                                                                                                       >  1:404065/1‑238 (MQ=255)
                                                                                                                                    gCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCCTGATGTAGTTGCA‑GCCTGCGGTGGCGAGCCGGCTCATCCGGGTTTTGTTGAGGTTCTGGTCAGCGTCAGCGGGGTTGGATTAAAAGCAGTCACACCCCCTAGAGATg                                                                                                      >  1:1034826/1‑238 (MQ=255)
                                                                                                                                             ttCCGGCCGGGGTTCAAACTCCGGGGGCGGGCGGAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  2:404065/239‑1 (MQ=255)
                                                                                                                                             ttCCAGGCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTTAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  1:357740/239‑1 (MQ=255)
                                                                                                                                               ccagccagCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCCGCGGCATCCCCATGATGTAGGTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCACCGGCAGCGTGGTTGGGTTAACACCAGGTCCAGCCCACAGAGCTggtccttattcaa                                                                                          >  1:431835/1‑227 (MQ=255)
                                                                                                                                                          tCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTg                                                                               <  1:546303/239‑1 (MQ=255)
                                                                                                                                                           catcgtCCGGTTACGGGCGCTGGTTGGGTTACTGAAGCACAGTGTCGGTCTCGTGAAATTCGGTGGTCTTGTAGTTGGGCA‑TGATGTCATCACCCAGCGGCATCCCCATGACGGAGTTGCA‑GCCTGCGGTGGCGAGCAGGTTCATCAGGTTTTCGTGGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGa                                                                              <  2:92612/234‑1 (MQ=255)
                                                                                                                                                                tCCGGTGGCGGGGTGAGGTTGGGTAACTGACGCACAGTGGCGTGATCGTGGAATGCGGTGGTCTGGGAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGAtcttc                                                                         <  2:436979/239‑1 (MQ=255)
                                                                                                                                                                cccggtggcggggtggggTTGAGTAACTGGCGCACAGTGGCGGTATCGTGGAATGCGGTGTTCTGGTAGTTGGGCA‑TGATGGCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGAtcttc                                                                         <  2:535617/223‑1 (MQ=255)
                                                                                                                                                                 ccGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGGCAGCGTCAGCGTGGTTGGGATAACAGCAGTCACAGCCCATAGAGATGCCGCCCAGCTTGCCCATAAAGGGAacttct                                                                        >  1:535617/1‑239 (MQ=255)
                                                                                                                                                                         ggggGCAGGTTGAGTAACTGACGCACGGTGGCGGTACCGTGGAATGCGGTGGTCTGGTAGTTGAGCT‑TGATGTCAACCCCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑c                                                               <  2:868925/239‑1 (MQ=255)
                                                                                                                                                                          ggCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGGGCA‑TGATGTCATCACCCAGCGGCATCCCCCTGAGGTTGTTGCA‑GCCCGCGGGGGCGGGCAGGAGCATCAGGGTTTTGTGGAGGTTCTGGTCAACCTCAGCGTGGTTGGGCTAACACCAGTACCACCCCCTAGGAAAGCCGCCCAGCTTTGCCCCAAAATGATCTTCTAAgacccaa                                                               >  2:1010452/1‑233 (MQ=255)
                                                                                                                                                                           gcgcAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGa                                                                              <  1:897382/223‑1 (MQ=255)
                                                                                                                                                                                     aGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCa                                                                                                                                                                                   <  1:58030/112‑1 (MQ=255)
                                                                                                                                                                                     aGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCt                                                                                            <  1:615842/199‑1 (MQ=255)
                                                                                                                                                                                            gccGCACAGGGGCTGTATCGTGGAATGCGGTGGGCTTGTAGTTGAGCA‑TGATGTCACCCCCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAAGCTGGCGGTCg                                            <  2:774207/237‑1 (MQ=255)
                                                                                                                                                                                               gCACAGTGGCGGTATCGTGGGATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCCCCCAGCGGCATCCCCCTGATGTAGTTGCA‑GCCCGCGGTGGCGCGCAGGGtcatcaa                                                                                                                                                                         >  2:58030/1‑111 (MQ=255)
                                                                                                                                                                                               gCACAGTGGCGGTATCGTGGAATGCGGTGGGCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCCGGATCATCCGGGTTTCGTTGAGGTTCTTGTCAGCGGCAGCGGGGTTGGTATAACATCAGTCCAAGCCCCTAGAGAAGCCGCCCAGCTTGCCCATAAAAGGATGTTCTAAGCCAGACAC‑GATAAACTGggggggctgt                                         >  2:223420/1‑233 (MQ=255)
                                                                                                                                                                                                gacaGTGGCGGTATCGTGGAAGTCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCCCCCAGCGGCATCCACATTATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCAGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGt                                        <  1:512210/238‑1 (MQ=255)
                                                                                                                                                                                                    gTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTTTagag                                    >  1:441185/1‑239 (MQ=255)
                                                                                                                                                                                                     tGGCGGTATCGTGGAATGCGGTGGTCTGGGAGTTGGGCA‑TGATGTCATCCCCCAGCGGCATCCCCCTGATGGAGTTGCA‑GCCTTCGGTGGCGAGCCGGGTCATCAGGTTTTCGTTGAGGGTTGGGTCAGCGTTCGCGTGGGTGGGTTAACAGCAGGCACAGACCACAGAAGTGCCGGCCAGCTTGCCCCTAAAGGGAGCATCTAACCCAG‑CACGGCAAAACTTGCgggggtttgaaaga                                   >  2:844771/1‑227 (MQ=255)
                                                                                                                                                                                                      ggCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAgaga                                   <  1:487376/238‑1 (MQ=255)
                                                                                                                                                                                                      ggCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAgaga                                   <  1:1017782/238‑1 (MQ=255)
                                                                                                                                                                                                          gTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCCTGGTGTAGTTGCA‑GCCTGCGGTGGCGAGGCGGGTCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCCGGGTGGTTGGTATAACACAAGTCACAGCCCACAGAAATTCCGCCCACCTTGCACATGAAGCGAGCTTCCAAACCAG‑CACGGGAAAACTGGCGGtggggggagagaaaga                               >  2:1085363/1‑223 (MQ=255)
                                                                                                                                                                                                          gTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTTGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCCGCTTGCCCATAAAGTGATCTTCTAATCCAG‑CCCGGTGAAACTGGGGGTCGGTGTAAAGATaat                               >  1:198060/1‑236 (MQ=255)
                                                                                                                                                                                                           tttCGTGGAATGCGGGGGTTTGGTAGTTGAGCA‑TGTTGTCATCCCCCAGGGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTc                              <  2:441185/236‑1 (MQ=255)
                                                                                                                                                                                                             tCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCgg                           <  1:243320/239‑1 (MQ=255)
                                                                                                                                                                                                              cggtggaTGCGGGGGTTTGGTAGTGGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCgg                           <  2:494260/232‑1 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCa                                                                                                                                                                                   >  2:1053564/1‑86 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCa                                                                                                                                                                                   <  1:1053564/86‑1 (MQ=255)
                                                                                                                                                                                                               gTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACACCGGcc                         >  1:213680/1‑239 (MQ=255)
                                                                                                                                                                                                                     tggggtggtTTGGTGGTTGAGGA‑TGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATa                    <  1:575989/235‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGc                                                 <  2:879654/209‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGc                                                 >  1:879654/1‑209 (MQ=255)
                                                                                                                                                                                                                       cGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGCCGTTGTAGAGATACTCCGGCCCAATaaa                  >  1:494260/1‑238 (MQ=255)
                                                                                                                                                                                                                              cTGGTAGTTGAGCA‑TGATGTCATCCCCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCa          <  1:844771/239‑1 (MQ=255)
                                                                                                                                                                                                                                ggTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACg        <  1:529111/239‑1 (MQ=255)
                                                                                                                                                                                                                                 gTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCCTAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATAATCCGGCCCAATAAAGCCGGCCACgg       >  1:715081/1‑239 (MQ=255)
                                                                                                                                                                                                                                    gTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCATTAAAGCCGACCACGgggt    >  1:100010/1‑239 (MQ=255)
                                                                                                                                                                                                                                      tGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGCATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGTTAATCTGCCGGTCGTTGTAAAGATACTCCGGCCCAATAAAGCCGACCAGGGGGTTg  >  1:619397/1‑239 (MQ=255)
                                                                                                                                                                                                                                       gAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCTTAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCACATAAAGCCGACCACGGGGTTg  >  1:1102453/1‑238 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
GTCCCATTGACGCCATCACGCTGCGTACAATTTCTTCAATCTGTTTTTGATCCATGATATGTTATCTCCGCGTCATCAGAAGAACAGTGACGGATCGCCCGCCCGTTTGGTCAGGCGACCGTTTGCCATAATGCCCATGCTTTCCAGCCAGCGTTCAAACTCCGGTGACGGGCGCAGGTTGAGTAACTGACGCACAGTGGCGGTATCGTGGAATGCGGTGGTCTGGTAGTTGAGCA‑TGATGTCATCACCCAGCGGCATCCCCATGATGTAGTTGCA‑GCCTGCGGTGGCGAGCAGGATCATCAGGTTTTCGTTGAGGTTCTGGTCAGCGTCAGCGTGGTTGGTATAACAGCAGTCACAGCCCATAGAGATGCCGCTCAGCTTGCCCATAAAGTGATCTTCTAAGCCAG‑CACGGATAATCTGGCGGTCGTTGTAGAGATACTCCGGCCCAATAAAGCCGACCACGGTGTTG  >  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: ‑