Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,218,928 G→A 40.9% 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,9280GA40.9% 227.9 / inf 281E424K (GAA→AAA) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (82/84);  new base A (60/55);  total (142/139)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 283 total aligned reads displayed.

CAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTT  >  CP009273/4218780‑4219076
                                                                                                                                                    |                                                                                                                                                    
 aaCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGAGATCGTCGAGCGCATGGGCAGCGaa                                                                                                                                                    >  2:576768/1‑150 (MQ=255)
    gTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAAtt                                                                                                                                                 >  1:286419/1‑150 (MQ=255)
    gTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAAtt                                                                                                                                                 >  2:4022707/1‑150 (MQ=255)
    gTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAAtt                                                                                                                                                 <  2:2284664/150‑1 (MQ=255)
    gTTCTCTACACCGACATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAAtt                                                                                                                                                 <  1:1125741/150‑1 (MQ=255)
     ttCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTg                                                                                                                                                <  1:3582501/150‑1 (MQ=255)
     ttCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTg                                                                                                                                                <  1:2468537/150‑1 (MQ=255)
        tctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTg                                                                                                                                             <  1:3343253/150‑1 (MQ=255)
        tctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTg                                                                                                                                             <  2:3582550/150‑1 (MQ=255)
          tACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGAc                                                                                                                                           >  2:2294705/1‑150 (MQ=255)
          tACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGAc                                                                                                                                           <  2:1032894/150‑1 (MQ=255)
           acacCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACa                                                                                                                                          <  2:3722014/150‑1 (MQ=255)
             acCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACaaa                                                                                                                                        >  2:2508089/1‑150 (MQ=255)
              ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                       <  1:393545/150‑1 (MQ=255)
              ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                       <  2:1468830/150‑1 (MQ=255)
              ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAAt                                                                                                                                       <  2:2450528/150‑1 (MQ=255)
               cGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATc                                                                                                                                      <  2:3953898/150‑1 (MQ=255)
                gCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCa                                                                                                                                     >  1:2098812/1‑150 (MQ=255)
                       aGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGgcag                                                                                                                              <  1:553158/150‑1 (MQ=255)
                          tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAc                                                                                                                           >  1:974822/1‑150 (MQ=255)
                               cTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcggcgg                                                                                                                      >  2:1463709/1‑150 (MQ=255)
                                  gcgcGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCa                                                                                                                   <  1:30785/150‑1 (MQ=255)
                                  gcgcGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCa                                                                                                                   <  1:1016984/150‑1 (MQ=255)
                                   cgcgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAt                                                                                                                  <  2:2025498/150‑1 (MQ=255)
                                     cgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAttt                                                                                                                <  2:766184/150‑1 (MQ=255)
                                        aTGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCg                                                                                                             >  2:2016164/1‑150 (MQ=255)
                                         tGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGc                                                                                                            >  2:3131958/1‑150 (MQ=255)
                                           ccAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCtt                                                                                                          >  2:2443206/1‑150 (MQ=255)
                                              aaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAt                                                                                                       >  1:2943274/1‑150 (MQ=255)
                                               aagaagAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATg                                                                                                      >  2:1760153/1‑150 (MQ=255)
                                                    gagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTg                                                                                                 <  2:154024/150‑1 (MQ=255)
                                                     agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGa                                                                                                >  1:810566/1‑150 (MQ=255)
                                                     agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGa                                                                                                <  1:2321120/150‑1 (MQ=255)
                                                        tGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAgg                                                                                             <  2:3363323/150‑1 (MQ=255)
                                                         ggCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGa                                                                                            <  2:3094165/150‑1 (MQ=255)
                                                           cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                          >  1:2440816/1‑150 (MQ=255)
                                                           cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                          <  2:1894941/150‑1 (MQ=255)
                                                             gaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTAGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCAATGGCAGCACGGCGGGCATTTTCGCTTGTTGGGTTGAAGGAACTg                                                                                        <  1:1206724/150‑1 (MQ=255)
                                                              aggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGACCTgg                                                                                       >  1:3886610/1‑150 (MQ=255)
                                                                  agagTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgc                                                                                   <  1:2804263/150‑1 (MQ=255)
                                                                  agagTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgc                                                                                   >  2:1978224/1‑150 (MQ=255)
                                                                    agTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgcgc                                                                                 <  1:4081548/150‑1 (MQ=255)
                                                                    agTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgcgc                                                                                 >  2:3469110/1‑150 (MQ=255)
                                                                       cgcgCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGctct                                                                              <  2:3910229/150‑1 (MQ=255)
                                                                       cgcgCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGctct                                                                              <  1:3273937/150‑1 (MQ=255)
                                                                         cgccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTAt                                                                            <  2:619174/150‑1 (MQ=255)
                                                                          gccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                           <  1:3538011/150‑1 (MQ=255)
                                                                           ccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGa                                                                          <  2:2537924/150‑1 (MQ=255)
                                                                            cgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAg                                                                         <  2:2110290/150‑1 (MQ=255)
                                                                              cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                       <  1:3461734/150‑1 (MQ=255)
                                                                               tGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                      <  2:2902993/150‑1 (MQ=255)
                                                                               tGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                      <  2:1858287/150‑1 (MQ=255)
                                                                                 ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                    <  1:1891064/150‑1 (MQ=255)
                                                                                 ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                    <  2:1763889/150‑1 (MQ=255)
                                                                                    agagATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCt                                                                 >  1:3835862/1‑150 (MQ=255)
                                                                                      agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCa                                                               <  1:2915154/150‑1 (MQ=255)
                                                                                      agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCa                                                               >  1:1225299/1‑150 (MQ=255)
                                                                                        atcatcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGt                                                             >  1:4154140/1‑150 (MQ=255)
                                                                                        atcatcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGt                                                             >  1:3427590/1‑150 (MQ=255)
                                                                                        atcatcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGt                                                             >  2:3816961/1‑150 (MQ=255)
                                                                                           atcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAt                                                          >  1:2028754/1‑150 (MQ=255)
                                                                                             cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAttt                                                        >  2:3616401/1‑150 (MQ=255)
                                                                                             cGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAAttt                                                        <  1:1758412/150‑1 (MQ=255)
                                                                                               aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                      <  2:3059254/150‑1 (MQ=255)
                                                                                                aaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaaa                                                     >  2:4102309/1‑150 (MQ=255)
                                                                                                 aTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAg                                                    <  1:1273087/150‑1 (MQ=255)
                                                                                                     cGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTgg                                                >  1:3153118/1‑150 (MQ=255)
                                                                                                         cAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                            >  1:270717/1‑150 (MQ=255)
                                                                                                             cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                        <  1:1667231/150‑1 (MQ=255)
                                                                                                             cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                        >  2:719313/1‑150 (MQ=255)
                                                                                                               tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCACTGTCGgt                                      >  2:2920113/1‑150 (MQ=255)
                                                                                                                gAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtg                                     >  2:1645974/1‑150 (MQ=255)
                                                                                                                gAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtg                                     >  2:486171/1‑150 (MQ=255)
                                                                                                                 aaCAGGCCTCCGATAACTTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTTATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTATgtgt                                    >  2:1150064/1‑150 (MQ=255)
                                                                                                                    aGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                 <  2:2408/150‑1 (MQ=255)
                                                                                                                    aGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                 <  2:2176596/150‑1 (MQ=255)
                                                                                                                     ggCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttctt                                >  1:927208/1‑150 (MQ=255)
                                                                                                                      gCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttcttc                               <  1:101906/150‑1 (MQ=255)
                                                                                                                      gCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttcttc                               >  1:1582507/1‑150 (MQ=255)
                                                                                                                       ccTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                              >  1:2166866/1‑150 (MQ=255)
                                                                                                                        cTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTtctc                             >  1:1392268/1‑150 (MQ=255)
                                                                                                                        cTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTtctc                             <  2:2653830/150‑1 (MQ=255)
                                                                                                                           cGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                          >  1:2494013/1‑150 (MQ=255)
                                                                                                                            gATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                         <  1:1103448/150‑1 (MQ=255)
                                                                                                                            gATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                         <  1:2755492/150‑1 (MQ=255)
                                                                                                                             atatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                        >  1:3104804/1‑150 (MQ=255)
                                                                                                                             atatCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                        >  1:2572462/1‑150 (MQ=255)
                                                                                                                              tatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGa                       >  2:4004559/1‑150 (MQ=255)
                                                                                                                               atCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGAt                      >  1:3546837/1‑150 (MQ=255)
                                                                                                                                  gtcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                   >  1:3167440/1‑150 (MQ=255)
                                                                                                                                   tcgAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCa                  <  2:1839736/150‑1 (MQ=255)
                                                                                                                                      aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCa               >  2:2612123/1‑150 (MQ=255)
                                                                                                                                       gcgcATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAg              <  2:3305031/150‑1 (MQ=255)
                                                                                                                                       gcgcATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAg              <  2:1405776/150‑1 (MQ=255)
                                                                                                                                       gcgcATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAg              <  1:1828830/150‑1 (MQ=255)
                                                                                                                                         gcATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgcg            >  2:3579851/1‑150 (MQ=255)
                                                                                                                                               gcagcaAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc      <  2:1229533/150‑1 (MQ=255)
                                                                                                                                               gCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc      <  2:1045620/150‑1 (MQ=255)
                                                                                                                                                  gCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGtt   <  2:1514077/150‑1 (MQ=255)
                                                                                                                                                   cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt  >  2:3206926/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTT  >  CP009273/4218780‑4219076

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: