Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,872,520 G→T 88.4% P183P (CCG→CCT yeaP → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,5200GT88.4% 759.3 / 91.2 267P183P (CCG→CCTyeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (20/11);  new base T (129/107);  total (149/118)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 268 total aligned reads displayed.

GGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCT  >  CP009273/1872375‑1872666
                                                                                                                                                 |                                                                                                                                                  
ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATc                                                                                                                                                <  1:342310/150‑1 (MQ=255)
 gAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcg                                                                                                                                               <  1:2977735/150‑1 (MQ=255)
  aGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcgg                                                                                                                                              >  2:3474491/1‑150 (MQ=255)
   gTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggc                                                                                                                                             <  1:187409/150‑1 (MQ=255)
           gcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtt                                                                                                                                     <  2:2440788/150‑1 (MQ=255)
           gcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGAtt                                                                                                                                     <  1:94883/150‑1 (MQ=255)
            cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                    >  1:4170879/1‑150 (MQ=255)
            cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                    >  2:1269/1‑150 (MQ=255)
             ggAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtttt                                                                                                                                   <  2:543474/150‑1 (MQ=255)
             ggAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtttt                                                                                                                                   <  2:3575958/150‑1 (MQ=255)
              gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTCCCGAATCGGCGGGCGAtgtgt                                                                                                                                  >  2:2367176/1‑146 (MQ=255)
              gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                  >  1:1776285/1‑150 (MQ=255)
                 cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaa                                                                                                                               >  1:3010609/1‑150 (MQ=255)
                   ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaaa                                                                                                                             >  1:2920385/1‑150 (MQ=255)
                   ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaaaa                                                                                                                             <  1:2527962/150‑1 (MQ=255)
                     ttttACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATc                                                                                                                           <  1:3122198/150‑1 (MQ=255)
                      tttACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCt                                                                                                                          <  1:3241888/150‑1 (MQ=255)
                        tACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTga                                                                                                                        <  1:235142/150‑1 (MQ=255)
                         aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgac                                                                                                                       <  2:730651/150‑1 (MQ=255)
                           aGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTgacga                                                                                                                     >  1:205829/1‑150 (MQ=255)
                            gTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCGACCTAATCGGCGGGCGATTTTTGAAAATCTgacgac                                                                                                                    <  2:623838/150‑1 (MQ=255)
                             ttATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGaca                                                                                                                   >  1:1656392/1‑150 (MQ=255)
                                      ggACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTcc                                                                                                          >  2:4335300/1‑150 (MQ=255)
                                        aCTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCt                                                                                                        >  1:1149220/1‑150 (MQ=255)
                                        aCTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCt                                                                                                        >  2:2946562/1‑150 (MQ=255)
                                         cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                       <  1:3611108/150‑1 (MQ=255)
                                          tGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                      >  1:2335460/1‑150 (MQ=255)
                                              tGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCg                                                                                                  <  2:1453524/150‑1 (MQ=255)
                                              tGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCg                                                                                                  >  2:525773/1‑150 (MQ=255)
                                                cacaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGc                                                                                                >  1:1254574/1‑150 (MQ=255)
                                                cacaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGc                                                                                                >  1:1254509/1‑150 (MQ=255)
                                                 acaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCa                                                                                               <  2:3492467/150‑1 (MQ=255)
                                                 acaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCa                                                                                               <  2:709127/150‑1 (MQ=255)
                                                  caATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCAt                                                                                              >  2:2686784/1‑150 (MQ=255)
                                                    atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                            >  2:796791/1‑150 (MQ=255)
                                                     tatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCtt                                                                                           <  2:3606447/150‑1 (MQ=255)
                                                      atatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTa                                                                                          >  2:3976654/1‑150 (MQ=255)
                                                       tatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTaa                                                                                         >  1:1364320/1‑150 (MQ=255)
                                                        atTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAAc                                                                                        >  1:3869185/1‑150 (MQ=255)
                                                        atTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAAc                                                                                        >  1:3202879/1‑150 (MQ=255)
                                                          tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                      <  2:2314138/150‑1 (MQ=255)
                                                          tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                      <  2:2742065/150‑1 (MQ=255)
                                                          tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                      <  2:209488/150‑1 (MQ=255)
                                                           cAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCTATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                     <  1:150464/150‑1 (MQ=255)
                                                           cAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                     <  1:2619256/150‑1 (MQ=255)
                                                           cAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                     <  1:356018/150‑1 (MQ=255)
                                                            aaaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                    >  1:832092/1‑150 (MQ=255)
                                                            aaaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                    >  1:2207561/1‑150 (MQ=255)
                                                             aaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATaa                                                                                   >  1:1982705/1‑150 (MQ=255)
                                                              aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAg                                                                                  >  2:1148961/1‑150 (MQ=255)
                                                               aaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGa                                                                                 >  2:193847/1‑150 (MQ=255)
                                                               aaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGa                                                                                 >  2:2493884/1‑150 (MQ=255)
                                                                agagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAt                                                                                <  1:100464/150‑1 (MQ=255)
                                                                 gagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                               <  1:907914/150‑1 (MQ=255)
                                                                   ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAAt                                                                             >  2:452735/1‑150 (MQ=255)
                                                                       cTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTAAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAATATAATGATc                                                                         >  2:4491152/1‑150 (MQ=255)
                                                                        tGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcg                                                                        >  1:2781626/1‑150 (MQ=255)
                                                                         ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                       >  1:4382102/1‑150 (MQ=255)
                                                                         ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                       >  2:2889440/1‑150 (MQ=255)
                                                                                cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                >  2:3883830/1‑150 (MQ=255)
                                                                                cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                <  2:2164368/150‑1 (MQ=255)
                                                                                 aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCg                                                               >  1:3493530/1‑150 (MQ=255)
                                                                                 aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCg                                                               <  1:3639672/150‑1 (MQ=255)
                                                                                      cgcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTg                                                          >  2:3690579/1‑150 (MQ=255)
                                                                                      cgcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGAGCTg                                                          >  1:3371251/1‑150 (MQ=255)
                                                                                       gcgACGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTgg                                                         <  2:3241530/150‑1 (MQ=255)
                                                                                        cgAAGCCAATGCTGCGCTGATTGCGCAATCATATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                        >  1:3827961/1‑150 (MQ=255)
                                                                                                    tgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAtt                                            >  2:4427658/1‑150 (MQ=255)
                                                                                                    tgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAtt                                            >  2:4358070/1‑150 (MQ=255)
                                                                                                     gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTa                                           <  2:2392590/150‑1 (MQ=255)
                                                                                                     gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGACAATGATCGCGTTTATCGATCTGGATAACTTCAAATTa                                           <  2:1707172/150‑1 (MQ=255)
                                                                                                        cTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATc                                        >  2:1247254/1‑150 (MQ=255)
                                                                                                         tGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCa                                       >  1:3973129/1‑150 (MQ=255)
                                                                                                         tGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCa                                       >  2:4216867/1‑150 (MQ=255)
                                                                                                          gATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCaa                                      >  1:1291128/1‑150 (MQ=255)
                                                                                                           aTTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAAt                                     <  1:554018/150‑1 (MQ=255)
                                                                                                                     cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                           <  2:1028440/150‑1 (MQ=255)
                                                                                                                      gTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTgg                          <  2:751582/150‑1 (MQ=255)
                                                                                                                        ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                        <  2:2263141/150‑1 (MQ=255)
                                                                                                                        ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                        >  1:716035/1‑150 (MQ=255)
                                                                                                                        ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                        >  2:2071420/1‑150 (MQ=255)
                                                                                                                          aCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCat                      <  2:11302/150‑1 (MQ=255)
                                                                                                                           ccGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCata                     <  2:2189380/150‑1 (MQ=255)
                                                                                                                            cGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGGCataa                    >  2:3053856/1‑150 (MQ=255)
                                                                                                                             gACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataat                   >  2:3144838/1‑150 (MQ=255)
                                                                                                                             gACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATAGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataat                   >  1:162786/1‑150 (MQ=255)
                                                                                                                              aCTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataata                  <  1:4389114/150‑1 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAg                 >  1:3728255/1‑150 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAg                 <  1:1542574/150‑1 (MQ=255)
                                                                                                                               cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTGTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAg                 <  1:2160021/150‑1 (MQ=255)
                                                                                                                                tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt                <  1:2681677/150‑1 (MQ=255)
                                                                                                                                 cGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCCTAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg               >  1:4259302/1‑150 (MQ=255)
                                                                                                                                 cGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg               >  1:4042489/1‑150 (MQ=255)
                                                                                                                                  gTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTgg              >  2:2767164/1‑150 (MQ=255)
                                                                                                                                   ttAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGc             <  2:137380/150‑1 (MQ=255)
                                                                                                                                    tAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCg            >  2:4472144/1‑150 (MQ=255)
                                                                                                                                        cGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCt        >  1:145127/1‑150 (MQ=255)
                                                                                                                                        cGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCt        >  2:1511946/1‑150 (MQ=255)
                                                                                                                                          ggCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGt      >  2:3397712/1‑150 (MQ=255)
                                                                                                                                              aCCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtct  <  1:120433/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCT  >  CP009273/1872375‑1872666

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: