Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,458,609 G→T 93.7% D390Y (GAT→TAT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,6090GT93.7% 920.3 / 48.9 301D390Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base T (124/158);  total (132/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 301 total aligned reads displayed.

GCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAA  >  CP009273/4458461‑4458755
                                                                                                                                                    |                                                                                                                                                  
 cTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTat                                                                                                                                                  <  2:708136/150‑1 (MQ=255)
   ggCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGTTAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTatat                                                                                                                                                >  1:2884301/1‑150 (MQ=255)
   ggCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTatat                                                                                                                                                <  2:3529811/150‑1 (MQ=255)
       cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACAGCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTcc                                                                                                                                            >  1:247746/1‑150 (MQ=255)
          ggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGt                                                                                                                                         >  2:3875235/1‑150 (MQ=255)
          ggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGt                                                                                                                                         >  2:3875068/1‑150 (MQ=255)
               aGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                    <  1:3492115/150‑1 (MQ=255)
                 tAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                  <  2:1443566/150‑1 (MQ=255)
                 tAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                  <  2:4435839/150‑1 (MQ=255)
                 tAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                  <  2:3721068/150‑1 (MQ=255)
                  aaaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAg                                                                                                                                 >  1:3718268/1‑150 (MQ=255)
                    aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAGTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCg                                                                                                                               >  1:184706/1‑150 (MQ=255)
                     gCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGa                                                                                                                              >  1:4332022/1‑150 (MQ=255)
                        gTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcg                                                                                                                           <  2:2320727/150‑1 (MQ=255)
                          cGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcgcg                                                                                                                         >  1:1713880/1‑150 (MQ=255)
                          cGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcgcg                                                                                                                         >  1:3170097/1‑150 (MQ=255)
                             aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                      >  1:3319777/1‑150 (MQ=255)
                             aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                      >  2:1504148/1‑150 (MQ=255)
                             aaCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTAGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgc                                                                                                                      <  2:339257/150‑1 (MQ=255)
                                    aGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAg                                                                                                               <  2:577214/150‑1 (MQ=255)
                                      tGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgcg                                                                                                             <  2:4150004/150‑1 (MQ=255)
                                        aTCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGcgcg                                                                                                           <  1:3010376/150‑1 (MQ=255)
                                         tCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGt                                                                                                          <  1:102784/150‑1 (MQ=255)
                                          cGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTg                                                                                                         >  1:4453511/1‑150 (MQ=255)
                                          cGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTg                                                                                                         >  1:1802245/1‑150 (MQ=255)
                                            tCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGt                                                                                                       <  1:3917993/150‑1 (MQ=255)
                                            tCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGt                                                                                                       <  1:3917953/150‑1 (MQ=255)
                                              aaaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTg                                                                                                     >  2:937106/1‑150 (MQ=255)
                                              aaaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTg                                                                                                     >  2:2964393/1‑150 (MQ=255)
                                               aaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                    <  2:1034136/150‑1 (MQ=255)
                                                aCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGaa                                                                                                   >  1:2485245/1‑150 (MQ=255)
                                                aCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGaa                                                                                                   >  1:1956391/1‑150 (MQ=255)
                                                     tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                              >  1:1486596/1‑150 (MQ=255)
                                                     tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                              >  1:1576563/1‑150 (MQ=255)
                                                     tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTc                                                                                              >  1:3296412/1‑150 (MQ=255)
                                                      ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCa                                                                                             >  1:1399074/1‑150 (MQ=255)
                                                      ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCa                                                                                             >  1:2658497/1‑150 (MQ=255)
                                                       gATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAt                                                                                            >  2:4013572/1‑150 (MQ=255)
                                                       gATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAt                                                                                            <  2:136601/150‑1 (MQ=255)
                                                       gATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAt                                                                                            <  2:3036148/150‑1 (MQ=255)
                                                        aTGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAtt                                                                                           <  2:4508071/150‑1 (MQ=255)
                                                          gCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTAt                                                                                         >  2:2144504/1‑150 (MQ=255)
                                                          gCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTAt                                                                                         <  1:3390930/150‑1 (MQ=255)
                                                           cTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATc                                                                                        >  2:2139381/1‑150 (MQ=255)
                                                            tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                       <  2:477718/150‑1 (MQ=255)
                                                            tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                       >  1:3372283/1‑150 (MQ=255)
                                                            tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCa                                                                                       <  2:590746/150‑1 (MQ=255)
                                                               tCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGAc                                                                                    <  2:4317336/150‑1 (MQ=255)
                                                                 aGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCg                                                                                  >  1:243625/1‑150 (MQ=255)
                                                                 aGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCg                                                                                  >  1:2041818/1‑150 (MQ=255)
                                                                      tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACtt                                                                             <  1:1065566/150‑1 (MQ=255)
                                                                      tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACtt                                                                             <  1:7992/150‑1 (MQ=255)
                                                                            gcAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                       <  1:2115859/150‑1 (MQ=255)
                                                                             cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAcc                                                                      <  1:2269505/150‑1 (MQ=255)
                                                                              taTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCt                                                                     <  1:2622377/149‑1 (MQ=255)
                                                                              aaTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCt                                                                     <  2:781521/150‑1 (MQ=255)
                                                                               aTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTg                                                                    <  1:839100/150‑1 (MQ=255)
                                                                                   atatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCg                                                                >  1:2409215/1‑150 (MQ=255)
                                                                                   atatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCg                                                                <  1:3243083/150‑1 (MQ=255)
                                                                                    tatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGa                                                               >  2:154686/1‑150 (MQ=255)
                                                                                    tatTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGa                                                               <  2:899668/150‑1 (MQ=255)
                                                                                     atTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGAt                                                              >  1:1293381/1‑150 (MQ=255)
                                                                                       tCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATAc                                                            <  2:204198/150‑1 (MQ=255)
                                                                                       tCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATAc                                                            <  2:2386515/150‑1 (MQ=255)
                                                                                       tCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATAc                                                            <  2:2171934/150‑1 (MQ=255)
                                                                                           tgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCg                                                        <  2:2506590/150‑1 (MQ=255)
                                                                                              aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTg                                                     >  2:1941702/1‑150 (MQ=255)
                                                                                              aCTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTg                                                     <  2:1989673/150‑1 (MQ=255)
                                                                                               cTGATACAACCGGCACCCTGACGCAGGATAATATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACATTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                    <  1:400235/150‑1 (MQ=255)
                                                                                               cTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                    >  1:2931457/1‑150 (MQ=255)
                                                                                               cTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                    <  2:3778927/150‑1 (MQ=255)
                                                                                               cTGAAAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                    <  1:400165/150‑1 (MQ=255)
                                                                                                       aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGt                                            <  1:4308071/150‑1 (MQ=255)
                                                                                                       aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGt                                            <  1:2206065/150‑1 (MQ=255)
                                                                                                          ggCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACg                                         >  2:2171477/1‑150 (MQ=255)
                                                                                                          ggCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACg                                         >  1:1648474/1‑150 (MQ=255)
                                                                                                           gCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACgg                                        >  2:3613280/1‑150 (MQ=255)
                                                                                                            cACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                       <  1:3584295/150‑1 (MQ=255)
                                                                                                             aCCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGAt                                      <  1:1359152/150‑1 (MQ=255)
                                                                                                             aCCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGAt                                      <  2:1692896/150‑1 (MQ=255)
                                                                                                              cccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTGCGGATg                                     >  1:2548378/1‑150 (MQ=255)
                                                                                                                 tGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAag                                  >  2:4278568/1‑150 (MQ=255)
                                                                                                                   aCGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAagag                                <  1:937106/150‑1 (MQ=255)
                                                                                                                           tAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgcgcg                        <  1:1116618/150‑1 (MQ=255)
                                                                                                                                 tgtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgct                  <  1:3875068/150‑1 (MQ=255)
                                                                                                                                 tgtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgct                  <  1:3875235/150‑1 (MQ=255)
                                                                                                                                 tgtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctcgct                  <  2:126067/150‑1 (MQ=255)
                                                                                                                                  gtgCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTg                 <  2:3299346/150‑1 (MQ=255)
                                                                                                                                    gCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGc               <  2:1188216/150‑1 (MQ=255)
                                                                                                                                     cTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGcc              >  2:1611000/1‑150 (MQ=255)
                                                                                                                                      tGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa             >  2:2584162/1‑150 (MQ=255)
                                                                                                                                      tGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCa             >  2:2583338/1‑150 (MQ=255)
                                                                                                                                       ggAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg            >  1:632959/1‑150 (MQ=255)
                                                                                                                                        gagaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGt           <  1:4360375/150‑1 (MQ=255)
                                                                                                                                          gaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg         <  2:1761085/150‑1 (MQ=255)
                                                                                                                                          gaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTGAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg         >  2:653626/1‑150 (MQ=255)
                                                                                                                                              cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGATCAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg     <  1:4225186/150‑1 (MQ=255)
                                                                                                                                              cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg     >  1:71197/1‑150 (MQ=255)
                                                                                                                                               aTACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc    <  2:3170097/150‑1 (MQ=255)
                                                                                                                                                 aCCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCaa  >  2:2797603/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTGGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGCAA  >  CP009273/4458461‑4458755

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: ‑