Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,139,910 G→T 66.6% R1055L (CGA→CTA)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,139,9100GT66.6% 324.1 / inf 300R1055L (CGA→CTA) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (48/52);  new base T (105/95);  total (153/147)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 301 total aligned reads displayed.

CAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATTA  >  CP009273/2139765‑2140059
                                                                                                                                                 |                                                                                                                                                      
  ggTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCgg                                                                                                                                                  <  2:3226158/150‑1 (MQ=255)
   gTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCggg                                                                                                                                                 <  2:2848732/150‑1 (MQ=255)
   gTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCggg                                                                                                                                                 <  1:2586367/150‑1 (MQ=255)
    tGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGa                                                                                                                                                <  1:2959958/150‑1 (MQ=255)
      ggCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATg                                                                                                                                              <  2:1951057/150‑1 (MQ=255)
       gCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGa                                                                                                                                             <  1:3951231/150‑1 (MQ=255)
          gcgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGaaaa                                                                                                                                          <  2:3912760/150‑1 (MQ=255)
           cgATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAAc                                                                                                                                         >  1:4626838/1‑150 (MQ=255)
             aTCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCa                                                                                                                                       >  1:3443413/1‑150 (MQ=255)
              tCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCAt                                                                                                                                      >  2:3374984/1‑150 (MQ=255)
                 gCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGc                                                                                                                                   >  2:762759/1‑150 (MQ=255)
                     aTCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCggg                                                                                                                               <  1:54300/150‑1 (MQ=255)
                        ttCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGGGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCa                                                                                                                            >  2:4609171/1‑150 (MQ=255)
                        ttCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCa                                                                                                                            >  2:522282/1‑150 (MQ=255)
                         tCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAg                                                                                                                           >  1:1249713/1‑150 (MQ=255)
                          cAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGt                                                                                                                          <  2:3703948/150‑1 (MQ=255)
                          cAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGt                                                                                                                          <  2:2368228/150‑1 (MQ=255)
                           aaCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                         <  1:2023032/150‑1 (MQ=255)
                           aaCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTc                                                                                                                         >  1:4694318/1‑150 (MQ=255)
                            aCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGCTGAAAACCATCGCCGGGCCAGTCg                                                                                                                        >  2:4651382/1‑150 (MQ=255)
                            aCTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCg                                                                                                                        >  1:1122498/1‑150 (MQ=255)
                                gCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGtttt                                                                                                                    <  2:79799/150‑1 (MQ=255)
                                 cTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTa                                                                                                                   <  2:1215532/150‑1 (MQ=255)
                                  tGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGttttta                                                                                                                  >  1:1330849/1‑148 (MQ=255)
                                      aGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCtt                                                                                                              >  1:2309669/1‑150 (MQ=255)
                                       gCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGTTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTa                                                                                                             <  1:4697426/150‑1 (MQ=255)
                                        cGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAg                                                                                                            <  1:2079154/150‑1 (MQ=255)
                                        cGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAg                                                                                                            >  1:1884264/1‑150 (MQ=255)
                                         gAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                           <  2:29970/150‑1 (MQ=255)
                                         gAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                           >  1:4330397/1‑150 (MQ=255)
                                         gAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGt                                                                                                           <  1:3983005/150‑1 (MQ=255)
                                          aaTATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTg                                                                                                          <  2:1585251/150‑1 (MQ=255)
                                           atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCGTAGTGa                                                                                                         >  1:842138/1‑150 (MQ=255)
                                           atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGa                                                                                                         >  1:2734939/1‑150 (MQ=255)
                                           atatGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGa                                                                                                         >  1:765419/1‑150 (MQ=255)
                                             agggCATATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATg                                                                                                       <  2:4250037/148‑1 (MQ=255)
                                               ggCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGATGGa                                                                                                     >  1:1136476/1‑150 (MQ=255)
                                                     ttACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCt                                                                                               <  2:2100873/150‑1 (MQ=255)
                                                        ccTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTc                                                                                            <  2:3458474/150‑1 (MQ=255)
                                                         cTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCt                                                                                           <  2:3953206/150‑1 (MQ=255)
                                                          tGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTg                                                                                          <  1:2825978/150‑1 (MQ=255)
                                                             tACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGGTGGATACGCTTTCTGGcc                                                                                       >  1:497576/1‑149 (MQ=255)
                                                               cTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATc                                                                                     <  2:140786/150‑1 (MQ=255)
                                                                ttGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCg                                                                                    >  1:4758866/1‑150 (MQ=255)
                                                                 tGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCgg                                                                                   >  1:3997477/1‑150 (MQ=255)
                                                                 tGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCgg                                                                                   <  1:983721/150‑1 (MQ=255)
                                                                   aTGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCg                                                                                 >  2:1980012/1‑150 (MQ=255)
                                                                   aTGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCg                                                                                 <  1:3195375/150‑1 (MQ=255)
                                                                    tGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCGTAGTGATGGATACGCGTTCTGGCATCGGCGt                                                                                >  1:1711864/1‑150 (MQ=255)
                                                                     ggTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATTAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTc                                                                               <  2:1289226/150‑1 (MQ=255)
                                                                     ggTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGtg                                                                               <  2:3733976/150‑2 (MQ=255)
                                                                     ggTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTc                                                                               <  2:3750626/150‑1 (MQ=255)
                                                                      gTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCg                                                                              >  1:737944/1‑150 (MQ=255)
                                                                      gTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCg                                                                              >  1:3171133/1‑150 (MQ=255)
                                                                       tGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGa                                                                             <  1:4077037/150‑1 (MQ=255)
                                                                        gAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGAt                                                                            <  1:1272472/150‑1 (MQ=255)
                                                                        gAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGAt                                                                            <  1:731552/150‑1 (MQ=255)
                                                                         aGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATc                                                                           <  2:2556588/150‑1 (MQ=255)
                                                                              tGCGCCAACGTGCAGGTTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGAtt                                                                      <  1:3721050/150‑1 (MQ=255)
                                                                                 gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                   >  2:4560683/1‑150 (MQ=255)
                                                                                 gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                   <  2:4241035/150‑1 (MQ=255)
                                                                                 gcCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTCTGGCATCGGCGTCGATCTGATTTAt                                                                   >  1:1754919/1‑150 (MQ=255)
                                                                                  ccAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTGTTACCCATAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTAtt                                                                  >  1:4547469/1‑149 (MQ=255)
                                                                                           aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtga                                                         <  2:1011188/150‑1 (MQ=255)
                                                                                           aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTTGCATCGGCGTCGATCAGATTTATtgtgaggtga                                                         >  1:4235041/1‑150 (MQ=255)
                                                                                           aGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATggtgaggtga                                                         >  2:4752990/1‑150 (MQ=255)
                                                                                            gggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAt                                                        <  2:2051892/150‑1 (MQ=255)
                                                                                             ggTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGCCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGAtt                                                       <  2:1861985/150‑1 (MQ=255)
                                                                                                    tGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATg                                               >  2:3484238/1‑150 (MQ=255)
                                                                                                       tggatgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGTTACGCTTTCTGGCATCGGCGGCGATCTGATTTATGGTGAGGTGATTGCCGATTccc                                            >  1:3836636/1‑150 (MQ=255)
                                                                                                           tgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTATTGCATCGGCGTCGATATGATTTATTGTGTGGTGTTTTGCGATGCCCAAcc                                        >  1:2857659/1‑150 (MQ=255)
                                                                                                           tgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAAcc                                        <  1:2633597/150‑1 (MQ=255)
                                                                                                           tgAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAAcc                                        <  1:205014/150‑1 (MQ=255)
                                                                                                            gagaTGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCATAGTGTTGGATAAGCGTTCTGGCATCGGCGTCGATCTGTTTTATTGTGTGGTGATTGCCGATGCCCAAccc                                       >  1:2528081/1‑149 (MQ=255)
                                                                                                             agaTGCTGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGGAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATAAGCTTTCTGGCATCGGCGTCGATATGTTTTATTGTGAGGTGATTGCCGATGCCCAACCGc                                      >  2:3789228/1‑150 (MQ=255)
                                                                                                             agaTGCTGATTACGAGTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCGTTCTGGCATCGGCGTCGATGGGATTTATGGTGAGGTGTTTGCCGATGCCCAACCGc                                      >  2:3264790/1‑150 (MQ=255)
                                                                                                                tGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTgg                                   <  1:1517694/150‑1 (MQ=255)
                                                                                                                   tGATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATGTGTTTTTTGGTGTGGTGTTTGGCGATTGCCAAAcccaggttt                                 >  1:1689033/1‑143 (MQ=255)
                                                                                                                    gATTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGAttt                               <  2:1356920/150‑1 (MQ=255)
                                                                                                                     aTTACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGCTTGGCGATGCCCAACCGCTGGATTTg                              >  2:4078972/1‑150 (MQ=255)
                                                                                                                     aTTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTg                              <  1:3632277/150‑1 (MQ=255)
                                                                                                                       tACGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCt                            <  2:3118273/150‑1 (MQ=255)
                                                                                                                        aCGATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTTGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAAACGCAGGATTTGCTg                           >  2:266759/1‑150 (MQ=255)
                                                                                                                        aCGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTg                           <  2:155838/150‑1 (MQ=255)
                                                                                                                          gATTATTCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGt                         <  1:346078/150‑1 (MQ=255)
                                                                                                                          gATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCGGGATTTGCTggg                         >  2:2545362/1‑149 (MQ=255)
                                                                                                                            ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGa                       >  1:3241959/1‑150 (MQ=255)
                                                                                                                            ttattCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGa                       <  1:4032188/150‑1 (MQ=255)
                                                                                                                                tCAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtag                   <  2:2225611/150‑1 (MQ=255)
                                                                                                                                 cAGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagt                  <  1:2730821/150‑1 (MQ=255)
                                                                                                                                  aGGGGCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAAtagta                 <  1:2367845/150‑1 (MQ=255)
                                                                                                                                    gggCACGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGTATagtagt               >  1:1590370/1‑150 (MQ=255)
                                                                                                                                         cGCTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTAttt          <  1:2177931/150‑1 (MQ=255)
                                                                                                                                         cGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGGTTGATAAGCTTTCTGGCATCGGCGTCGATCTGATTTTTGGTTTGGTGATTTGCGATGACCAAACCCTGGCTTTGCTGGGGTATAGTAGtttttt          >  1:3767985/1‑150 (MQ=255)
                                                                                                                                          gCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc         <  2:4064827/150‑1 (MQ=255)
                                                                                                                                          gCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTc         <  2:1608215/150‑1 (MQ=255)
                                                                                                                                           cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcg        <  2:714413/150‑1 (MQ=255)
                                                                                                                                           cTCAGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTcg        <  2:2425418/150‑1 (MQ=255)
                                                                                                                                            tCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCGTAGTGTTGGATACGCATTCTGGCCACGGCGTCGATCTGGTTTGTGGGGTGGGGGTTGGCGCTGCCCAACCGCTGGGTTTGGTGGGGGATATTAGTTATTTcgc       >  1:2123839/1‑150 (MQ=255)
                                                                                                                                              aGCTACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTAACCTTAGTGATGGATACGCTTTCTGGCATCGGCGGCGATCGGGTTTATTGTGAGGTGTTTTCCGATTCCCAAACCCTGGTTTTGGTGGGGTATATTTTTTAtttatcga     >  1:3019272/1‑145 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCGGGATGAAAACCATCGCCGGGCCAGTCGTTTTACCCTTAGTGATGGATACGCTTTCTGGCATCGGCGTCGATCTGATTTATGGTGAGGTGATTGCCGATGCCCAACCGCTGGATTTGCTGGTGAATAGTAGTTATTTCGCGATTA  >  CP009273/2139765‑2140059

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