Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 32,527 C→A 29.5% R571S (CGT→AGT)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927332,5270CA29.5% 589.5 / 219.8 293R571S (CGT→AGT) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (97/109);  new base A (35/52);  total (132/161)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CTGCACCCGGTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCC  >  CP009273/32380‑32675
                                                                                                                                                   |                                                                                                                                                    
   cACCCGGTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTAt                                                                                                                                                  <  1:2717385/149‑1 (MQ=255)
        ggTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTc                                                                                                                                             <  2:528708/149‑1 (MQ=255)
         gTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCa                                                                                                                                            >  2:4391837/1‑149 (MQ=255)
         gTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCa                                                                                                                                            >  2:2749851/1‑149 (MQ=255)
          tttATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTCTCGGTCAg                                                                                                                                           >  2:3932207/1‑149 (MQ=255)
          tttATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAg                                                                                                                                           <  2:2331004/149‑1 (MQ=255)
               aaGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTAt                                                                                                                                      <  1:969066/149‑1 (MQ=255)
               aaGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTAt                                                                                                                                      <  1:2639700/149‑1 (MQ=255)
               aaGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTAt                                                                                                                                      <  2:968691/149‑1 (MQ=255)
                   gcgTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCAGGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGaa                                                                                                                                  >  2:1505424/1‑149 (MQ=255)
                      tGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTc                                                                                                                               <  2:946512/149‑1 (MQ=255)
                       ggATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTGTCGGTCAGGGTATCGAATTCg                                                                                                                              <  1:566711/149‑1 (MQ=255)
                          tACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGact                                                                                                                           <  2:2009542/149‑1 (MQ=255)
                            ccTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGCCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGactac                                                                                                                         <  1:1329441/149‑1 (MQ=255)
                              tgtgCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGTGTATCGAATTCGACTACTg                                                                                                                       >  2:1894627/1‑149 (MQ=255)
                              tgtgCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTg                                                                                                                       >  2:2451222/1‑149 (MQ=255)
                              tgtgCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTg                                                                                                                       >  2:46795/1‑149 (MQ=255)
                                tgCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGtt                                                                                                                     <  1:1937759/149‑1 (MQ=255)
                                    gCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGt                                                                                                                 >  1:2245772/1‑149 (MQ=255)
                                    gCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGt                                                                                                                 >  1:3399091/1‑149 (MQ=255)
                                      agagTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTac                                                                                                               <  1:1190806/149‑1 (MQ=255)
                                      agagTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTac                                                                                                               <  1:2533032/149‑1 (MQ=255)
                                       gagTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTaca                                                                                                              >  2:2194017/1‑149 (MQ=255)
                                        agTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTacac                                                                                                             <  1:3783497/149‑1 (MQ=255)
                                           tCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGcc                                                                                                          <  2:3743391/149‑1 (MQ=255)
                                           tCGCCACAGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCTGGGTATCGAATTCGACTAATGTTGCGTACACGcc                                                                                                          <  2:1354050/149‑1 (MQ=255)
                                            cGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCt                                                                                                         <  1:751522/149‑1 (MQ=255)
                                             gCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTc                                                                                                        <  1:4226395/149‑1 (MQ=255)
                                             gCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTc                                                                                                        <  1:4609560/149‑1 (MQ=255)
                                             gCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTc                                                                                                        <  2:2625473/149‑1 (MQ=255)
                                               cACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGc                                                                                                      <  2:523966/149‑1 (MQ=255)
                                                aCCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCt                                                                                                     >  1:1404155/1‑149 (MQ=255)
                                                aCCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCt                                                                                                     <  1:616719/149‑1 (MQ=255)
                                                  cGACACCGCTTACATGTACTCCACTTATGACGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTgg                                                                                                   <  2:630826/149‑1 (MQ=255)
                                                  cGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTgg                                                                                                   <  2:2945583/149‑1 (MQ=255)
                                                  cGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTgg                                                                                                   <  1:1162339/149‑1 (MQ=255)
                                                    acacCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGgcg                                                                                                 >  2:2069316/1‑149 (MQ=255)
                                                         gCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTgcg                                                                                            >  2:4090498/1‑149 (MQ=255)
                                                         gCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTgcg                                                                                            >  1:3517018/1‑149 (MQ=255)
                                                           ttACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGcgcg                                                                                          >  2:1532512/1‑149 (MQ=255)
                                                           ttACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGAGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGcgcg                                                                                          >  2:4721104/1‑149 (MQ=255)
                                                             aCATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGaa                                                                                        >  1:4374455/1‑149 (MQ=255)
                                                               aTGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGa                                                                                      >  1:504413/1‑149 (MQ=255)
                                                               aTGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGa                                                                                      <  2:3357415/149‑1 (MQ=255)
                                                               aTGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGa                                                                                      <  2:1180175/149‑1 (MQ=255)
                                                                  tACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACgg                                                                                   >  1:2060771/1‑149 (MQ=255)
                                                                   aCTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGt                                                                                  <  1:1692390/149‑1 (MQ=255)
                                                                      ccACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTAc                                                                               <  1:41963/149‑1 (MQ=255)
                                                                       cACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACg                                                                              >  1:2925912/1‑149 (MQ=255)
                                                                        aCTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGa                                                                             <  2:2571003/149‑1 (MQ=255)
                                                                        aCTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGa                                                                             >  1:3647353/1‑149 (MQ=255)
                                                                            aTGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCACGAAGACGGTTACGAAAcc                                                                         >  2:2767126/1‑149 (MQ=255)
                                                                             tGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCa                                                                        >  1:2734600/1‑149 (MQ=255)
                                                                                agaagaGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTa                                                                     <  2:1236890/149‑1 (MQ=255)
                                                                                 gaagaGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTAt                                                                    >  1:2992543/1‑149 (MQ=255)
                                                                                 gaagaGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTAt                                                                    >  1:2896455/1‑149 (MQ=255)
                                                                                   agagTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATgg                                                                  <  1:2482451/149‑1 (MQ=255)
                                                                                    gagTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGt                                                                 >  1:1551760/1‑149 (MQ=255)
                                                                                     agTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGtt                                                                <  2:1924253/149‑1 (MQ=255)
                                                                                     agTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGtt                                                                <  2:2921319/149‑1 (MQ=255)
                                                                                     agTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGtt                                                                <  2:4288865/149‑1 (MQ=255)
                                                                                     agTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGtt                                                                <  2:1948025/149‑1 (MQ=255)
                                                                                      gTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTa                                                               >  1:4262963/1‑149 (MQ=255)
                                                                                        gcgaagcgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAAc                                                             <  2:4536495/149‑1 (MQ=255)
                                                                                          gaagcgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTg                                                           >  1:4074418/1‑149 (MQ=255)
                                                                                          gaagcgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTg                                                           <  1:4403197/149‑1 (MQ=255)
                                                                                             gcgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCTGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTaa                                                        <  2:3170890/149‑1 (MQ=255)
                                                                                             gcgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTaa                                                        <  2:634216/149‑1 (MQ=255)
                                                                                              cgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAc                                                       <  2:788069/149‑1 (MQ=255)
                                                                                              cgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAc                                                       <  2:2529337/149‑1 (MQ=255)
                                                                                              cgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAc                                                       <  1:698186/149‑1 (MQ=255)
                                                                                              cgaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAc                                                       <  2:3524379/149‑1 (MQ=255)
                                                                                               gaaTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAAcc                                                      >  2:3637160/1‑149 (MQ=255)
                                                                                                 aTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCg                                                    >  1:1814906/1‑149 (MQ=255)
                                                                                                 aTCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCg                                                    <  2:3838997/149‑1 (MQ=255)
                                                                                                   ccGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGa                                                  >  2:3455068/1‑149 (MQ=255)
                                                                                                     gTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGaaa                                                >  2:853832/1‑149 (MQ=255)
                                                                                                     gTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGaaa                                                <  1:122866/149‑1 (MQ=255)
                                                                                                          ccgaccgTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtc                                           >  2:828719/1‑149 (MQ=255)
                                                                                                           cgaccgTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCAGTATGTTTAACTGTAACCCGGAAACCGtct                                          >  2:2257896/1‑149 (MQ=255)
                                                                                                            gaccgTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGtctc                                         >  2:2419056/1‑149 (MQ=255)
                                                                                                             accgTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTcc                                        <  2:3030121/149‑1 (MQ=255)
                                                                                                               cgTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCAc                                      <  2:1595049/149‑1 (MQ=255)
                                                                                                                 tGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCg                                    >  2:4313507/1‑149 (MQ=255)
                                                                                                                   aaaaaaTCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGAc                                  <  2:2824471/149‑1 (MQ=255)
                                                                                                                   aaaaaaTCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGAc                                  >  1:3582778/1‑149 (MQ=255)
                                                                                                                   aaaaaaTCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGAc                                  <  1:4571641/149‑1 (MQ=255)
                                                                                                                       aaTCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACg                              >  1:434370/1‑149 (MQ=255)
                                                                                                                            tGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACt                         <  2:4479882/149‑1 (MQ=255)
                                                                                                                             ggTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACtt                        <  2:4235482/149‑1 (MQ=255)
                                                                                                                               tGCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTcc                      >  2:2459187/1‑149 (MQ=255)
                                                                                                                                gCTTGGCGGCGGCCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCg                     <  1:2818068/149‑1 (MQ=255)
                                                                                                                                  ttGGCGGCGGCCCGAACCGTATCGGTCAGGGTCTCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCTTCTCCACCGACTACGACACTTCCGAc                   >  1:2300025/1‑149 (MQ=255)
                                                                                                                                       ggcggcCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCct              >  1:3168877/1‑149 (MQ=255)
                                                                                                                                        gcggcCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctc             >  1:4659799/1‑149 (MQ=255)
                                                                                                                                         cggcCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctct            >  1:3486126/1‑149 (MQ=255)
                                                                                                                                         cggcCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCctct            >  1:1098727/1‑149 (MQ=255)
                                                                                                                                           gcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTAc          >  2:1696356/1‑149 (MQ=255)
                                                                                                                                           gcCCGAACAGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTAc          >  1:450026/1‑149 (MQ=255)
                                                                                                                                                aaCCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACAGCCTCTACTTCGa     >  1:4473422/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCACCCGGTTTATAAGCGCGTGGATACCTGTGCGGCAGAGTTCGCCACCGACACCGCTTACATGTACTCCACTTATGAAGAAGAGTGCGAAGCGAATCCGTCTACCGACCGTGAAAAAATCATGGTGCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCC  >  CP009273/32380‑32675

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: