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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,5200GT92.1% 953.8 / 71.1 329P183P (CCG→CCTyeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base T (119/184);  total (134/195)
Fisher's exact test for biased strand distribution p-value = 9.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 329 total aligned reads displayed.

TGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTC  >  CP009273/1872374‑1872667
                                                                                                                                                  |                                                                                                                                                   
 ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATc                                                                                                                                                 >  2:2172694/1‑150 (MQ=255)
 ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATc                                                                                                                                                 >  2:1061752/1‑150 (MQ=255)
 ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATc                                                                                                                                                 <  2:4252870/150‑1 (MQ=255)
  gAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcg                                                                                                                                                <  2:3612882/150‑1 (MQ=255)
     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                             <  1:3540042/150‑1 (MQ=255)
     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                             <  2:2058853/150‑1 (MQ=255)
       aaCGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCggcggg                                                                                                                                           >  1:1811747/1‑150 (MQ=255)
       aaCGCGCGGAACAGGTTGTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCggcggg                                                                                                                                           <  2:1791536/150‑1 (MQ=255)
         cgcgcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCggcgggcg                                                                                                                                         >  2:3264711/1‑150 (MQ=255)
           cgcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAt                                                                                                                                       <  1:244107/150‑1 (MQ=255)
               gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                   >  1:2553163/1‑150 (MQ=255)
               gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                   >  2:2068078/1‑150 (MQ=255)
                  cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaa                                                                                                                                >  1:3292791/1‑150 (MQ=255)
                  cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaa                                                                                                                                >  1:4164945/1‑150 (MQ=255)
                  cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACAGGGCTACCTAATCGGCGGGAGATTTTAGaa                                                                                                                                >  2:3191263/1‑150 (MQ=255)
                   aGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaa                                                                                                                               <  1:1025855/150‑1 (MQ=255)
                   aGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaa                                                                                                                               <  1:513724/150‑1 (MQ=255)
                     gTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAAt                                                                                                                             <  1:3621546/150‑1 (MQ=255)
                      ttttACAGTTATTCGCCGAACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATc                                                                                                                            <  1:562602/150‑1 (MQ=255)
                          aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgac                                                                                                                        <  1:1932862/150‑1 (MQ=255)
                           cAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacg                                                                                                                       <  1:836998/150‑1 (MQ=255)
                            aGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacga                                                                                                                      <  1:1381629/150‑1 (MQ=255)
                            aGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacga                                                                                                                      >  1:3412673/1‑150 (MQ=255)
                             gTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacgac                                                                                                                     >  2:1980286/1‑150 (MQ=255)
                               tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGacac                                                                                                                   >  1:3984512/1‑150 (MQ=255)
                                aTTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACt                                                                                                                  <  1:1703116/150‑1 (MQ=255)
                                 ttCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTg                                                                                                                 <  1:940649/150‑1 (MQ=255)
                                 ttCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTg                                                                                                                 <  1:940637/150‑1 (MQ=255)
                                  tCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGt                                                                                                                >  2:1400925/1‑150 (MQ=255)
                                    gCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGttt                                                                                                              <  1:2791792/150‑1 (MQ=255)
                                     ccGGACTGATTGCACAATATATTCAAAAAGAGGCACGGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGtttt                                                                                                             <  1:2915003/150‑1 (MQ=255)
                                       ggACTGAGTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTcc                                                                                                           <  2:2989371/150‑1 (MQ=255)
                                          cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                        <  2:3502012/150‑1 (MQ=255)
                                          cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                        <  1:1785581/150‑1 (MQ=255)
                                            gATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGc                                                                                                      >  2:3051807/1‑150 (MQ=255)
                                              tttcacaTTATATTCAAAAAGAGTCACTTTTTGAATAGATTCGCGAAGCATATGCTGCGCTGATTGCGAAATCGTATACCGACTCGTTTAACGGGCTAACTAATAGGCGGGCGATTTTTTAAAATCTGACGACACTGTTTTCCCTCTccc                                                                                                    <  2:2962207/147‑1 (MQ=255)
                                              ttGCACTATATATTCACAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGccc                                                                                                    <  2:4415921/150‑1 (MQ=255)
                                               tGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCg                                                                                                   <  2:11803/150‑1 (MQ=255)
                                                  acaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCa                                                                                                >  1:3783812/1‑150 (MQ=255)
                                                  acaATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCa                                                                                                >  2:2762758/1‑150 (MQ=255)
                                                    aaTATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATc                                                                                              >  1:1540286/1‑150 (MQ=255)
                                                      tatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCtt                                                                                            <  2:373552/150‑1 (MQ=255)
                                                           tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                       >  1:751825/1‑150 (MQ=255)
                                                           tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                       <  1:1664195/150‑1 (MQ=255)
                                                            cAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                      <  2:4065720/150‑1 (MQ=255)
                                                             aaaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                     <  1:3017694/150‑1 (MQ=255)
                                                             aaaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                     >  2:2326357/1‑150 (MQ=255)
                                                              aaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATaa                                                                                    >  2:275074/1‑150 (MQ=255)
                                                                aaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGa                                                                                  <  2:708132/150‑1 (MQ=255)
                                                                 agagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAt                                                                                 <  1:1528981/150‑1 (MQ=255)
                                                                  gagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                                <  1:62322/150‑1 (MQ=255)
                                                                   agGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATaa                                                                               >  1:1246871/1‑150 (MQ=255)
                                                                   agGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATaa                                                                               <  2:3607902/150‑1 (MQ=255)
                                                                      cACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGa                                                                            >  1:2372798/1‑150 (MQ=255)
                                                                        cTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATc                                                                          <  1:3061351/150‑1 (MQ=255)
                                                                          ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                        <  1:2034175/150‑1 (MQ=255)
                                                                          ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                        >  2:2207950/1‑150 (MQ=255)
                                                                              gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                    >  1:325438/1‑150 (MQ=255)
                                                                              gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                    >  1:266984/1‑150 (MQ=255)
                                                                                aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                                  <  2:2615311/150‑1 (MQ=255)
                                                                                 cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                 >  1:4509691/1‑150 (MQ=255)
                                                                                  aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCg                                                                >  2:1500994/1‑150 (MQ=255)
                                                                                      gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                            <  2:2009481/150‑1 (MQ=255)
                                                                                      gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGATACCTAATCGGCGGGCTATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATTATCGCGTTTATCGATCt                                                            <  2:4389290/150‑1 (MQ=255)
                                                                                       cgcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCTACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTg                                                           <  2:1234539/150‑1 (MQ=255)
                                                                                          gAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                        >  1:332885/1‑150 (MQ=255)
                                                                                          gAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                        >  2:266729/1‑150 (MQ=255)
                                                                                           aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                       <  1:485904/150‑1 (MQ=255)
                                                                                           aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                       <  2:2626271/150‑1 (MQ=255)
                                                                                           aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                       >  2:2550332/1‑150 (MQ=255)
                                                                                                aaTGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTc                                                  >  2:4204008/1‑150 (MQ=255)
                                                                                                  tgctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCaa                                                <  2:3122795/150‑1 (MQ=255)
                                                                                                  tgctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCaa                                                <  2:2190059/150‑1 (MQ=255)
                                                                                                     tgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAtt                                             <  1:1552848/150‑1 (MQ=255)
                                                                                                        gcTGATTGCGAAATCGTATTCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAAt                                          <  1:4265245/150‑1 (MQ=255)
                                                                                                         cTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATc                                         >  2:1354809/1‑150 (MQ=255)
                                                                                                          tGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCa                                        <  2:308063/150‑1 (MQ=255)
                                                                                                             ttGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATg                                     <  2:4245478/150‑1 (MQ=255)
                                                                                                                cgcAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATc                                  <  1:3075865/150‑1 (MQ=255)
                                                                                                                   aaTCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                               <  1:1646843/150‑1 (MQ=255)
                                                                                                                   aaTCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                               <  2:3919998/150‑1 (MQ=255)
                                                                                                                     tCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtttt                             <  2:751825/150‑1 (MQ=255)
                                                                                                                     tCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtttt                             <  1:1992908/150‑1 (MQ=255)
                                                                                                                     tCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtttt                             <  2:2571313/150‑1 (MQ=255)
                                                                                                                      cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                            >  1:4075218/1‑150 (MQ=255)
                                                                                                                      cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                            <  1:3928120/150‑1 (MQ=255)
                                                                                                                      cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                            >  1:3481608/1‑150 (MQ=255)
                                                                                                                        tataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGt                          >  2:3329820/1‑150 (MQ=255)
                                                                                                                        tataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGt                          <  2:1766710/150‑1 (MQ=255)
                                                                                                                         ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                         <  2:3412450/150‑1 (MQ=255)
                                                                                                                               aCTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCataata                   <  2:1318206/150‑1 (MQ=255)
                                                                                                                                cTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAg                  >  2:582065/1‑150 (MQ=255)
                                                                                                                                 tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt                 >  2:2054374/1‑150 (MQ=255)
                                                                                                                                  cGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg                >  1:2221572/1‑150 (MQ=255)
                                                                                                                                    ttcaCCGGGCGACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGc              <  2:4447041/147‑1 (MQ=255)
                                                                                                                                    ttAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGc              <  2:3093690/150‑1 (MQ=255)
                                                                                                                                       aCCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGAt           >  2:1055696/1‑150 (MQ=255)
                                                                                                                                          gggCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTg        >  2:482325/1‑150 (MQ=255)
                                                                                                                                            gCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGtt      >  1:3941255/1‑150 (MQ=255)
                                                                                                                                                ccGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtctc  <  2:1334758/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTC  >  CP009273/1872374‑1872667

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: