Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,218,928 G→A 22.3% 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,9280GA22.3% 814.9 / 243.2 350E424K (GAA→AAA) yjbBputative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (130/142);  new base A (37/41);  total (167/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 357 total aligned reads displayed.

AACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTG  >  CP009273/4218781‑4219077
                                                                                                                                                   |                                                                                                                                                     
aaCGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCaaa                                                                                                                                                     >  2:714714/1‑150 (MQ=255)
  cGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAAt                                                                                                                                                   >  2:667755/1‑150 (MQ=255)
            acCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACaaa                                                                                                                                         <  2:4906697/150‑1 (MQ=255)
             ccGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAAt                                                                                                                                        <  1:4134044/150‑1 (MQ=255)
              cGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATc                                                                                                                                       >  1:4185977/1‑150 (MQ=255)
                ccATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCAc                                                                                                                                     >  2:3671997/1‑150 (MQ=255)
                 cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    >  1:4728492/1‑150 (MQ=255)
                 cATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCAGGGGCAGCGAAATTGCTGACAAATCACt                                                                                                                                    >  2:2757733/1‑150 (MQ=255)
                  aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTg                                                                                                                                   <  2:3680811/150‑1 (MQ=255)
                  aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTg                                                                                                                                   <  2:2966608/150‑1 (MQ=255)
                  aTTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTg                                                                                                                                   >  2:1957489/1‑150 (MQ=255)
                     aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                >  2:1568628/1‑150 (MQ=255)
                     aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                <  2:1058237/150‑1 (MQ=255)
                     aaGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgca                                                                                                                                >  1:1991185/1‑150 (MQ=255)
                        cTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGgcagca                                                                                                                             >  2:1562824/1‑150 (MQ=255)
                         tGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAc                                                                                                                            >  2:2889155/1‑150 (MQ=255)
                          gTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcg                                                                                                                           <  2:2106611/150‑1 (MQ=255)
                          gTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcg                                                                                                                           <  1:4331639/150‑1 (MQ=255)
                           tATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCAcgg                                                                                                                          <  1:147069/150‑1 (MQ=255)
                             tCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcggcg                                                                                                                        <  1:2995638/150‑1 (MQ=255)
                              cTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCAcggcgg                                                                                                                       >  1:1919537/1‑150 (MQ=255)
                               tGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCggg                                                                                                                      >  2:2436204/1‑150 (MQ=255)
                                 gcgcGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCa                                                                                                                    <  2:126411/150‑1 (MQ=255)
                                  cgcgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAt                                                                                                                   >  2:2289218/1‑150 (MQ=255)
                                   gcgGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAtt                                                                                                                  <  1:5118753/150‑1 (MQ=255)
                                     ggATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCAtttt                                                                                                                >  1:3174503/1‑150 (MQ=255)
                                      gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  2:2716553/150‑1 (MQ=255)
                                      gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               >  1:821557/1‑150 (MQ=255)
                                      gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  2:74897/150‑1 (MQ=255)
                                      gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  1:1530524/150‑1 (MQ=255)
                                      gATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTc                                                                                                               <  2:75035/150‑1 (MQ=255)
                                       aTGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCg                                                                                                              >  1:1820178/1‑150 (MQ=255)
                                        tGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGc                                                                                                             <  1:1220774/150‑1 (MQ=255)
                                          ccAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCtt                                                                                                           >  2:3879294/1‑150 (MQ=255)
                                             aaaGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAt                                                                                                        <  2:4180951/150‑1 (MQ=255)
                                              aagaagAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATg                                                                                                       <  1:3500368/150‑1 (MQ=255)
                                              aagaagAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATg                                                                                                       <  2:51228/150‑1 (MQ=255)
                                                   gagCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGCTCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTg                                                                                                  <  2:4287949/150‑1 (MQ=255)
                                                       tGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAgg                                                                                              <  2:1814713/150‑1 (MQ=255)
                                                        ggCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGa                                                                                             <  1:2054790/150‑1 (MQ=255)
                                                        ggCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGa                                                                                             <  1:2314939/150‑1 (MQ=255)
                                                          cGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAAc                                                                                           >  2:3997388/1‑150 (MQ=255)
                                                           ggaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACt                                                                                          >  1:5073720/1‑150 (MQ=255)
                                                            gaggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTg                                                                                         >  1:4954780/1‑150 (MQ=255)
                                                              ggaAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGa                                                                                       <  2:4238530/150‑1 (MQ=255)
                                                                 agagTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATgc                                                                                    >  2:3799163/1‑150 (MQ=255)
                                                                     tCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGctc                                                                                >  2:4106638/1‑150 (MQ=255)
                                                                     tCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGctc                                                                                <  1:2243887/150‑1 (MQ=255)
                                                                      cgcgCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGAGATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCCTTTTCGCTTGATGGGGTGAAGGACCTGGATGCGctct                                                                               >  2:4475055/1‑150 (MQ=255)
                                                                       gcgcCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTa                                                                              <  1:4219759/150‑1 (MQ=255)
                                                                         gccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                            <  2:4579604/150‑1 (MQ=255)
                                                                         gccgcTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATg                                                                            <  2:4856898/150‑1 (MQ=255)
                                                                             cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                        <  2:3433736/150‑1 (MQ=255)
                                                                             cTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCa                                                                        >  1:1637198/1‑150 (MQ=255)
                                                                              tGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                       >  2:3506670/1‑150 (MQ=255)
                                                                              tGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCaa                                                                       >  2:457657/1‑150 (MQ=255)
                                                                               gggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAt                                                                      <  2:5048684/150‑1 (MQ=255)
                                                                               gggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAt                                                                      <  1:3103020/150‑1 (MQ=255)
                                                                                ggCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAtt                                                                     >  1:3384425/1‑150 (MQ=255)
                                                                                  cAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGc                                                                   >  1:1494129/1‑150 (MQ=255)
                                                                                   agagATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCt                                                                  >  1:2149480/1‑150 (MQ=255)
                                                                                              aaaTGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTaa                                                       >  2:3556273/1‑150 (MQ=255)
                                                                                                  gTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCt                                                   >  2:86908/1‑150 (MQ=255)
                                                                                                   tCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTg                                                  >  1:1318193/1‑150 (MQ=255)
                                                                                                    cGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTgg                                                 >  1:4476150/1‑150 (MQ=255)
                                                                                                     gCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGc                                                >  1:1358051/1‑150 (MQ=255)
                                                                                                      cTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCa                                               <  1:359649/150‑1 (MQ=255)
                                                                                                      cTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCa                                               <  1:105938/150‑1 (MQ=255)
                                                                                                      cTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCa                                               <  2:2597241/150‑1 (MQ=255)
                                                                                                        cAACCTTGAACAGGCCTCCGATATCGTAGAGCGCATGGGAATCGAAATTGCTGACAAATCACTGGCAGCACGTCGGGCATTTTCGCTTGATGTGTTGAAGGAACTGGATGCGCTCTATGAGCAAGTGCTCAGTAATTTAAAGCTGGCAAt                                             <  1:701351/150‑1 (MQ=255)
                                                                                                            cTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCg                                         <  2:4705754/150‑1 (MQ=255)
                                                                                                              tGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGgt                                       <  2:2255022/150‑1 (MQ=255)
                                                                                                                  cAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGttc                                   >  1:2468056/1‑150 (MQ=255)
                                                                                                                      ccTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTtcttct                               <  1:4907230/150‑1 (MQ=255)
                                                                                                                        tCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTctct                             >  1:905850/1‑150 (MQ=255)
                                                                                                                          cGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTgg                           <  2:1455541/150‑1 (MQ=255)
                                                                                                                           gATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGc                          >  1:2341659/1‑150 (MQ=255)
                                                                                                                            atatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                         >  2:3631446/1‑150 (MQ=255)
                                                                                                                            atatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGAGGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                         >  2:565278/1‑150 (MQ=255)
                                                                                                                            atatCGTCGAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCg                         <  1:1586824/150‑1 (MQ=255)
                                                                                                                             tatCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAAATGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGa                        >  1:3775516/1‑150 (MQ=255)
                                                                                                                              atCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGAt                       <  1:4096253/150‑1 (MQ=255)
                                                                                                                               tcgtcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATg                      <  1:4533994/150‑1 (MQ=255)
                                                                                                                                cgtcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGt                     >  2:141838/1‑150 (MQ=255)
                                                                                                                                cgtcgAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGt                     >  2:3303349/1‑150 (MQ=255)
                                                                                                                                 gtcgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                    >  1:1106349/1‑150 (MQ=255)
                                                                                                                                 gtcgAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                    <  2:4172975/150‑1 (MQ=255)
                                                                                                                                 gtcgAGCGCATGGGCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTc                    >  1:1271133/1‑150 (MQ=255)
                                                                                                                                   cgAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAc                  <  1:3814438/150‑1 (MQ=255)
                                                                                                                                    gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAcc                 >  2:1994882/1‑150 (MQ=255)
                                                                                                                                    gAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCAcc                 >  1:2900292/1‑150 (MQ=255)
                                                                                                                                     aGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCa                >  2:3273486/1‑150 (MQ=255)
                                                                                                                                       cgcATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAgc              >  2:2316640/1‑150 (MQ=255)
                                                                                                                                           tGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtc          <  1:143432/150‑1 (MQ=255)
                                                                                                                                            gggCAGCAAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCtcg         <  2:3338086/150‑1 (MQ=255)
                                                                                                                                                aGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTcgtcgt     >  1:3701402/1‑150 (MQ=255)
                                                                                                                                                  cGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGttt   <  1:2796565/150‑1 (MQ=255)
                                                                                                                                                   gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTACTTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTg  >  1:4042985/1‑150 (MQ=255)
                                                                                                                                                   gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTg  >  2:1168916/1‑150 (MQ=255)
                                                                                                                                                   gAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTg  >  1:2677311/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACGTTCTCTACACCGCCATTAAGCTGTATCTGGCGCGGATGCCAAAAGAAGAGCTGGCGGAGGAAGAGTCGCGCCGCTGGGCAGAGATCATCGAAATGTCGCTCAACCTTGAACAGGCCTCCGATATCGTCGAGCGCATGGGCAGCGAAATTGCTGACAAATCACTGGCAGCACGGCGGGCATTTTCGCTTGATGGGTTGAAGGAACTGGATGCGCTCTATGAGCAATTGCTCAGTAATTTAAAGCTGGCAATGTCGGTGTTCTTCTCTGGCGATGTCACCAGCGCTCGTCGTTTG  >  CP009273/4218781‑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: