Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,483,772 C→A 7.7% A139A (GCG→GCT pabA ← aminodeoxychorismate synthase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,483,7720CA 7.7% 1028.5 / 36.4 273A139A (GCG→GCTpabAaminodeoxychorismate synthase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (133/119);  new base A (14/7);  total (147/126)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

TCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACTTTACC  >  CP009273/3483625‑3483918
                                                                                                                                                   |                                                                                                                                                  
tCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACg                                                                                                                                                   >  1:4752958/1‑149 (MQ=255)
  aGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCt                                                                                                                                                 >  2:4319767/1‑149 (MQ=255)
  aGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCt                                                                                                                                                 >  2:1993587/1‑149 (MQ=255)
   gCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTg                                                                                                                                                >  1:1045707/1‑149 (MQ=255)
   gCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTg                                                                                                                                                >  1:2377694/1‑149 (MQ=255)
    cGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTgg                                                                                                                                               <  2:4449885/149‑1 (MQ=255)
      aTGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTa                                                                                                                                             >  2:131546/1‑149 (MQ=255)
      aTGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTa                                                                                                                                             >  2:131669/1‑149 (MQ=255)
         cAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATg                                                                                                                                          <  1:4468201/149‑1 (MQ=255)
           ggAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAg                                                                                                                                        <  1:1651105/149‑1 (MQ=255)
            gAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAAAACTTTCTGGATGGAACTGCACACCTTCAAGATCCCACTGGCGATGGCGAATCCCCATAATCACTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGt                                                                                                                                       >  2:2244975/1‑149 (MQ=255)
                  aGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGtt                                                                                                                                 >  2:3249708/1‑149 (MQ=255)
                  aGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTGTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGtt                                                                                                                                 <  2:3544385/149‑1 (MQ=255)
                    ccAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTcc                                                                                                                               <  2:3875360/149‑1 (MQ=255)
                     cagcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTcca                                                                                                                              >  2:2631695/1‑149 (MQ=255)
                     cagcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTcca                                                                                                                              >  2:3783154/1‑149 (MQ=255)
                       gcagTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacc                                                                                                                            >  2:3700955/1‑149 (MQ=255)
                          gTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacc                                                                                                                         >  2:3764412/1‑149 (MQ=255)
                           ttGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacca                                                                                                                        <  1:3932643/149‑1 (MQ=255)
                           ttGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacca                                                                                                                        <  2:32578/149‑1 (MQ=255)
                           ttGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccaccacca                                                                                                                        <  2:3773201/149‑1 (MQ=255)
                            tGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAg                                                                                                                       >  1:935152/1‑149 (MQ=255)
                              aTGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCg                                                                                                                     <  2:2312416/149‑1 (MQ=255)
                              aTGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCg                                                                                                                     >  1:439718/1‑149 (MQ=255)
                               tGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGa                                                                                                                    >  1:2946710/1‑149 (MQ=255)
                                gTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGaa                                                                                                                   <  2:3336093/149‑1 (MQ=255)
                                 tCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAAt                                                                                                                  >  1:3421621/1‑149 (MQ=255)
                                     tGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTa                                                                                                              <  1:1693621/149‑1 (MQ=255)
                                      gTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAg                                                                                                             <  1:3838604/149‑1 (MQ=255)
                                       ttCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGc                                                                                                            >  2:3569785/1‑149 (MQ=255)
                                          gCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCgtg                                                                                                         >  1:736950/1‑149 (MQ=255)
                                           cTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTCATGAGTCAGGTTCCACCACCAGCGAATGGTAGCgtgt                                                                                                        >  2:1861356/1‑149 (MQ=255)
                                           cTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCgtgt                                                                                                        >  2:2538829/1‑149 (MQ=255)
                                            tAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTc                                                                                                       <  2:734563/149‑1 (MQ=255)
                                             aaGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCa                                                                                                      >  2:2066336/1‑149 (MQ=255)
                                             aaGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCa                                                                                                      >  2:1202055/1‑149 (MQ=255)
                                              aGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCAc                                                                                                     <  1:4559190/149‑1 (MQ=255)
                                              aGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCAc                                                                                                     >  1:3637955/1‑149 (MQ=255)
                                               gAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACg                                                                                                    >  2:2580857/1‑149 (MQ=255)
                                               gAATACTGTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACg                                                                                                    <  2:4224376/149‑1 (MQ=255)
                                                 aTACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGt                                                                                                  <  2:3384879/149‑1 (MQ=255)
                                                   aCTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTaa                                                                                                >  1:3141666/1‑149 (MQ=255)
                                                   aCTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTaa                                                                                                >  1:750713/1‑149 (MQ=255)
                                                    cTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAg                                                                                               >  2:453840/1‑149 (MQ=255)
                                                     tttCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGt                                                                                              <  2:4203682/149‑1 (MQ=255)
                                                          ggatggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGAtt                                                                                         <  2:3268861/149‑1 (MQ=255)
                                                          ggatggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGAtt                                                                                         <  2:2516499/149‑1 (MQ=255)
                                                           gatggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGAttt                                                                                        <  1:3120438/149‑1 (MQ=255)
                                                            atggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTg                                                                                       >  2:520309/1‑149 (MQ=255)
                                                             tggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGc                                                                                      >  1:2432691/1‑149 (MQ=255)
                                                              ggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGcc                                                                                     >  1:3880359/1‑149 (MQ=255)
                                                              ggaACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGcc                                                                                     <  1:2602963/149‑1 (MQ=255)
                                                                aaCTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAg                                                                                   >  1:1911853/1‑149 (MQ=255)
                                                                 aCTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGc                                                                                  >  2:2663666/1‑149 (MQ=255)
                                                                       cacCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGa                                                                            <  1:3341509/149‑1 (MQ=255)
                                                                        acCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCCGCGACTGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGaa                                                                           >  1:3695481/1‑149 (MQ=255)
                                                                         ccTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGaaa                                                                          <  1:3251532/149‑1 (MQ=255)
                                                                              cAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGc                                                                     >  2:4360538/1‑149 (MQ=255)
                                                                                gATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGACCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGccc                                                                   <  1:2755925/149‑1 (MQ=255)
                                                                                 aTCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCt                                                                  >  1:1344713/1‑149 (MQ=255)
                                                                                 aTCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCt                                                                  >  1:1107645/1‑149 (MQ=255)
                                                                                    ccACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCAc                                                               <  1:2301947/149‑1 (MQ=255)
                                                                                       cTGGCGATGGCTAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGt                                                            <  2:4262085/149‑1 (MQ=255)
                                                                                         ggcgatggcgaATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTa                                                          <  1:3493863/149‑1 (MQ=255)
                                                                                             atggcgaATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTAtgtg                                                      >  2:234150/1‑149 (MQ=255)
                                                                                                 cgaATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAAt                                                  <  2:2946710/149‑1 (MQ=255)
                                                                                                  gaATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAACGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATc                                                 <  2:3170918/149‑1 (MQ=255)
                                                                                                  gaATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATc                                                 >  1:3058928/1‑149 (MQ=255)
                                                                                                    aTCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCgg                                               <  1:1062107/149‑1 (MQ=255)
                                                                                                    aTCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGGGGATGTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCgg                                               >  2:3581676/1‑149 (MQ=255)
                                                                                                       cccATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCAAGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGa                                            <  1:378692/149‑1 (MQ=255)
                                                                                                        ccATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAg                                           <  2:2864857/149‑1 (MQ=255)
                                                                                                        ccATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAg                                           <  2:154288/149‑1 (MQ=255)
                                                                                                          aTAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGt                                         >  2:323128/1‑149 (MQ=255)
                                                                                                          aTAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGt                                         <  1:2269617/149‑1 (MQ=255)
                                                                                                            aaTCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGttt                                       <  2:2377694/149‑1 (MQ=255)
                                                                                                                 ctcGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCg                                  >  1:4068894/1‑149 (MQ=255)
                                                                                                                 ctcGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCg                                  <  1:4517191/149‑1 (MQ=255)
                                                                                                                   cGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTg                                <  2:2034546/149‑1 (MQ=255)
                                                                                                                    gggTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGc                               >  1:633270/1‑149 (MQ=255)
                                                                                                                     ggTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCa                              >  1:1924155/1‑149 (MQ=255)
                                                                                                                     ggTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCa                              <  1:2153472/149‑1 (MQ=255)
                                                                                                                         tCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACACCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGAc                          <  2:4489105/149‑1 (MQ=255)
                                                                                                                             tCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCttt                      <  1:2646394/149‑1 (MQ=255)
                                                                                                                             tCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCttt                      <  2:254771/149‑1 (MQ=255)
                                                                                                                              ccAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTg                     <  2:4042090/149‑1 (MQ=255)
                                                                                                                              ccAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTg                     <  2:3355458/149‑1 (MQ=255)
                                                                                                                               cAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTgc                    >  1:4140161/1‑149 (MQ=255)
                                                                                                                                  gCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTgcggc                 >  2:1440014/1‑149 (MQ=255)
                                                                                                                                        acgtcaAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGcac           >  1:923448/1‑149 (MQ=255)
                                                                                                                                           tcaAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAAc        >  1:4663603/1‑149 (MQ=255)
                                                                                                                                            caAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACt       >  1:2015140/1‑149 (MQ=255)
                                                                                                                                            caAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACt       >  2:4007147/1‑149 (MQ=255)
                                                                                                                                             aaaGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGACTGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACtt      >  1:2725366/1‑149 (MQ=255)
                                                                                                                                             aaaGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACtt      >  1:134046/1‑149 (MQ=255)
                                                                                                                                             aaaGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACtt      >  1:2526822/1‑149 (MQ=255)
                                                                                                                                             aaaGCAAGCTGGTAATGAGTCAGGTTCGACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACtt      >  2:2016036/1‑149 (MQ=255)
                                                                                                                                               aGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACTTTa    >  2:3616350/1‑149 (MQ=255)
                                                                                                                                               aGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACTTTa    <  1:3700955/149‑1 (MQ=255)
                                                                                                                                                 cACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACTTTAcc  <  2:3243810/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACAACTTTACC  >  CP009273/3483625‑3483918

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: