Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 509,737 G→A 87.6% E97K (GAG→AAG)  cueR → copper‑responsive regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273509,7370GA87.6% 601.5 / 68.8 210E97K (GAG→AAG) cueRcopper‑responsive regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base A (66/118);  total (76/134)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 214 total aligned reads displayed.

CCTTACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATG  >  CP009273/509591‑509880
                                                                                                                                                  |                                                                                                                                               
ccgtACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTgagga                                                                                                                                                <  2:1486618/145‑1 (MQ=255)
ccTTACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTgagga                                                                                                                                                <  2:3227743/148‑1 (MQ=255)
 cTTACGGCGCCCGGCACGGCTGGTGGGTTTTAACCTGGAAGAGCGCGGCGAGCTGGTGAATGTGTTTAACGACCCGCAGCGGCACAGCCCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTgagaag                                                                                                                                               <  1:576257/148‑1 (MQ=255)
  gtACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGc                                                                                                                                              <  2:1916756/147‑1 (MQ=255)
     cTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGc                                                                                                                                           >  1:1778970/1‑148 (MQ=255)
     cTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGc                                                                                                                                           >  1:1778743/1‑148 (MQ=255)
      tGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCa                                                                                                                                          <  1:1480321/148‑1 (MQ=255)
      ggcgcCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGATGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCa                                                                                                                                          <  1:2924290/147‑1 (MQ=255)
       gcgcCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCaa                                                                                                                                         <  1:224830/148‑1 (MQ=255)
          ccAGGCCCGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATcc                                                                                                                                      <  1:2858483/148‑1 (MQ=255)
           cAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCa                                                                                                                                     <  2:1755841/148‑1 (MQ=255)
            aGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCAt                                                                                                                                    <  1:130067/148‑1 (MQ=255)
             ggCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATg                                                                                                                                   >  2:1790856/1‑148 (MQ=255)
                   gCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGAc                                                                                                                             <  2:2840566/148‑1 (MQ=255)
                         gggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctg                                                                                                                       >  2:1055401/1‑148 (MQ=255)
                          gtcTTTACCCTGGAAGAGAGCGGCGCGCTGTTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctgc                                                                                                                      <  1:864546/146‑1 (MQ=255)
                          ggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctgc                                                                                                                      <  1:1030104/148‑1 (MQ=255)
                          ggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctgc                                                                                                                      <  2:1589162/148‑1 (MQ=255)
                           gCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATACGCGACCAgctgct                                                                                                                     >  2:488796/1‑148 (MQ=255)
                              ttAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGc                                                                                                                  <  2:2524461/148‑1 (MQ=255)
                                aaCCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCAc                                                                                                                <  2:24457/148‑1 (MQ=255)
                                 aCCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACt                                                                                                               >  2:622315/1‑148 (MQ=255)
                                 aCCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACt                                                                                                               <  1:1088146/148‑1 (MQ=255)
                                    tGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGc                                                                                                            >  1:2751503/1‑148 (MQ=255)
                                    tGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATAGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGc                                                                                                            >  2:541067/1‑148 (MQ=255)
                                     ggAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCg                                                                                                           >  1:1256743/1‑148 (MQ=255)
                                      gAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGa                                                                                                          >  1:1537701/1‑148 (MQ=255)
                                         gagagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtg                                                                                                       <  2:2079253/148‑1 (MQ=255)
                                          agagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                      <  2:1164928/148‑1 (MQ=255)
                                              cggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgc                                                                                                  >  2:2655958/1‑148 (MQ=255)
                                              cggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgc                                                                                                  >  2:1006106/1‑148 (MQ=255)
                                               ggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgcc                                                                                                 <  1:3132713/148‑1 (MQ=255)
                                               ggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgcc                                                                                                 <  1:1316439/148‑1 (MQ=255)
                                                    gCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGc                                                                                            >  1:236768/1‑148 (MQ=255)
                                                     cTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCg                                                                                           <  2:2870069/148‑1 (MQ=255)
                                                          gAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGAc                                                                                      <  2:717157/148‑1 (MQ=255)
                                                            aTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAg                                                                                    <  1:1406280/148‑1 (MQ=255)
                                                              cTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcg                                                                                  >  2:417721/1‑148 (MQ=255)
                                                               tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcgc                                                                                 <  2:949651/148‑1 (MQ=255)
                                                               tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcgc                                                                                 >  1:2931695/1‑148 (MQ=255)
                                                               tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcgc                                                                                 <  1:2448533/148‑1 (MQ=255)
                                                                  ttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGa                                                                              <  2:677026/148‑1 (MQ=255)
                                                                   tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGAc                                                                             >  1:611458/1‑148 (MQ=255)
                                                                     aCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                           <  2:1447647/148‑1 (MQ=255)
                                                                     aCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                           >  2:1248530/1‑148 (MQ=255)
                                                                      cGACCCGCAGCGGCACAGCGCCGCCGTCAAACGGCGCACGCTGGCGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCAGGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGc                                                                          <  2:2876701/148‑1 (MQ=255)
                                                                      cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGc                                                                          >  2:1740532/1‑148 (MQ=255)
                                                                           cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAt                                                                     >  1:2878742/1‑148 (MQ=255)
                                                                           cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAt                                                                     <  2:2092023/148‑1 (MQ=255)
                                                                           cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAt                                                                     >  1:2037108/1‑148 (MQ=255)
                                                                            gCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAtt                                                                    <  2:1806656/148‑1 (MQ=255)
                                                                             cAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                   <  2:315893/148‑1 (MQ=255)
                                                                             cAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                   <  2:1219775/148‑1 (MQ=255)
                                                                             cAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                   >  1:2494164/1‑148 (MQ=255)
                                                                              aGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTAt                                                                  <  1:3221690/148‑1 (MQ=255)
                                                                              aGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTAt                                                                  <  2:1654104/148‑1 (MQ=255)
                                                                               gCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATc                                                                 <  2:2275313/148‑1 (MQ=255)
                                                                                cGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCg                                                                >  2:2372834/1‑148 (MQ=255)
                                                                                   cacaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                             <  2:159244/148‑1 (MQ=255)
                                                                                     caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAt                                                           <  1:650448/148‑1 (MQ=255)
                                                                                       gcgcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtct                                                         <  2:2947732/148‑1 (MQ=255)
                                                                                        cgcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctc                                                        >  1:1748104/1‑148 (MQ=255)
                                                                                         gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctct                                                       <  2:1536607/148‑1 (MQ=255)
                                                                                         gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctct                                                       <  1:38646/148‑1 (MQ=255)
                                                                                         gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctct                                                       <  2:794297/148‑1 (MQ=255)
                                                                                          ccGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctctc                                                      <  1:3309303/148‑1 (MQ=255)
                                                                                                 cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgctgc                                               >  1:1015339/1‑148 (MQ=255)
                                                                                                     cGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAAACCATGCGCGACCAGCAGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtc                                           >  2:2192666/1‑148 (MQ=255)
                                                                                                       gcgcACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtcat                                         <  1:622315/148‑1 (MQ=255)
                                                                                                          cACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTcatcat                                      <  1:1002762/148‑1 (MQ=255)
                                                                                                          cACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTcatcat                                      <  1:2655958/148‑1 (MQ=255)
                                                                                                              cTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCggg                                  <  1:776792/148‑1 (MQ=255)
                                                                                                                 gagaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAg                               >  1:3049582/1‑148 (MQ=255)
                                                                                                                   gaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCACTCCGGCTGCTGTCATCATCGGGCAggg                             >  2:2756219/1‑148 (MQ=255)
                                                                                                                    aaGGTGGCGGAGATCGCACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGt                            <  2:3050776/148‑1 (MQ=255)
                                                                                                                     aGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTg                           <  2:2318556/148‑1 (MQ=255)
                                                                                                                      ggTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGa                          <  2:2493831/148‑1 (MQ=255)
                                                                                                                      ggTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGa                          >  2:1359158/1‑148 (MQ=255)
                                                                                                                      ggTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGa                          >  2:1147208/1‑148 (MQ=255)
                                                                                                                       gTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAt                         <  2:824178/148‑1 (MQ=255)
                                                                                                                       gTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAt                         <  1:1884699/148‑1 (MQ=255)
                                                                                                                           cGGAGATCGAACGACACATTGAGAAGCTGCACTCCCTGCGCGACCAGCTGCTGGCACGGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                     <  2:1865033/148‑1 (MQ=255)
                                                                                                                           cGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                     <  1:104343/148‑1 (MQ=255)
                                                                                                                               gaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacga                 <  1:1790856/148‑1 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                >  2:861187/1‑148 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                <  2:1778959/148‑1 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                <  1:936684/148‑1 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                >  1:1963441/1‑148 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                <  2:1819879/148‑1 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                >  1:2021671/1‑148 (MQ=255)
                                                                                                                                aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                >  1:2987278/1‑148 (MQ=255)
                                                                                                                                      cGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCt          >  1:2863631/1‑148 (MQ=255)
                                                                                                                                      cGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCt          <  1:497712/148‑1 (MQ=255)
                                                                                                                                        acacATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCg        >  1:3305696/1‑148 (MQ=255)
                                                                                                                                        acacATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCg        >  2:746077/1‑148 (MQ=255)
                                                                                                                                         cacaTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCgg       <  1:2604175/148‑1 (MQ=255)
                                                                                                                                         cacaTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCgg       <  2:1203168/148‑1 (MQ=255)
                                                                                                                                          acaTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGa      <  2:3282424/148‑1 (MQ=255)
                                                                                                                                           caTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGat     >  2:127333/1‑148 (MQ=255)
                                                                                                                                           caTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGat     <  1:2267526/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCTTACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATG  >  CP009273/509591‑509880

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: