Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,457,499 C→A 100% R20S (CGT→AGT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,457,4990CA99.3% 926.4 / ‑2.1 282R20S (CGT→AGT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (135/145);  total (135/147)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 282 total aligned reads displayed.

GTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAA  >  CP009273/4457352‑4457645
                                                                                                                                                   |                                                                                                                                                  
 tCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTg                                                                                                                                                 >  1:3814588/1‑150 (MQ=255)
    cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATc                                                                                                                                              >  1:4271256/1‑150 (MQ=255)
    cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATc                                                                                                                                              >  1:2419680/1‑150 (MQ=255)
     tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATcc                                                                                                                                             <  1:723838/150‑1 (MQ=255)
      gCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCa                                                                                                                                            >  2:4896953/1‑150 (MQ=255)
       cTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAt                                                                                                                                           >  1:1606959/1‑150 (MQ=255)
        tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAtt                                                                                                                                          <  2:3372585/150‑1 (MQ=255)
         cAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTg                                                                                                                                         >  2:586555/1‑150 (MQ=255)
         cAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTg                                                                                                                                         <  2:5124412/150‑1 (MQ=255)
           gCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGcc                                                                                                                                       <  2:442082/150‑1 (MQ=255)
            cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCt                                                                                                                                      <  1:3433648/150‑1 (MQ=255)
             tttATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTg                                                                                                                                     <  2:5081555/150‑1 (MQ=255)
             tttATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTg                                                                                                                                     <  2:2524356/150‑1 (MQ=255)
               tattaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgc                                                                                                                                   <  1:3746755/150‑1 (MQ=255)
                 ttaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcgc                                                                                                                                 <  2:4965636/150‑1 (MQ=255)
                        cAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGcagcag                                                                                                                          >  2:3324361/1‑150 (MQ=255)
                           gTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACa                                                                                                                       <  1:2221640/150‑1 (MQ=255)
                            tAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAg                                                                                                                      <  1:575675/150‑1 (MQ=255)
                               ggCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGTGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                   <  1:1134322/150‑1 (MQ=255)
                               ggCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                   >  2:2626732/1‑150 (MQ=255)
                               ggCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                   <  1:4522003/150‑1 (MQ=255)
                                 cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                                 >  2:693569/1‑150 (MQ=255)
                                  ttCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATa                                                                                                                >  1:739773/1‑150 (MQ=255)
                                   tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATAc                                                                                                               >  2:3632452/1‑150 (MQ=255)
                                   tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATAc                                                                                                               >  2:804160/1‑150 (MQ=255)
                                    cGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACg                                                                                                              >  1:4585300/1‑150 (MQ=255)
                                    cGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACg                                                                                                              <  1:4197464/150‑1 (MQ=255)
                                      ccACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggt                                                                                                            >  2:4476590/1‑150 (MQ=255)
                                       cACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtg                                                                                                           >  1:1394715/1‑150 (MQ=255)
                                            cTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTccc                                                                                                      >  2:3434626/1‑150 (MQ=255)
                                              gAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGAAGACAGTGAATACGGTGGTCccgc                                                                                                    <  1:4062477/150‑1 (MQ=255)
                                               aaGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgcc                                                                                                   >  1:5080105/1‑150 (MQ=255)
                                                 gACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                                 <  2:900590/150‑1 (MQ=255)
                                                   cATTTATTTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTcc                                                                                               <  2:297579/150‑1 (MQ=255)
                                                    aTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCt                                                                                              <  1:4284973/150‑1 (MQ=255)
                                                    aTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCt                                                                                              <  2:4633153/150‑1 (MQ=255)
                                                      ttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                            <  1:4962781/150‑1 (MQ=255)
                                                       tCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTaa                                                                                           <  1:4614875/150‑1 (MQ=255)
                                                          gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                        >  2:3866282/1‑150 (MQ=255)
                                                          gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                        >  1:2841214/1‑150 (MQ=255)
                                                              tGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                    >  2:3672343/1‑150 (MQ=255)
                                                              tGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                    >  2:2773771/1‑150 (MQ=255)
                                                               gTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAt                                                                                   <  1:1840404/150‑1 (MQ=255)
                                                                tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAtt                                                                                  <  1:2145141/150‑1 (MQ=255)
                                                                 ttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTg                                                                                 <  2:1394715/150‑1 (MQ=255)
                                                                   cAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGcc                                                                               >  2:4138690/1‑150 (MQ=255)
                                                                      aCAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGa                                                                            <  2:1547284/150‑1 (MQ=255)
                                                                         gTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaaa                                                                         >  2:1655515/1‑150 (MQ=255)
                                                                             ggAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGc                                                                     >  2:1793687/1‑150 (MQ=255)
                                                                             ggAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGc                                                                     <  1:4281299/150‑1 (MQ=255)
                                                                              gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCg                                                                    <  1:3270854/150‑1 (MQ=255)
                                                                              gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCg                                                                    >  2:3066007/1‑150 (MQ=255)
                                                                              gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCg                                                                    >  2:4248602/1‑150 (MQ=255)
                                                                                 ggACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTg                                                                 <  1:758940/150‑1 (MQ=255)
                                                                                  gACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGa                                                                <  1:917841/150‑1 (MQ=255)
                                                                                   aCTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGAt                                                               <  1:2552070/150‑1 (MQ=255)
                                                                                    cTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                              >  2:3085768/1‑150 (MQ=255)
                                                                                     tCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCATGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCATGAAGTGCGCATTGACTGAAAATGGCGGTGATGc                                                             >  1:1318103/1‑150 (MQ=255)
                                                                                      ccTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                            <  1:1525017/150‑1 (MQ=255)
                                                                                         tATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                         <  2:3339719/150‑1 (MQ=255)
                                                                                         tATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                         <  1:3286828/150‑1 (MQ=255)
                                                                                           tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaag                                                       <  2:3814588/150‑1 (MQ=255)
                                                                                           tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaag                                                       <  2:3161850/150‑1 (MQ=255)
                                                                                              ttAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGAGCTGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                    >  2:1332609/1‑150 (MQ=255)
                                                                                                      aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                            >  2:1170344/1‑150 (MQ=255)
                                                                                                       tttttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGa                                           >  1:303867/1‑150 (MQ=255)
                                                                                                            aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACg                                      <  2:2475906/150‑1 (MQ=255)
                                                                                                            aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACg                                      >  1:2511895/1‑150 (MQ=255)
                                                                                                            aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACg                                      <  2:4789764/150‑1 (MQ=255)
                                                                                                             cccGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGt                                     >  2:4171336/1‑150 (MQ=255)
                                                                                                              ccGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGtt                                    >  2:3894931/1‑150 (MQ=255)
                                                                                                               cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTc                                   <  1:3999010/150‑1 (MQ=255)
                                                                                                                ggCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                  <  2:780565/150‑1 (MQ=255)
                                                                                                                 gCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGa                                 <  1:4896953/150‑1 (MQ=255)
                                                                                                                 gCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGa                                 <  1:2336921/150‑1 (MQ=255)
                                                                                                                  cTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGac                                <  1:693569/150‑1 (MQ=255)
                                                                                                                   tCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGaca                               <  1:1940076/150‑1 (MQ=255)
                                                                                                                    cATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGacac                              >  2:1159693/1‑150 (MQ=255)
                                                                                                                     aTTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                             >  1:476919/1‑150 (MQ=255)
                                                                                                                      ttCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAccc                            >  1:1097221/1‑150 (MQ=255)
                                                                                                                       tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCa                           >  1:793820/1‑150 (MQ=255)
                                                                                                                        cGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                          <  2:884940/150‑1 (MQ=255)
                                                                                                                        cGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                          >  1:357374/1‑150 (MQ=255)
                                                                                                                         gCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATc                         <  2:2087216/150‑1 (MQ=255)
                                                                                                                          ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATcc                        >  1:16926/1‑150 (MQ=255)
                                                                                                                          ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATcc                        <  1:4622520/150‑1 (MQ=255)
                                                                                                                           cATGTACCTTCCCGGCTGGTTCAAAGTGATACATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGGGCCTGAAAATGGCGGAGATGCCCGTAGAAGAATTGTGGAAAACGTTCGACACCCATcag                       <  2:2312800/150‑3 (MQ=255)
                                                                                                                            aTTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCgg                      >  2:4066835/1‑150 (MQ=255)
                                                                                                                             tttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                     <  1:586555/150‑1 (MQ=255)
                                                                                                                                  cTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGt                >  2:3577528/1‑150 (MQ=255)
                                                                                                                                   ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGtt               <  1:1293444/150‑1 (MQ=255)
                                                                                                                                    tCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa              >  1:4378855/1‑150 (MQ=255)
                                                                                                                                    tCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa              >  1:4175968/1‑150 (MQ=255)
                                                                                                                                     cccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaa             >  1:3346730/1‑150 (MQ=255)
                                                                                                                                          cTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAg        >  1:2542549/1‑150 (MQ=255)
                                                                                                                                           tGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg       <  1:2687856/150‑1 (MQ=255)
                                                                                                                                            ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGc      <  2:737397/150‑1 (MQ=255)
                                                                                                                                              ttCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCgg    >  1:2204716/1‑150 (MQ=255)
                                                                                                                                              ttCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCgg    <  2:266192/150‑1 (MQ=255)
                                                                                                                                                cATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGaa  <  2:779364/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGAA  >  CP009273/4457352‑4457645

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: