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,4990CA96.5% 694.0 / 14.3 230R20S (CGT→AGT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (129/93);  total (131/99)
Fisher's exact test for biased strand distribution p-value = 7.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 230 total aligned reads displayed.

CCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGG  >  CP009273/4457354‑4457643
                                                                                                                                                 |                                                                                                                                                
cccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGt                                                                                                                                                >  2:1413432/1‑148 (MQ=255)
cccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGt                                                                                                                                                <  2:1223419/148‑1 (MQ=255)
cccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGt                                                                                                                                                >  1:36812/1‑148 (MQ=255)
 ccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTg                                                                                                                                               <  1:3205562/148‑1 (MQ=255)
  cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGa                                                                                                                                              >  2:307687/1‑148 (MQ=255)
  cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGa                                                                                                                                              <  2:238125/148‑1 (MQ=255)
  cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGa                                                                                                                                              >  1:3119753/1‑148 (MQ=255)
  cTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGa                                                                                                                                              <  1:2506161/148‑1 (MQ=255)
   tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGAt                                                                                                                                             <  1:208180/148‑1 (MQ=255)
   tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACACGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGAt                                                                                                                                             >  1:101137/1‑148 (MQ=255)
      tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCa                                                                                                                                          <  2:2830927/148‑1 (MQ=255)
       cAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAt                                                                                                                                         <  2:2139248/148‑1 (MQ=255)
        aGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAtt                                                                                                                                        >  1:1247021/1‑148 (MQ=255)
        aGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAtt                                                                                                                                        <  2:646722/148‑1 (MQ=255)
         gCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTg                                                                                                                                       >  2:1210906/1‑148 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      >  2:1686941/1‑148 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      >  2:918801/1‑148 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      >  2:2514736/1‑148 (MQ=255)
                taCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcg                                                                                                                                >  1:2772873/1‑148 (MQ=255)
                     tCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCgcag                                                                                                                           <  2:2995806/148‑1 (MQ=255)
                     tCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCgcag                                                                                                                           <  2:819220/148‑1 (MQ=255)
                         gTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGa                                                                                                                       <  2:726308/148‑1 (MQ=255)
                          tAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGAc                                                                                                                      >  2:2093874/1‑148 (MQ=255)
                          tAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGAc                                                                                                                      >  1:1195946/1‑148 (MQ=255)
                           aaGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACa                                                                                                                     >  1:2809796/1‑148 (MQ=255)
                            aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAg                                                                                                                    <  2:1463035/148‑1 (MQ=255)
                             ggCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                   >  2:383988/1‑148 (MQ=255)
                              gCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTg                                                                                                                  >  2:281256/1‑148 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               >  1:2519947/1‑148 (MQ=255)
                                 tCGCCACCCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               >  2:2268180/1‑148 (MQ=255)
                                  cGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATa                                                                                                              <  1:2095427/148‑1 (MQ=255)
                                   gCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATAc                                                                                                             >  1:2174965/1‑148 (MQ=255)
                                         ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtggt                                                                                                       >  1:560116/1‑148 (MQ=255)
                                         ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCGTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtggt                                                                                                       <  1:3022887/148‑1 (MQ=255)
                                            gAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTccc                                                                                                    >  1:1931591/1‑148 (MQ=255)
                                             aaGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccg                                                                                                   >  2:2040907/1‑148 (MQ=255)
                                             aaGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccg                                                                                                   >  1:3409100/1‑148 (MQ=255)
                                                 cATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                               >  2:1241714/1‑148 (MQ=255)
                                                  aTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                              >  1:1894742/1‑148 (MQ=255)
                                                   tttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTcc                                                                                             >  2:1307392/1‑148 (MQ=255)
                                                    ttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCt                                                                                            >  1:1462607/1‑148 (MQ=255)
                                                      cTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                          <  2:647039/148‑1 (MQ=255)
                                                      cTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTa                                                                                          <  2:884543/148‑1 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  2:2096973/148‑1 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        >  1:2880326/1‑148 (MQ=255)
                                                         tACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGt                                                                                       >  1:1952907/1‑148 (MQ=255)
                                                          aCTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCCTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                      >  2:1331760/1‑148 (MQ=255)
                                                            tGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTgcg                                                                                    >  2:1999199/1‑148 (MQ=255)
                                                              tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGAGCa                                                                                  >  2:1528749/1‑148 (MQ=255)
                                                                 cAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTg                                                                               <  1:3393522/148‑1 (MQ=255)
                                                                   gACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGcc                                                                             <  2:2151997/148‑1 (MQ=255)
                                                                   gACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGcc                                                                             <  2:1838670/148‑1 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  1:787318/1‑148 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  2:2340956/1‑148 (MQ=255)
                                                                            gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATgg                                                                    >  1:3337953/1‑148 (MQ=255)
                                                                            gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATgg                                                                    <  1:2606346/148‑1 (MQ=255)
                                                                              gggACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCg                                                                  <  2:566114/148‑1 (MQ=255)
                                                                                gACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGt                                                                <  1:918801/148‑1 (MQ=255)
                                                                                gACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGt                                                                >  1:1291756/1‑148 (MQ=255)
                                                                                          gTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaa                                                      >  2:2950919/1‑148 (MQ=255)
                                                                                          gTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaa                                                      >  2:314281/1‑148 (MQ=255)
                                                                                              aaaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                  <  1:1750145/148‑1 (MQ=255)
                                                                                              aaaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                  >  2:359202/1‑148 (MQ=255)
                                                                                               aaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaaga                                                 <  2:189535/148‑1 (MQ=255)
                                                                                               aaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaaga                                                 <  2:1539820/148‑1 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                <  1:1566928/148‑1 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                >  1:3085785/1‑148 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                <  1:351345/148‑1 (MQ=255)
                                                                                                 gAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                               <  2:1575562/148‑1 (MQ=255)
                                                                                                  aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                              >  1:3496447/1‑148 (MQ=255)
                                                                                                    aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtgt                                            >  2:2551915/1‑148 (MQ=255)
                                                                                                    aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtgt                                            <  2:1011153/148‑1 (MQ=255)
                                                                                                       tttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGa                                         <  2:2966095/148‑1 (MQ=255)
                                                                                                        ttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaa                                        <  1:2956815/148‑1 (MQ=255)
                                                                                                         tACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaa                                       <  1:566818/148‑1 (MQ=255)
                                                                                                          aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaaa                                      >  2:1639002/1‑148 (MQ=255)
                                                                                                          aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaaa                                      <  1:1413432/148‑1 (MQ=255)
                                                                                                           cccGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAAc                                     <  1:2550713/148‑1 (MQ=255)
                                                                                                            ccGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACg                                    <  2:243597/148‑1 (MQ=255)
                                                                                                            ccGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACg                                    >  2:2775703/1‑148 (MQ=255)
                                                                                                             cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGt                                   <  1:2942662/148‑1 (MQ=255)
                                                                                                             cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGt                                   >  1:1822054/1‑148 (MQ=255)
                                                                                                             cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGt                                   <  2:3334099/148‑1 (MQ=255)
                                                                                                               gCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTc                                 <  1:2056532/148‑1 (MQ=255)
                                                                                                                  cATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGac                              >  1:1142764/1‑148 (MQ=255)
                                                                                                                    ttCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGacac                            >  2:155732/1‑148 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           >  2:442963/1‑148 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           >  2:114944/1‑148 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           <  2:1582987/148‑1 (MQ=255)
                                                                                                                        ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                        >  1:2120353/1‑148 (MQ=255)
                                                                                                                             tACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                   >  1:3042958/1‑148 (MQ=255)
                                                                                                                                  tCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGt              <  2:1932131/148‑1 (MQ=255)
                                                                                                                                   cccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATGGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGtt             <  1:2093621/148‑1 (MQ=255)
                                                                                                                                    ccGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa            <  1:3286781/148‑1 (MQ=255)
                                                                                                                                    ccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa            <  2:2772873/148‑1 (MQ=255)
                                                                                                                                     cGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaa           <  1:234918/148‑1 (MQ=255)
                                                                                                                                      gTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaaa          <  1:2325118/148‑1 (MQ=255)
                                                                                                                                        cTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATc        >  1:2692123/1‑148 (MQ=255)
                                                                                                                                            ttCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGc    >  1:3259219/1‑148 (MQ=255)
                                                                                                                                              cATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCgg  >  1:1022244/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: ‑