Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,867,537 C→A 10.2% R414S (CGT→AGT)  yeaJ → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,867,5370CA10.2% 849.2 / 52.0 245R414S (CGT→AGT) yeaJputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (121/99);  new base A (16/9);  total (137/108)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 252 total aligned reads displayed.

AGCGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCC  >  CP009273/1867388‑1867686
                                                                                                                                                     |                                                                                                                                                     
   ggTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTa                                                                                                                                                    <  1:3176785/150‑1 (MQ=255)
    gTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTaa                                                                                                                                                   <  2:650102/150‑1 (MQ=255)
        cAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGt                                                                                                                                               >  2:1010257/1‑150 (MQ=255)
         aGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTg                                                                                                                                              >  2:3177335/1‑150 (MQ=255)
               tGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATg                                                                                                                                        <  1:801435/150‑1 (MQ=255)
                     aaTCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATcc                                                                                                                                  >  2:2502376/1‑150 (MQ=255)
                           gTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCg                                                                                                                            <  2:3507215/150‑1 (MQ=255)
                             tcttcGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGt                                                                                                                          >  2:393113/1‑150 (MQ=255)
                              cttcGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTg                                                                                                                         <  1:4524222/150‑1 (MQ=255)
                                tcGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGc                                                                                                                       >  2:3696747/1‑150 (MQ=255)
                                  ggTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGa                                                                                                                     <  1:2533477/150‑1 (MQ=255)
                                        gTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAAtt                                                                                                               >  1:3760242/1‑150 (MQ=255)
                                           tATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTg                                                                                                            <  1:2682459/150‑1 (MQ=255)
                                           tATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTg                                                                                                            >  2:430332/1‑150 (MQ=255)
                                             ttGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGca                                                                                                          <  2:4232188/150‑1 (MQ=255)
                                                 tATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGcatcat                                                                                                      >  2:2362072/1‑150 (MQ=255)
                                                  aTTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCAtt                                                                                                     <  2:3642844/150‑1 (MQ=255)
                                                  aTTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCAtt                                                                                                     >  2:2589257/1‑150 (MQ=255)
                                                  aTTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCAtt                                                                                                     <  2:3642225/150‑1 (MQ=255)
                                                   ttGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTc                                                                                                    >  1:3062606/1‑150 (MQ=255)
                                                     gACATGGACAAGTTAACGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCtt                                                                                                  <  2:3433576/150‑1 (MQ=255)
                                                     gACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCtt                                                                                                  <  2:3435353/150‑1 (MQ=255)
                                                     gACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCtt                                                                                                  <  2:3294904/150‑1 (MQ=255)
                                                         tGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCg                                                                                              >  1:2464007/1‑150 (MQ=255)
                                                         tGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCg                                                                                              >  1:3297585/1‑150 (MQ=255)
                                                          ggACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGa                                                                                             >  1:2338625/1‑150 (MQ=255)
                                                           gACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGAt                                                                                            <  2:263386/150‑1 (MQ=255)
                                                            aCAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGAtt                                                                                           >  1:1228825/1‑150 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          <  2:2583381/150‑1 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          <  1:4216730/150‑1 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          >  2:2314430/1‑150 (MQ=255)
                                                             cAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTc                                                                                          >  2:3865952/1‑150 (MQ=255)
                                                               aGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGa                                                                                        >  1:3546103/1‑150 (MQ=255)
                                                                gTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAc                                                                                       <  1:1390107/150‑1 (MQ=255)
                                                                 ttAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCTAcg                                                                                      >  2:4079841/1‑150 (MQ=255)
                                                                 ttAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcg                                                                                      >  1:1723436/1‑150 (MQ=255)
                                                                  tAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgc                                                                                     >  1:3580364/1‑150 (MQ=255)
                                                                   aaaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgcc                                                                                    >  2:4108728/1‑150 (MQ=255)
                                                                   aaaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgcc                                                                                    >  2:52937/1‑150 (MQ=255)
                                                                    aaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgccg                                                                                   >  1:3298210/1‑150 (MQ=255)
                                                                    aaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgccg                                                                                   <  2:1411272/150‑1 (MQ=255)
                                                                    aaGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGAcgccg                                                                                   <  2:1135316/150‑1 (MQ=255)
                                                                      gCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCa                                                                                 >  1:47074/1‑150 (MQ=255)
                                                                      gCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCa                                                                                 >  1:87289/1‑150 (MQ=255)
                                                                      gCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCa                                                                                 >  2:642813/1‑150 (MQ=255)
                                                                       cAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCaa                                                                                >  2:2718903/1‑150 (MQ=255)
                                                                           taaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTg                                                                            <  1:1665304/150‑1 (MQ=255)
                                                                             aatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTgca                                                                          >  2:55539/1‑150 (MQ=255)
                                                                             aatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTgca                                                                          <  2:2347066/150‑1 (MQ=255)
                                                                              atGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCACATTgcag                                                                         >  2:1363921/1‑150 (MQ=255)
                                                                               tGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCTTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTgcagc                                                                        <  1:400689/150‑1 (MQ=255)
                                                                                 acacCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGcac                                                                      >  2:30474/1‑150 (MQ=255)
                                                                                  cacCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGcaca                                                                     >  1:2445193/1‑150 (MQ=255)
                                                                                   acCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGACGATTAGTAAAAGTGATTATGCCATCCGACTCGGTTGCGATGAATTATGAATCATTCTTGTCGATTCGACGCCGCAAATTTCAGCACaa                                                                    >  1:3012543/1‑150 (MQ=255)
                                                                                    cccTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAAc                                                                   <  2:1089392/150‑1 (MQ=255)
                                                                                    cccTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAAc                                                                   <  1:227425/150‑1 (MQ=255)
                                                                                      cTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTg                                                                 <  2:1100814/150‑1 (MQ=255)
                                                                                      cTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTg                                                                 >  2:2029612/1‑150 (MQ=255)
                                                                                        cGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGcc                                                               <  2:2754837/150‑1 (MQ=255)
                                                                                         ggTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCt                                                              >  2:2445333/1‑150 (MQ=255)
                                                                                         ggTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCt                                                              >  1:3135872/1‑150 (MQ=255)
                                                                                           tcatcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGa                                                            <  1:2653087/150‑1 (MQ=255)
                                                                                            catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGaa                                                           >  2:2677106/1‑150 (MQ=255)
                                                                                             atcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAAc                                                          >  2:2243334/1‑150 (MQ=255)
                                                                                             atcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAAc                                                          >  1:4468429/1‑150 (MQ=255)
                                                                                              tcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACg                                                         >  1:1502243/1‑150 (MQ=255)
                                                                                              tcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACg                                                         <  1:3836480/150‑1 (MQ=255)
                                                                                                            tAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTc                                           >  2:1235320/1‑150 (MQ=255)
                                                                                                            tAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTc                                           >  1:4019312/1‑150 (MQ=255)
                                                                                                                  ttACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAAc                                     <  1:1099634/150‑1 (MQ=255)
                                                                                                                   tACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACa                                    <  1:3402242/150‑1 (MQ=255)
                                                                                                                    actttattaGCTCAGGCTATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGAAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACat                                   <  1:2813412/147‑1 (MQ=255)
                                                                                                                    aCGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACat                                   <  2:3240001/150‑1 (MQ=255)
                                                                                                                      gTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACatat                                 <  1:2105022/150‑1 (MQ=255)
                                                                                                                      gTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGAAACatat                                 >  1:340052/1‑150 (MQ=255)
                                                                                                                      gTTATTAGCTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACatat                                 <  2:876735/150‑1 (MQ=255)
                                                                                                                       ttattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATc                                <  2:1262229/150‑1 (MQ=255)
                                                                                                                       ttattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATc                                <  2:3542250/150‑1 (MQ=255)
                                                                                                                        tattaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcg                               <  2:281353/150‑1 (MQ=255)
                                                                                                                         attaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcgc                              >  1:2953369/1‑150 (MQ=255)
                                                                                                                          ttaGCTCAGGCGATTAAACAGTCGATTCGTAAATGTGATTATGCCTTCCGACTGGGTGGCGTTGATTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATGGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATcgcg                             <  1:4301819/150‑1 (MQ=255)
                                                                                                                              cTCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGcgccgc                         <  1:3635806/150‑1 (MQ=255)
                                                                                                                               tCAGGCGATTAAACAGTCGATTAGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCa                        >  2:3803161/1‑150 (MQ=255)
                                                                                                                                     gATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAaga                  >  2:4522962/1‑150 (MQ=255)
                                                                                                                                      aTTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAagag                 <  1:54828/150‑1 (MQ=255)
                                                                                                                                        tacaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTATTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAAGGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGAt               <  2:4172488/147‑1 (MQ=255)
                                                                                                                                        tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGAt               <  2:1610781/150‑1 (MQ=255)
                                                                                                                                        tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGAt               <  2:4442599/150‑1 (MQ=255)
                                                                                                                                         aaaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATc              <  2:759745/150‑1 (MQ=255)
                                                                                                                                          aaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCg             >  2:2224905/1‑150 (MQ=255)
                                                                                                                                            cAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGc           <  1:3305593/150‑1 (MQ=255)
                                                                                                                                            cAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGc           >  1:507342/1‑150 (MQ=255)
                                                                                                                                               tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc        <  1:2502376/150‑1 (MQ=255)
                                                                                                                                               tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc        <  1:236214/150‑1 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       <  1:958556/150‑1 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       >  1:3855915/1‑150 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       >  1:4344533/1‑150 (MQ=255)
                                                                                                                                                cGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtct       <  2:165698/150‑1 (MQ=255)
                                                                                                                                                     cGTAAAAGTGATTATGCCATCCGACTCGGTGGCTATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc  <  2:3555510/150‑1 (MQ=255)
                                                                                                                                                     cGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc  <  2:3557043/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGGTTGCAGAAACTGGTGCAATCCGGTTCTTCGGTGATGTTTATTGCTATTGACATGGACAAGTTAAAGCAAATAAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCC  >  CP009273/1867388‑1867686

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