Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 328,436 A→T 7.6% N119Y (AAC→TAC)  yahA → c‑di‑GMP‑specific phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273328,4360AT 7.6% 801.0 / 45.5 250N119Y (AAC→TAC) yahAc‑di‑GMP‑specific phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (143/88);  new base T (11/8);  total (154/96)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 260 total aligned reads displayed.

CAGAGCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTT  >  CP009273/328289‑328585
                                                                                                                                                   |                                                                                                                                                     
cAGAGCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAAc                                                                                                                                                     <  2:508514/150‑1 (MQ=255)
  gagCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCa                                                                                                                                                   >  1:1285046/1‑150 (MQ=255)
   agCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCAc                                                                                                                                                  <  1:2188209/150‑1 (MQ=255)
    gCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACg                                                                                                                                                 <  1:1828988/150‑1 (MQ=255)
     cGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGa                                                                                                                                                >  2:1381757/1‑150 (MQ=255)
       atatTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCCTATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAAt                                                                                                                                              >  1:3022680/1‑150 (MQ=255)
       atatTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAATACCACGAAt                                                                                                                                              <  2:1575056/150‑1 (MQ=255)
       atatTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAAt                                                                                                                                              >  1:1727696/1‑150 (MQ=255)
       atatTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAAt                                                                                                                                              <  1:639197/150‑1 (MQ=255)
             cTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCGTATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAAc                                                                                                                                        >  2:4345545/1‑150 (MQ=255)
                tCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAATACCACGAATTCAAACCGt                                                                                                                                     >  1:1892887/1‑150 (MQ=255)
                 cTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTg                                                                                                                                    >  1:1232059/1‑150 (MQ=255)
                      gcgATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATcc                                                                                                                               <  1:4268330/150‑1 (MQ=255)
                       cgATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCa                                                                                                                              >  1:4165474/1‑150 (MQ=255)
                         atatTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAAc                                                                                                                            <  2:1538636/150‑1 (MQ=255)
                         atatTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAAc                                                                                                                            >  1:2153703/1‑150 (MQ=255)
                          tatTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAAcc                                                                                                                           <  1:1214431/150‑1 (MQ=255)
                           atTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAATACCACGAATTCAAACCGTGGATCCAACCg                                                                                                                          <  2:398277/150‑1 (MQ=255)
                           atTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCg                                                                                                                          >  1:995339/1‑150 (MQ=255)
                           atTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCg                                                                                                                          >  1:2020682/1‑150 (MQ=255)
                           atTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCg                                                                                                                          >  2:2025289/1‑150 (MQ=255)
                               tctttcAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGttt                                                                                                                      >  2:865022/1‑150 (MQ=255)
                                  ttcAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCt                                                                                                                   <  1:3690103/150‑1 (MQ=255)
                                  ttcAGTAAAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATAATAATGATAATCACTATCACCATATCGTCACGCCTGAAGCAATCAGACTGCAGTTGGAATACCACGAATTCAAACCGTGGATCCAACCGGTTTTCt                                                                                                                   >  2:2810227/1‑150 (MQ=255)
                                    cAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTgc                                                                                                                 >  1:574555/1‑150 (MQ=255)
                                    cAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTgc                                                                                                                 <  2:4259650/150‑1 (MQ=255)
                                     aGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTgcg                                                                                                                <  2:1143446/150‑1 (MQ=255)
                                     aGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTgcg                                                                                                                >  1:1080831/1‑150 (MQ=255)
                                       tACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTgcgcg                                                                                                              >  2:576816/1‑150 (MQ=255)
                                         cAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCa                                                                                                            <  2:1136530/150‑1 (MQ=255)
                                            tCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGAc                                                                                                         >  1:2385765/1‑150 (MQ=255)
                                            tCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGAc                                                                                                         <  2:4068315/150‑1 (MQ=255)
                                             ccGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACt                                                                                                        >  1:498994/1‑150 (MQ=255)
                                              cGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTg                                                                                                       >  2:241281/1‑150 (MQ=255)
                                               ggAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTgg                                                                                                      >  2:468087/1‑150 (MQ=255)
                                                  gaTCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGt                                                                                                   <  1:2055222/150‑1 (MQ=255)
                                                  gaTCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGt                                                                                                   <  1:3961333/150‑1 (MQ=255)
                                                      atatCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTg                                                                                               <  1:3296760/150‑1 (MQ=255)
                                                       tatCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGa                                                                                              >  1:3520135/1‑150 (MQ=255)
                                                        atCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGAc                                                                                             <  1:3772286/150‑1 (MQ=255)
                                                              cACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCtg                                                                                       >  1:4450176/1‑150 (MQ=255)
                                                                cGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCtgtg                                                                                     <  1:2096418/150‑1 (MQ=255)
                                                                cGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCtgtg                                                                                     >  1:931616/1‑150 (MQ=255)
                                                                  gggAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAg                                                                                   <  1:4077370/150‑1 (MQ=255)
                                                                    gAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGt                                                                                 <  2:1301537/150‑1 (MQ=255)
                                                                      gTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGc                                                                               >  2:3002349/1‑150 (MQ=255)
                                                                      gTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGc                                                                               >  2:682884/1‑150 (MQ=255)
                                                                       tAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCt                                                                              >  1:4292161/1‑150 (MQ=255)
                                                                          tAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGt                                                                           <  1:3223027/150‑1 (MQ=255)
                                                                             tCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCg                                                                        <  2:3859050/150‑1 (MQ=255)
                                                                             tCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCg                                                                        <  2:3475848/150‑1 (MQ=255)
                                                                              cacaGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGc                                                                       >  2:3197706/1‑150 (MQ=255)
                                                                               acaGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCt                                                                      >  2:1096924/1‑150 (MQ=255)
                                                                                       atTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACAt                                                              <  1:1617159/150‑1 (MQ=255)
                                                                                       atTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACAt                                                              >  2:2936820/1‑150 (MQ=255)
                                                                                           aTGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCcac                                                          <  2:914707/150‑1 (MQ=255)
                                                                                             gATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACaa                                                        >  2:1571805/1‑150 (MQ=255)
                                                                                               tAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAAc                                                      >  1:2645027/1‑150 (MQ=255)
                                                                                               tAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAAc                                                      <  2:3273574/150‑1 (MQ=255)
                                                                                                     cTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAAt                                                >  2:4139305/1‑150 (MQ=255)
                                                                                                      tATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGTAACATCCACAAACGGGAAtt                                               >  1:1337289/1‑150 (MQ=255)
                                                                                                      tATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                               >  1:4384710/1‑150 (MQ=255)
                                                                                                         cACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACTGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                            >  2:234961/1‑150 (MQ=255)
                                                                                                         cACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                            >  1:1224815/1‑150 (MQ=255)
                                                                                                          aCCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                           >  1:809762/1‑150 (MQ=255)
                                                                                                           ccATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAATACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                          >  2:3528850/1‑150 (MQ=255)
                                                                                                           ccATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                          >  1:1771949/1‑150 (MQ=255)
                                                                                                            cATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCa                                         >  2:1259619/1‑150 (MQ=255)
                                                                                                             atatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                        <  2:722181/150‑1 (MQ=255)
                                                                                                             atatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                        >  1:2826385/1‑150 (MQ=255)
                                                                                                             atatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                        >  2:3614162/1‑150 (MQ=255)
                                                                                                             atatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                        >  2:4186452/1‑150 (MQ=255)
                                                                                                              tatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                       <  1:411800/150‑1 (MQ=255)
                                                                                                              tatCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                       >  2:891966/1‑150 (MQ=255)
                                                                                                               atCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGCTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCg                                      <  2:3022680/150‑1 (MQ=255)
                                                                                                               atCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCg                                      <  1:3923415/150‑1 (MQ=255)
                                                                                                                 cGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAATACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGa                                    >  2:4054950/1‑150 (MQ=255)
                                                                                                                  gTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGAt                                   <  2:1283838/150‑1 (MQ=255)
                                                                                                                    cACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGACTTATCCCACCGGATCa                                 >  2:3134721/1‑150 (MQ=255)
                                                                                                                    cACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCa                                 <  1:1012016/150‑1 (MQ=255)
                                                                                                                        ccTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                             >  1:2260348/1‑150 (MQ=255)
                                                                                                                        ccTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                             >  1:1125551/1‑150 (MQ=255)
                                                                                                                        ccTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                             >  2:2624635/1‑150 (MQ=255)
                                                                                                                        ccTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                             <  2:1727814/150‑1 (MQ=255)
                                                                                                                          tGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAt                           <  2:2797004/150‑1 (MQ=255)
                                                                                                                           gAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                          >  1:1102537/1‑150 (MQ=255)
                                                                                                                           gAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                          >  2:1948703/1‑150 (MQ=255)
                                                                                                                                cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                     <  2:438361/150‑1 (MQ=255)
                                                                                                                                cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                     >  1:3346544/1‑150 (MQ=255)
                                                                                                                                cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                     <  2:2153703/150‑1 (MQ=255)
                                                                                                                                cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                     <  2:2137077/150‑1 (MQ=255)
                                                                                                                                 aTCAGTCTGGCGTTGGAATACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTg                    >  2:2709776/1‑150 (MQ=255)
                                                                                                                                 aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTg                    <  2:2061850/150‑1 (MQ=255)
                                                                                                                                  tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTgg                   <  1:897594/150‑1 (MQ=255)
                                                                                                                                       cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg              >  2:2715680/1‑150 (MQ=255)
                                                                                                                                        tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGt             >  1:2805229/1‑150 (MQ=255)
                                                                                                                                         ggCGTTGGAATACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc            <  1:619558/150‑1 (MQ=255)
                                                                                                                                         ggCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc            >  1:917945/1‑150 (MQ=255)
                                                                                                                                                 aaaaCCACGAATTCACACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTAGCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc    <  1:3255783/150‑1 (MQ=255)
                                                                                                                                                   aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCtt  <  1:4107287/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGAGCGATATTACTTTCTGGCGCGATATTTTCTTTCAGTACAATCCGGAGATCATATCCGCCACGGGGAGTAATAGTCACAGATATATTAATGATAATCACTATCACCATATCGTCACGCCTGAAGCCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTT  >  CP009273/328289‑328585

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: ‑