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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,872,5200GT67.7% 369.5 / inf 316P183P (CCG→CCTyeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (41/61);  new base T (109/105);  total (150/166)
Fisher's exact test for biased strand distribution p-value = 9.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 320 total aligned reads displayed.

CAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCAT  >  CP009273/1872371‑1872669
                                                                                                                                                     |                                                                                                                                                     
cAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTTATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCg                                                                                                                                                       >  1:3164841/1‑150 (MQ=255)
cAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTTATTGAACAATATATTCAAAAAGAGGCACTGGTTGAACAGCAGCGCGAAACCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCg                                                                                                                                                       >  1:3847770/1‑150 (MQ=255)
     gAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATcg                                                                                                                                                  <  2:2153634/150‑1 (MQ=255)
       gTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggc                                                                                                                                                <  1:1338714/150‑1 (MQ=255)
               gcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGAtt                                                                                                                                        <  1:3185427/150‑1 (MQ=255)
               gcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGAtt                                                                                                                                        <  1:3024532/150‑1 (MQ=255)
                cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                       <  2:2179921/150‑1 (MQ=255)
                  gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                     <  2:2937813/150‑1 (MQ=255)
                   aaCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTg                                                                                                                                    <  2:385655/150‑1 (MQ=255)
                     cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaa                                                                                                                                  <  2:4173168/150‑1 (MQ=255)
                     cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaa                                                                                                                                  >  2:93442/1‑150 (MQ=255)
                     cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaa                                                                                                                                  >  1:4501575/1‑150 (MQ=255)
                      aGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaa                                                                                                                                 >  2:4427800/1‑150 (MQ=255)
                      aGGTTTTACAGTTATTCGCCGGACTGAGTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGaaa                                                                                                                                 <  2:3994653/150‑1 (MQ=255)
                       ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaaa                                                                                                                                >  1:2863350/1‑150 (MQ=255)
                        gTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAAt                                                                                                                               >  1:4139967/1‑150 (MQ=255)
                          tttACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCt                                                                                                                             <  2:1962424/150‑1 (MQ=255)
                             aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgac                                                                                                                          >  1:1759497/1‑150 (MQ=255)
                             aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgac                                                                                                                          <  1:145191/150‑1 (MQ=255)
                              cAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacg                                                                                                                         <  2:1385699/150‑1 (MQ=255)
                              cAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacg                                                                                                                         >  2:724554/1‑150 (MQ=255)
                               aGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTgacga                                                                                                                        >  2:2417889/1‑150 (MQ=255)
                                 ttATTCGCCGGACTGATTGCACAATATATTCAAACAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGaca                                                                                                                      <  1:3362749/150‑1 (MQ=255)
                                  tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGacac                                                                                                                     <  1:3785131/150‑1 (MQ=255)
                                   agtCGCCGGACTGCTTGCACAATATATGCCAATAGAGGCACTGGTTGAACCTCGGCGCGAAGCCAATGCTGCGCTGATTGCGCTATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGacacc                                                                                                                    <  2:3191113/148‑2 (MQ=255)
                                   aTTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGccacc                                                                                                                    <  2:560314/150‑2 (MQ=255)
                                    ttCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGCTTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTg                                                                                                                   >  1:2110076/1‑150 (MQ=255)
                                        ccGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGtttt                                                                                                               <  1:1012698/150‑1 (MQ=255)
                                         cGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTc                                                                                                              >  2:2404771/1‑150 (MQ=255)
                                             cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                          <  1:4955685/150‑1 (MQ=255)
                                             cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                          <  2:5012102/150‑1 (MQ=255)
                                             cTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAAACGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTc                                                                                                          >  2:3604394/1‑150 (MQ=255)
                                              tGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                         >  2:2455439/1‑150 (MQ=255)
                                              tGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                         >  2:3549483/1‑150 (MQ=255)
                                               gATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGc                                                                                                        >  1:467246/1‑150 (MQ=255)
                                                aTTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGcc                                                                                                       >  1:1455605/1‑150 (MQ=255)
                                                 ttGCACTATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGccc                                                                                                      <  2:2660904/150‑1 (MQ=255)
                                                 ttGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGccc                                                                                                      <  2:3203846/150‑1 (MQ=255)
                                                 ttGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGccc                                                                                                      >  2:481911/1‑150 (MQ=255)
                                                  tGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCg                                                                                                     <  2:1137848/150‑1 (MQ=255)
                                                   gCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCgg                                                                                                    >  2:1487874/1‑150 (MQ=255)
                                                        atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                               <  2:1545919/150‑1 (MQ=255)
                                                        atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                               >  2:881995/1‑150 (MQ=255)
                                                        atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                               >  1:805015/1‑150 (MQ=255)
                                                         tatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCtt                                                                                              >  1:218778/1‑150 (MQ=255)
                                                             ttCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAAcc                                                                                          <  2:5106169/150‑1 (MQ=255)
                                                              tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCa                                                                                         <  1:4025804/150‑1 (MQ=255)
                                                                aaaaaGAGGCACTGGTTGAACAGCTGAGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                       >  2:650813/1‑150 (MQ=255)
                                                                  aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAg                                                                                     >  2:4329325/1‑150 (MQ=255)
                                                                  aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAg                                                                                     <  1:5113415/150‑1 (MQ=255)
                                                                   aaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGa                                                                                    >  1:4238576/1‑150 (MQ=255)
                                                                   aaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGa                                                                                    >  1:1120624/1‑150 (MQ=255)
                                                                    agagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAt                                                                                   >  2:4881874/1‑150 (MQ=255)
                                                                     gagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                                  >  2:1696994/1‑150 (MQ=255)
                                                                     gagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                                  >  2:370685/1‑150 (MQ=255)
                                                                       ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAAt                                                                                >  1:2004363/1‑150 (MQ=255)
                                                                       ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAAt                                                                                >  2:5193347/1‑150 (MQ=255)
                                                                       ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAAt                                                                                <  1:3888261/150‑1 (MQ=255)
                                                                       ggCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAAt                                                                                <  1:3127552/150‑1 (MQ=255)
                                                                         cACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGa                                                                              >  2:2102626/1‑150 (MQ=255)
                                                                          aCTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGGTTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGAt                                                                             >  1:54419/1‑150 (MQ=255)
                                                                          aCTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGAt                                                                             <  1:2925307/150‑1 (MQ=255)
                                                                          aCTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGAt                                                                             <  1:3390158/150‑1 (MQ=255)
                                                                            tGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcg                                                                           >  1:2814405/1‑150 (MQ=255)
                                                                             ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                          >  1:1255612/1‑150 (MQ=255)
                                                                             ggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgc                                                                          <  1:730825/150‑1 (MQ=255)
                                                                              gTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgcg                                                                         >  1:4828902/1‑150 (MQ=255)
                                                                              gTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATGTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgcg                                                                         >  1:980471/1‑150 (MQ=255)
                                                                              gTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCGACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcgcg                                                                         <  2:5200454/150‑1 (MQ=255)
                                                                                   aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                                    <  1:4811420/150‑1 (MQ=255)
                                                                                    cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                   >  1:1556671/1‑150 (MQ=255)
                                                                                    cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                   >  1:2565913/1‑150 (MQ=255)
                                                                                     aGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCg                                                                  <  2:149718/150‑1 (MQ=255)
                                                                                      gCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGa                                                                 <  1:1520655/150‑1 (MQ=255)
                                                                                         gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                              <  1:5027415/150‑1 (MQ=255)
                                                                                            cgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                           >  1:143411/1‑150 (MQ=255)
                                                                                             gAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                          <  2:223249/150‑1 (MQ=255)
                                                                                              aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                         <  2:4208908/150‑1 (MQ=255)
                                                                                              aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                         >  1:1811679/1‑150 (MQ=255)
                                                                                                  cAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACtt                                                     <  1:104072/150‑1 (MQ=255)
                                                                                                   aaTGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTc                                                    >  1:4711395/1‑150 (MQ=255)
                                                                                                          cgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTaa                                             >  1:2665426/1‑150 (MQ=255)
                                                                                                              gATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCaa                                         >  2:3429452/1‑150 (MQ=255)
                                                                                                                ttGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATg                                       <  2:1768000/150‑1 (MQ=255)
                                                                                                                   cgcACTCTTATTCCGTCTCGTTAACCGGGCTACCTAATCGCCGGGCGATTTTTGAATATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATATGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATc                                    <  2:2825342/150‑1 (MQ=255)
                                                                                                                   cgcAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATc                                    >  2:827447/1‑150 (MQ=255)
                                                                                                                    gcAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCg                                   <  1:1694057/150‑1 (MQ=255)
                                                                                                                         cGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                              <  1:1410462/150‑1 (MQ=255)
                                                                                                                          gTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTgg                             <  2:1195029/150‑1 (MQ=255)
                                                                                                                            ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                           >  1:1191223/1‑150 (MQ=255)
                                                                                                                            ataCCGACTCGTTAACCGGGCTACCGACTCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                           <  1:2237548/150‑1 (MQ=255)
                                                                                                                            ataCCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                           <  2:4038141/150‑1 (MQ=255)
                                                                                                                                     cGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg                  >  1:95704/1‑150 (MQ=255)
                                                                                                                                      gTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTgg                 <  1:4771272/150‑1 (MQ=255)
                                                                                                                                                cTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCTTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTGTGGTCATAAGAGTGGCGATCTGttt       >  2:2574418/1‑150 (MQ=255)
                                                                                                                                                cTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGttt       <  1:2277005/150‑1 (MQ=255)
                                                                                                                                                 tACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCt           >  1:422045/1‑145 (MQ=255)
                                                                                                                                                  aCCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtct     >  2:4559121/1‑150 (MQ=255)
                                                                                                                                                  aCCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtct     <  1:3458507/150‑1 (MQ=255)
                                                                                                                                                  aCCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTGTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtct     <  2:3970257/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCAT  >  CP009273/1872371‑1872669

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: