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,4990CA97.0% 714.7 / 11.5 232R20S (CGT→AGT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (118/107);  total (123/109)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 232 total aligned reads displayed.

CCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGG  >  CP009273/4457354‑4457643
                                                                                                                                                 |                                                                                                                                                
 ccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTg                                                                                                                                               >  2:838506/1‑148 (MQ=255)
  cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGa                                                                                                                                              <  1:3287196/148‑1 (MQ=255)
   tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGAt                                                                                                                                             >  2:1885853/1‑148 (MQ=255)
    gCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTTTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTTCCTTCCCGTCTGGTTCATAGTGATc                                                                                                                                            <  1:1359809/148‑1 (MQ=255)
    gCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATc                                                                                                                                            <  1:3076027/148‑1 (MQ=255)
      tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCa                                                                                                                                          >  1:883620/1‑148 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      <  2:2931989/148‑1 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCGTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      <  1:637563/148‑1 (MQ=255)
              attaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTgg                                                                                                                                  <  1:1720507/148‑1 (MQ=255)
              attaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTgg                                                                                                                                  <  1:2474332/148‑1 (MQ=255)
                taCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcg                                                                                                                                >  1:3247405/1‑148 (MQ=255)
                taCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcg                                                                                                                                >  1:2013345/1‑148 (MQ=255)
                  ccTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcgcg                                                                                                                              <  1:1952761/148‑1 (MQ=255)
                     tCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCgcag                                                                                                                           <  2:2751746/148‑1 (MQ=255)
                      cAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCgcagc                                                                                                                          <  2:2887008/148‑1 (MQ=255)
                         gTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGa                                                                                                                       >  2:1046077/1‑148 (MQ=255)
                            aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAg                                                                                                                    <  1:1662786/148‑1 (MQ=255)
                             ggCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                   <  2:1166922/148‑1 (MQ=255)
                              gCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTg                                                                                                                  >  1:2020709/1‑148 (MQ=255)
                               cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                 <  2:2781729/148‑1 (MQ=255)
                               cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                 <  1:1228270/148‑1 (MQ=255)
                                ttCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCGTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGaa                                                                                                                <  2:2862462/148‑1 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               <  1:1103840/148‑1 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               <  1:2826437/148‑1 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTGAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               <  2:2129467/148‑1 (MQ=255)
                                    ccACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACg                                                                                                            >  2:3066896/1‑148 (MQ=255)
                                     cACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACgg                                                                                                           <  2:2296638/148‑1 (MQ=255)
                                      aCGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggt                                                                                                          >  2:1198088/1‑148 (MQ=255)
                                          cTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTc                                                                                                      <  2:691318/148‑1 (MQ=255)
                                          cTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGATTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTc                                                                                                      <  1:2804502/148‑1 (MQ=255)
                                           tGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTcc                                                                                                     <  1:1719260/148‑1 (MQ=255)
                                               gACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgcc                                                                                                 >  2:171955/1‑148 (MQ=255)
                                               gACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgcc                                                                                                 >  2:888759/1‑148 (MQ=255)
                                                aCATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgccg                                                                                                <  2:205689/148‑1 (MQ=255)
                                                aCAGTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgccg                                                                                                <  2:1250433/148‑1 (MQ=255)
                                                 cATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                               >  2:1470435/1‑148 (MQ=255)
                                                 cATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                               >  2:3277903/1‑148 (MQ=255)
                                                   tttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTcc                                                                                             >  1:3445913/1‑148 (MQ=255)
                                                    ttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCt                                                                                            >  1:546796/1‑148 (MQ=255)
                                                      cTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                          >  1:1436947/1‑148 (MQ=255)
                                                      cTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                          >  1:127116/1‑148 (MQ=255)
                                                      cTGTACTGGTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                          <  2:2461431/148‑1 (MQ=255)
                                                       tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTaa                                                                                         <  2:403703/148‑1 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  1:2223564/148‑1 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  1:3002759/148‑1 (MQ=255)
                                                         tACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGt                                                                                       <  2:235026/148‑1 (MQ=255)
                                                          aCTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                      >  1:439413/1‑148 (MQ=255)
                                                          aCTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                      <  2:614665/148‑1 (MQ=255)
                                                               ttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAt                                                                                 <  1:2811977/148‑1 (MQ=255)
                                                                  aGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGc                                                                              <  2:2613724/148‑1 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  1:2584752/1‑148 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  1:2683071/1‑148 (MQ=255)
                                                                         gCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaaa                                                                       <  2:87220/148‑1 (MQ=255)
                                                                          cGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAAt                                                                      <  2:1611324/148‑1 (MQ=255)
                                                                             aGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGc                                                                   >  2:1918824/1‑148 (MQ=255)
                                                                                gACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGt                                                                >  2:785856/1‑148 (MQ=255)
                                                                                    ccTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                            >  1:3183896/1‑148 (MQ=255)
                                                                                      ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                          <  2:2822009/148‑1 (MQ=255)
                                                                                      ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                          <  2:3510424/148‑1 (MQ=255)
                                                                                       tATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGccc                                                         <  2:168023/148‑1 (MQ=255)
                                                                                         tgtggaaaGAAATTTTTACCCCGCTCAGTCGCCATTGACCTTCCCGTCTGGTTCATAGTGCTCCATGGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATGGCCTGAAAATGGCGGTGATGCCCga                                                       <  2:518981/143‑1 (MQ=255)
                                                                                         tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                       >  2:2105628/1‑148 (MQ=255)
                                                                                         tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                       <  2:1873176/148‑1 (MQ=255)
                                                                                         tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                       <  1:344761/148‑1 (MQ=255)
                                                                                          gTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgac                                                      >  2:3167397/1‑147 (MQ=255)
                                                                                            ttAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaaga                                                    >  1:2395891/1‑148 (MQ=255)
                                                                                               aaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaaga                                                 <  1:18266/148‑1 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                <  2:566326/148‑1 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              >  1:268432/1‑148 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              >  1:2521273/1‑148 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              >  1:3332657/1‑148 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              >  1:3275717/1‑148 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              <  2:1928422/148‑1 (MQ=255)
                                                                                                   aaTTTTTACCTGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtg                                             <  1:1219850/148‑1 (MQ=255)
                                                                                                   aaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtg                                             <  1:3235210/148‑1 (MQ=255)
                                                                                                      ttttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATGGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                          <  2:1270900/148‑1 (MQ=255)
                                                                                                          aCCCGGATCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaaa                                      >  2:50790/1‑148 (MQ=255)
                                                                                                                 tCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGa                               >  2:2651050/1‑148 (MQ=255)
                                                                                                                  cATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGGCCTTAAGTGCGCATTGCCTGAAAATTGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGac                              >  2:140811/1‑148 (MQ=255)
                                                                                                                   aTTCGCCATTTACCTTCCCGTCTGGTACATAGTGATCCATTGCCTGGCGCGCATCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGTACAACGTACGaca                             >  1:3406597/1‑148 (MQ=255)
                                                                                                                    ttCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAACACGTTCGacac                            >  2:2256187/1‑148 (MQ=255)
                                                                                                                      cGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAccc                          <  2:1954437/148‑1 (MQ=255)
                                                                                                                       gCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCa                         >  1:2357868/1‑148 (MQ=255)
                                                                                                                        ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                        <  2:164651/148‑1 (MQ=255)
                                                                                                                           tttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCg                     >  1:380112/1‑148 (MQ=255)
                                                                                                                           tttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCg                     >  2:3248013/1‑148 (MQ=255)
                                                                                                                            ttACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCgg                    <  1:3209659/148‑1 (MQ=255)
                                                                                                                            ttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCgg                    >  2:821129/1‑148 (MQ=255)
                                                                                                                             tACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                   <  2:477289/148‑1 (MQ=255)
                                                                                                                              aCCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGaa                  >  1:2380728/1‑148 (MQ=255)
                                                                                                                                 ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg               <  1:2386477/148‑1 (MQ=255)
                                                                                                                                 ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg               >  1:2829743/1‑148 (MQ=255)
                                                                                                                                 ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg               <  1:620683/148‑1 (MQ=255)
                                                                                                                                  tCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGt              >  1:247731/1‑148 (MQ=255)
                                                                                                                                   cccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGACGAATTGTGGAAAACGTTCGACACCCATCCGGACGGGtt             >  2:3262614/1‑148 (MQ=255)
                                                                                                                                   cccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGtt             <  2:2557620/148‑1 (MQ=255)
                                                                                                                                    ccgacaggTTCATAGTGATCCATTGCCAGGCGCGCAGCAGACAGTGAATACGGAGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCTGTGATGACCGAAGAAGAATTGTGGAAAACGTTCGACACACATCCGGAAGGGTTa            >  2:298424/7‑148 (MQ=255)
                                                                                                                                    ccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa            <  1:2890826/148‑1 (MQ=255)
                                                                                                                                        cTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATc        >  2:853155/1‑148 (MQ=255)
                                                                                                                                              cATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCgg  >  2:853660/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGG  >  CP009273/4457354‑4457643

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: