Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 509,737 G→A 100% 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,7370GA97.3% 590.9 / 5.5 182E97K (GAG→AAG) cueRcopper‑responsive regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base A (65/112);  total (65/117)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 183 total aligned reads displayed.

TTACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGT  >  CP009273/509593‑509881
                                                                                                                                                |                                                                                                                                                
ttACTGCGCCCGGCACGGCTGGTGGGATTTAACCTGGAAGCGAGCGGCGCGCTGGTGAATCTTTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGc                                                                                                                                               <  1:30727/148‑1 (MQ=255)
ttACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGc                                                                                                                                               >  2:2629139/1‑148 (MQ=255)
ttACGGCGCCCGGCCCGGCTGGTGGGCTTTACCCTGGACGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGc                                                                                                                                               <  1:30226/148‑1 (MQ=255)
    tGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCa                                                                                                                                           >  2:121806/1‑148 (MQ=255)
     gcgcCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCaa                                                                                                                                          <  2:3363254/148‑1 (MQ=255)
       gcgaGGCACGGCAGGTGGGCTTTAAGCTGGTTGAGAGGGGCGAGCTGGTGAATCTGTTTTACGAGCCGCAGCGGCACCGCGACGTCGTCAAACGGCGCACTCTGGAGAAGGTGGCGGAGATCGAACTACACATTTAGGTGCTGCAATc                                                                                                                                        <  2:2525841/145‑1 (MQ=255)
           ggCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATg                                                                                                                                    <  2:2771252/148‑1 (MQ=255)
            gCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATgc                                                                                                                                   >  2:499857/1‑148 (MQ=255)
                ggCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGa                                                                                                                               <  2:1960256/148‑1 (MQ=255)
                       gggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctg                                                                                                                        <  2:2293859/148‑1 (MQ=255)
                        ggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctgc                                                                                                                       >  1:1324297/1‑148 (MQ=255)
                         gCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAgctgct                                                                                                                      >  2:3266198/1‑148 (MQ=255)
                           tttAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTgg                                                                                                                    >  2:2044577/1‑148 (MQ=255)
                             tAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCa                                                                                                                  <  1:2762438/148‑1 (MQ=255)
                                ccTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTg                                                                                                               <  1:276285/148‑1 (MQ=255)
                                 cTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCACACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTgg                                                                                                              >  2:3058352/1‑148 (MQ=255)
                                   ggAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCg                                                                                                            >  2:1146545/1‑148 (MQ=255)
                                     aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACTACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGaa                                                                                                          >  1:1024495/1‑148 (MQ=255)
                                     aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGaa                                                                                                          >  1:492447/1‑148 (MQ=255)
                                     aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGaa                                                                                                          >  1:2256808/1‑148 (MQ=255)
                                     aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGaa                                                                                                          >  1:491287/1‑148 (MQ=255)
                                        agagCGGCGAGCTGGTTAATCTTTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                       <  1:958718/148‑1 (MQ=255)
                                        agagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGTGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                       >  1:1759556/1‑148 (MQ=255)
                                        agagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                       <  1:664771/148‑1 (MQ=255)
                                         gagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcc                                                                                                      >  1:2878007/1‑148 (MQ=255)
                                          agCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcct                                                                                                     >  1:2055204/1‑148 (MQ=255)
                                            cggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgc                                                                                                   >  2:2443525/1‑148 (MQ=255)
                                            cggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgcctgc                                                                                                   <  2:1772411/148‑1 (MQ=255)
                                              gcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGccc                                                                                                 <  2:738862/148‑1 (MQ=255)
                                              gcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGccc                                                                                                 >  1:1918525/1‑148 (MQ=255)
                                               cgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCt                                                                                                <  1:156271/148‑1 (MQ=255)
                                                 aGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTgg                                                                                              >  1:614117/1‑148 (MQ=255)
                                                   cTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCg                                                                                            >  2:417547/1‑148 (MQ=255)
                                                     ggTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                          >  2:2074245/1‑148 (MQ=255)
                                                     ggTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                          >  1:1689503/1‑148 (MQ=255)
                                                       tGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGa                                                                                        >  2:761213/1‑148 (MQ=255)
                                                         aaTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACa                                                                                      <  2:2994327/148‑1 (MQ=255)
                                                         aaTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACa                                                                                      <  2:1548754/148‑1 (MQ=255)
                                                          aTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAg                                                                                     <  1:2076425/148‑1 (MQ=255)
                                                          aTCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAg                                                                                     >  1:2655613/1‑148 (MQ=255)
                                                           tCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgc                                                                                    >  1:1875506/1‑148 (MQ=255)
                                                             tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcgc                                                                                  <  2:1734860/148‑1 (MQ=255)
                                                              gTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGcc                                                                                 >  1:2113186/1‑148 (MQ=255)
                                                                 tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGAc                                                                              >  2:2766022/1‑148 (MQ=255)
                                                                   aCGACCCGCAGCGGCACAGCGCCGACGTCAATCGGCGCACGCTGGAGAATGTGGCGGAGATCGATCGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                            <  1:2589556/148‑1 (MQ=255)
                                                                    cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGc                                                                           >  2:2727190/1‑148 (MQ=255)
                                                                       cccGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCg                                                                        <  2:2234893/148‑1 (MQ=255)
                                                                       cccGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCg                                                                        <  2:2687181/148‑1 (MQ=255)
                                                                         cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAt                                                                      <  2:3117206/148‑1 (MQ=255)
                                                                          gCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAtt                                                                     <  1:1534380/148‑1 (MQ=255)
                                                                           cAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                    <  2:2851702/148‑1 (MQ=255)
                                                                                 cacaGCGCCGACGTCAAACGGCGCACGCTGTATAAGGTGGCGGATATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                              <  1:612005/148‑1 (MQ=255)
                                                                                 cacaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                              <  2:1052888/148‑1 (MQ=255)
                                                                                 cacaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                              >  2:360674/1‑148 (MQ=255)
                                                                                 cacaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                              <  1:2727190/148‑1 (MQ=255)
                                                                                   caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAt                                                            <  1:1805624/148‑1 (MQ=255)
                                                                                   caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAt                                                            <  1:1804555/148‑1 (MQ=255)
                                                                                   caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAt                                                            <  1:1686430/148‑1 (MQ=255)
                                                                                      cgcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctc                                                         >  2:2347881/1‑148 (MQ=255)
                                                                                       gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctct                                                        <  1:1117054/148‑1 (MQ=255)
                                                                                       gcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctct                                                        <  1:2714563/148‑1 (MQ=255)
                                                                                        ccGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctctc                                                       <  1:343720/148‑1 (MQ=255)
                                                                                        ccGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctctc                                                       <  1:135306/148‑1 (MQ=255)
                                                                                        ccGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtctctc                                                       <  1:134747/148‑1 (MQ=255)
                                                                                            cGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgc                                                   <  1:2044577/148‑1 (MQ=255)
                                                                                               cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgctgc                                                >  2:1621909/1‑148 (MQ=255)
                                                                                               cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgctgc                                                >  2:2171162/1‑148 (MQ=255)
                                                                                                  aCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGt                                             <  2:1242889/148‑1 (MQ=255)
                                                                                                  aCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGt                                             <  2:2992198/148‑1 (MQ=255)
                                                                                                         aCGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCatcatc                                      <  1:1875480/148‑1 (MQ=255)
                                                                                                          cGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCg                                     >  1:1673308/1‑148 (MQ=255)
                                                                                                            cTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCggg                                   <  2:1557123/148‑1 (MQ=255)
                                                                                                               gagaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAg                                >  2:969346/1‑148 (MQ=255)
                                                                                                               gagaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAg                                >  1:1145486/1‑148 (MQ=255)
                                                                                                               gagaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAg                                >  1:2405389/1‑148 (MQ=255)
                                                                                                               gagaAGGTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAg                                <  1:2467630/148‑1 (MQ=255)
                                                                                                                    ggTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCCGGGTGa                           >  2:3148529/1‑148 (MQ=255)
                                                                                                                    ggTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGa                           <  2:391737/148‑1 (MQ=255)
                                                                                                                     gTGGCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGAt                          <  2:1635528/148‑1 (MQ=255)
                                                                                                                        gCGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTaa                       <  1:2731612/148‑1 (MQ=255)
                                                                                                                         cGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                      <  2:1659870/148‑1 (MQ=255)
                                                                                                                         cGGAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAag                      <  2:1790495/148‑1 (MQ=255)
                                                                                                                          ggAGATCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAaga                     <  1:732254/148‑1 (MQ=255)
                                                                                                                            agaTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacg                   >  1:3279316/1‑148 (MQ=255)
                                                                                                                            agaTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacg                   <  1:121806/148‑1 (MQ=255)
                                                                                                                              aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                 >  2:2272803/1‑148 (MQ=255)
                                                                                                                              aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                 <  1:1547408/148‑1 (MQ=255)
                                                                                                                              aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                 <  2:2113186/148‑1 (MQ=255)
                                                                                                                              aTCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgag                 <  2:789158/148‑1 (MQ=255)
                                                                                                                               tCGAACGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgaga                <  2:2404320/148‑1 (MQ=255)
                                                                                                                                  aaCGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACgg             <  1:751340/148‑1 (MQ=255)
                                                                                                                                  aaCGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACgg             <  1:751417/148‑1 (MQ=255)
                                                                                                                                  aaCGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACgg             <  2:1296835/148‑1 (MQ=255)
                                                                                                                                   aCGACACATTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGc            <  2:1339185/148‑1 (MQ=255)
                                                                                                                                       cacaTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCgg        <  1:2649349/148‑1 (MQ=255)
                                                                                                                                        acaTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGa       <  2:2126571/148‑1 (MQ=255)
                                                                                                                                         caTTGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGat      >  2:2263874/1‑148 (MQ=255)
                                                                                                                                            tGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATg   <  1:917239/148‑1 (MQ=255)
                                                                                                                                            tGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATg   <  1:1400246/148‑1 (MQ=255)
                                                                                                                                            tGAGAAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATg   <  1:2930258/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTACTGCGCCAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGT  >  CP009273/509593‑509881

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: