Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,218,928 G→A 9.4% E424K (GAA→AAA)  yjbB → putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,218,9280GA 9.4% 1001.1 / 66.4 299E424K (GAA→AAA) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (138/132);  new base A (12/16);  total (150/149)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 305 total aligned reads displayed.

TCAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTG  >  CP009273/4218779‑4219077
                                                                                                                                                     |                                                                                                                                                     
tCAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCg                                                                                                                                                       <  1:851710/150‑1 (MQ=255)
 cAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGa                                                                                                                                                      <  1:1164946/150‑1 (MQ=255)
  aaCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                     >  1:1263077/1‑150 (MQ=255)
  aaCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                     >  1:1652956/1‑150 (MQ=255)
    cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAAt                                                                                                                                                   <  1:2281679/150‑1 (MQ=255)
    cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAAt                                                                                                                                                   <  2:1619015/150‑1 (MQ=255)
        ctctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCt                                                                                                                                               >  2:3862641/1‑150 (MQ=255)
           tACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGAc                                                                                                                                            <  2:1383324/150‑1 (MQ=255)
             cacCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACaa                                                                                                                                          <  2:1034200/150‑1 (MQ=255)
               ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATCTCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                        >  1:3410245/1‑150 (MQ=255)
               ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                        <  2:2920886/150‑1 (MQ=255)
                cGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATc                                                                                                                                       <  1:864831/150‑1 (MQ=255)
                 gCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCa                                                                                                                                      <  2:2548986/150‑1 (MQ=255)
                  ccATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCAc                                                                                                                                     >  2:2401439/1‑150 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    >  1:1140667/1‑150 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    >  2:3650572/1‑150 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    >  1:2684225/1‑150 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACt                                                                                                                                    <  2:1276481/150‑1 (MQ=255)
                     ttAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTgg                                                                                                                                  <  2:3273536/150‑1 (MQ=255)
                     ttAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTgg                                                                                                                                  <  2:2829153/150‑1 (MQ=255)
                       aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                >  1:3503319/1‑150 (MQ=255)
                       aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                <  1:2836006/150‑1 (MQ=255)
                       aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                >  1:498605/1‑150 (MQ=255)
                        aGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGgcag                                                                                                                               <  1:2182041/150‑1 (MQ=255)
                           tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAc                                                                                                                            >  2:1342341/1‑150 (MQ=255)
                           tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAc                                                                                                                            >  2:1799291/1‑150 (MQ=255)
                            gTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcg                                                                                                                           >  2:1978583/1‑150 (MQ=255)
                                  ggCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGc                                                                                                                     >  2:245931/1‑150 (MQ=255)
                                    cgcgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAt                                                                                                                   <  1:1393085/150‑1 (MQ=255)
                                      cgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAttt                                                                                                                 <  1:1348730/150‑1 (MQ=255)
                                        gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  1:2788769/150‑1 (MQ=255)
                                             cAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTg                                                                                                          >  2:4107061/1‑150 (MQ=255)
                                              aaaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGa                                                                                                         >  1:1096545/1‑150 (MQ=255)
                                               aaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAt                                                                                                        >  1:2401069/1‑150 (MQ=255)
                                               aaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAt                                                                                                        <  2:707457/150‑1 (MQ=255)
                                                 agaagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAggg                                                                                                      >  1:3801145/1‑150 (MQ=255)
                                                   aagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGt                                                                                                    >  2:2325283/1‑150 (MQ=255)
                                                     gagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTg                                                                                                  >  1:3066698/1‑150 (MQ=255)
                                                      agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGa                                                                                                 >  1:681332/1‑150 (MQ=255)
                                                         tGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAGGGGTTGAAgg                                                                                              >  1:1917307/1‑150 (MQ=255)
                                                         tGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCAGGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAgg                                                                                              >  1:3786875/1‑150 (MQ=255)
                                                                ggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGa                                                                                       <  2:3387349/150‑1 (MQ=255)
                                                                 gaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGAt                                                                                      <  1:1451716/150‑1 (MQ=255)
                                                                        cgcgCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGctct                                                                               <  1:4253404/150‑1 (MQ=255)
                                                                           gccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                            <  2:3763071/150‑1 (MQ=255)
                                                                           gccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                            <  1:102056/150‑1 (MQ=255)
                                                                            ccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGa                                                                           >  1:1379543/1‑150 (MQ=255)
                                                                             cgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAg                                                                          <  1:1013140/150‑1 (MQ=255)
                                                                              gcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGc                                                                         >  1:1733254/1‑150 (MQ=255)
                                                                               cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAGGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                        >  2:117001/1‑150 (MQ=255)
                                                                                tGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                       >  2:3248617/1‑150 (MQ=255)
                                                                                  ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                     <  1:2381260/150‑1 (MQ=255)
                                                                                  ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                     <  1:3707096/150‑1 (MQ=255)
                                                                                     agagATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCt                                                                  >  2:2866922/1‑150 (MQ=255)
                                                                                     agagATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTc                                                                 >  2:509058/1‑150 (MQ=255)
                                                                                       agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCa                                                                <  2:807621/150‑1 (MQ=255)
                                                                                            atcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAt                                                           <  1:4090024/150‑1 (MQ=255)
                                                                                              cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAttt                                                         <  2:2705715/150‑1 (MQ=255)
                                                                                               gAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTa                                                        <  1:3565941/150‑1 (MQ=255)
                                                                                               gAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTa                                                        <  2:489417/150‑1 (MQ=255)
                                                                                                aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                       >  1:2055326/1‑150 (MQ=255)
                                                                                                aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                       >  1:650149/1‑150 (MQ=255)
                                                                                                 aaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaaa                                                      >  1:3834609/1‑150 (MQ=255)
                                                                                                    gTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCt                                                   <  1:917533/150‑1 (MQ=255)
                                                                                                     tCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTg                                                  <  2:2073515/150‑1 (MQ=255)
                                                                                                      cGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTgg                                                 <  1:819570/150‑1 (MQ=255)
                                                                                                       gCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                <  1:849042/150‑1 (MQ=255)
                                                                                                       gCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                <  2:1116171/150‑1 (MQ=255)
                                                                                                        cTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCa                                               <  1:1675906/150‑1 (MQ=255)
                                                                                                         tCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCaa                                              >  1:224718/1‑150 (MQ=255)
                                                                                                          cAACCTTGAACAGGCCTCCGATCTCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                             <  1:1963517/150‑1 (MQ=255)
                                                                                                          cAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                             <  2:1652862/150‑1 (MQ=255)
                                                                                                             ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTc                                          >  2:989011/1‑150 (MQ=255)
                                                                                                             ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTc                                          <  1:1887747/150‑1 (MQ=255)
                                                                                                             ccTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTc                                          <  1:2132139/150‑1 (MQ=255)
                                                                                                              cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                         >  1:2552604/1‑150 (MQ=255)
                                                                                                              cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                         <  1:1008052/150‑1 (MQ=255)
                                                                                                               ttGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCgg                                        <  2:4288875/150‑1 (MQ=255)
                                                                                                                tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgt                                       >  1:3304429/1‑150 (MQ=255)
                                                                                                                tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgt                                       <  2:2862296/150‑1 (MQ=255)
                                                                                                                 gAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtg                                      <  1:4028985/150‑1 (MQ=255)
                                                                                                                  aaCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtgt                                     >  2:3051407/1‑150 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  1:776440/1‑150 (MQ=255)
                                                                                                                    cAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttc                                   <  1:1384299/150‑1 (MQ=255)
                                                                                                                        ccTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                               >  2:2086746/1‑150 (MQ=255)
                                                                                                                        ccTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                               >  1:3199256/1‑150 (MQ=255)
                                                                                                                        ccTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                               <  2:2632336/150‑1 (MQ=255)
                                                                                                                         cTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTtctc                              <  2:1470278/150‑1 (MQ=255)
                                                                                                                           ccGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTg                            <  2:1669962/150‑1 (MQ=255)
                                                                                                                               tatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGa                        <  2:2343089/150‑1 (MQ=255)
                                                                                                                                     cgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                  >  2:2202800/1‑150 (MQ=255)
                                                                                                                                          gcATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgcg             <  1:3646908/150‑1 (MQ=255)
                                                                                                                                            aTGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCt           >  2:2804635/1‑150 (MQ=255)
                                                                                                                                            aTGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCt           >  2:1688505/1‑150 (MQ=255)
                                                                                                                                            aTGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCt           <  2:888656/150‑1 (MQ=255)
                                                                                                                                            aTGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCt           >  2:3044412/1‑150 (MQ=255)
                                                                                                                                             tGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtc          >  2:3153274/1‑150 (MQ=255)
                                                                                                                                             tGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtc          <  2:681332/150‑1 (MQ=255)
                                                                                                                                                    cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt   <  2:142349/150‑1 (MQ=255)
                                                                                                                                                     gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTg  <  2:2158774/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTG  >  CP009273/4218779‑4219077

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: