Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,218,928 G→A 32.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,9280GA32.9% 439.8 / 299.5 301E424K (GAA→AAA) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (99/103);  new base A (37/62);  total (136/165)
Fisher's exact test for biased strand distribution p-value = 6.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

TCAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTG  >  CP009273/4218779‑4219077
                                                                                                                                                     |                                                                                                                                                     
tCAACGGTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCg                                                                                                                                                       <  2:3498214/150‑1 (MQ=255)
 cAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGa                                                                                                                                                      >  2:1913965/1‑150 (MQ=255)
 cAACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGa                                                                                                                                                      >  2:4299184/1‑150 (MQ=255)
        ctctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCt                                                                                                                                               >  2:3335950/1‑150 (MQ=255)
        ctctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCt                                                                                                                                               >  2:2084243/1‑150 (MQ=255)
        ctctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCt                                                                                                                                               <  1:1605780/150‑1 (MQ=255)
         tctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTg                                                                                                                                              <  1:4172858/150‑1 (MQ=255)
          ctACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGa                                                                                                                                             <  2:734436/150‑1 (MQ=255)
           tACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCACATGGGCAGCGAAATTGCTGAc                                                                                                                                            <  2:1510848/150‑1 (MQ=255)
              acCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGAAGAGATCATCGATATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGTGCAGCAAAATTGCTGACatt                                                                                                                                         <  2:794997/150‑3 (MQ=255)
               ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                        <  2:3943591/150‑1 (MQ=255)
               ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                        <  1:2291353/150‑1 (MQ=255)
                  ccATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCAc                                                                                                                                     >  1:1344658/1‑150 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    <  1:4344146/150‑1 (MQ=255)
                   cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACt                                                                                                                                    >  1:3221776/1‑150 (MQ=255)
                    aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTg                                                                                                                                   >  2:2189184/1‑150 (MQ=255)
                    aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTg                                                                                                                                   >  1:1531160/1‑150 (MQ=255)
                     ttAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTgg                                                                                                                                  >  2:2875155/1‑150 (MQ=255)
                        aGCTCTATCAGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGACAGAGCTCAACGAAATGTCGCACAACCATGAACAGGCCTCCGAGATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgcag                                                                                                                               <  1:3057441/150‑1 (MQ=255)
                           tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAc                                                                                                                            >  1:219137/1‑150 (MQ=255)
                            gTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcg                                                                                                                           <  1:704670/150‑1 (MQ=255)
                              aTCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcggc                                                                                                                         <  2:3617311/150‑1 (MQ=255)
                              aTCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcggc                                                                                                                         <  1:2928840/150‑1 (MQ=255)
                              aTCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcggc                                                                                                                         >  1:4192464/1‑150 (MQ=255)
                              aTCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcggc                                                                                                                         >  1:1740093/1‑150 (MQ=255)
                                  ggCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGc                                                                                                                     <  1:652787/150‑1 (MQ=255)
                                        gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  1:3944947/150‑1 (MQ=255)
                                        gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               >  1:4585783/1‑150 (MQ=255)
                                        gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  1:786114/150‑1 (MQ=255)
                                         aTGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCg                                                                                                              <  1:3842561/150‑1 (MQ=255)
                                                 agaagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATgg                                                                                                      >  2:2701342/1‑150 (MQ=255)
                                                  gaagaGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATggg                                                                                                     <  2:2157582/150‑1 (MQ=255)
                                                    agagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGtt                                                                                                   <  1:1997437/150‑1 (MQ=255)
                                                      agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGa                                                                                                 >  2:2780276/1‑150 (MQ=255)
                                                      agCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGa                                                                                                 <  1:3182643/150‑1 (MQ=255)
                                                       gCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGaa                                                                                                <  2:2799476/150‑1 (MQ=255)
                                                       gCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGaa                                                                                                <  2:1682719/150‑1 (MQ=255)
                                                          ggCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATAGAAAAGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTAGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGa                                                                                             >  1:238519/1‑150 (MQ=255)
                                                           gCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGaa                                                                                            <  2:2671308/150‑1 (MQ=255)
                                                            cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                           <  1:2872875/150‑1 (MQ=255)
                                                            cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                           >  2:2278522/1‑150 (MQ=255)
                                                            cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                           >  2:823306/1‑150 (MQ=255)
                                                             ggaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACt                                                                                          >  1:1022765/1‑150 (MQ=255)
                                                              gaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTg                                                                                         >  2:4388222/1‑150 (MQ=255)
                                                               aggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTgg                                                                                        <  2:4205843/150‑1 (MQ=255)
                                                               aggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTgg                                                                                        <  2:4504005/150‑1 (MQ=255)
                                                                  aaGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATg                                                                                     <  2:1223349/150‑1 (MQ=255)
                                                                  aaGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATg                                                                                     >  1:80232/1‑150 (MQ=255)
                                                                     agTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgcgc                                                                                  >  1:1083534/1‑150 (MQ=255)
                                                                     agTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgcgc                                                                                  >  1:2723017/1‑150 (MQ=255)
                                                                             cgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAg                                                                          >  2:1755697/1‑150 (MQ=255)
                                                                               cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTTAACAGTCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                        >  1:4774956/1‑150 (MQ=255)
                                                                               cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                        >  2:1084356/1‑150 (MQ=255)
                                                                                 gggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAGATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAt                                                                      <  2:2742666/150‑1 (MQ=255)
                                                                                  ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                     >  2:1122337/1‑150 (MQ=255)
                                                                                    cAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                   <  1:992798/150‑1 (MQ=255)
                                                                                       agaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCa                                                                >  2:1423305/1‑150 (MQ=255)
                                                                                        gaTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAg                                                               <  1:1252261/150‑1 (MQ=255)
                                                                                          tcatcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTa                                                             <  1:3295917/150‑1 (MQ=255)
                                                                                           catcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTaa                                                            >  1:3724314/1‑150 (MQ=255)
                                                                                           catcGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTaa                                                            <  1:4126282/150‑1 (MQ=255)
                                                                                               gAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTa                                                        <  2:3656721/150‑1 (MQ=255)
                                                                                                aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                       <  2:4019868/150‑1 (MQ=255)
                                                                                                   tGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGc                                                    <  1:3715267/150‑1 (MQ=255)
                                                                                                    gTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACGGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGATCAATTGCTCAGTAATTTAAAGCt                                                   <  2:2619190/150‑1 (MQ=255)
                                                                                                     tCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTg                                                  <  1:2489796/150‑1 (MQ=255)
                                                                                                      cGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTgg                                                 <  1:2939867/150‑1 (MQ=255)
                                                                                                        cTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCa                                               <  2:3016735/150‑1 (MQ=255)
                                                                                                         tCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCaa                                              >  1:3756403/1‑150 (MQ=255)
                                                                                                          cAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAAt                                             >  1:2997948/1‑150 (MQ=255)
                                                                                                                tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgt                                       >  1:2010616/1‑150 (MQ=255)
                                                                                                                 gAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtg                                      <  1:1868063/150‑1 (MQ=255)
                                                                                                                  aaCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtgt                                     <  2:2997948/150‑1 (MQ=255)
                                                                                                                  aaCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgtgt                                     <  2:4191913/150‑1 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  2:4669819/1‑150 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  1:3142213/1‑150 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  1:1521873/1‑150 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  2:1067240/1‑150 (MQ=255)
                                                                                                                   aCAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGtt                                    >  2:3191414/1‑150 (MQ=255)
                                                                                                                     aGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttct                                  >  1:930528/1‑150 (MQ=255)
                                                                                                                        ccTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                               <  1:2681556/150‑1 (MQ=255)
                                                                                                                         cTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTtctc                              >  1:3409818/1‑150 (MQ=255)
                                                                                                                            cGATCTCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTAGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                           <  2:238519/150‑1 (MQ=255)
                                                                                                                            cGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                           <  2:577346/150‑1 (MQ=255)
                                                                                                                            cGAGATCGTCGAGCGCATGGGCAGCAACATTGCTGACAAATCACTGGCAGCACGGCGGGCATGTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                           <  2:2475518/150‑1 (MQ=255)
                                                                                                                             gATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                          <  1:2918730/150‑1 (MQ=255)
                                                                                                                             gATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                          >  1:2156160/1‑150 (MQ=255)
                                                                                                                              atatCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                         >  2:2824481/1‑150 (MQ=255)
                                                                                                                                  cgtcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGt                     >  1:1679864/1‑150 (MQ=255)
                                                                                                                                      gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAcc                 >  1:821956/1‑150 (MQ=255)
                                                                                                                                      gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAcc                 >  2:1628371/1‑150 (MQ=255)
                                                                                                                                       aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCa                >  2:47935/1‑150 (MQ=255)
                                                                                                                                         cgcATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgc              >  2:1486957/1‑150 (MQ=255)
                                                                                                                                         cgcATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgc              >  1:3914895/1‑150 (MQ=255)
                                                                                                                                             tGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtc          >  1:4383959/1‑150 (MQ=255)
                                                                                                                                              gggCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg         <  2:3850407/150‑1 (MQ=255)
                                                                                                                                                gcagcaAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcgtc       <  1:4710103/150‑1 (MQ=255)
                                                                                                                                                  aGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTcgtcgt     >  2:1358146/1‑150 (MQ=255)
                                                                                                                                                   gCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGtt    <  1:490710/150‑1 (MQ=255)
                                                                                                                                                     gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGTGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTg  >  1:3519878/1‑150 (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: