Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,176,051 C→A 6.3% S48* (TCA→TAA)  zupT → zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,176,0510CA 6.3% 1072.9 / 30.8 286S48* (TCA→TAA) zupTzinc transporter
Reads supporting (aligned to +/- strand):  ref base C (131/137);  new base A (3/15);  total (134/152)
Fisher's exact test for biased strand distribution p-value = 1.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 311 total aligned reads displayed.

GTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCG  >  CP009273/3175902‑3176199
                                                                                                                                                     |                                                                                                                                                    
gTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatctc                                                                                                                                                      <  1:1656214/150‑1 (MQ=255)
 ttGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatctca                                                                                                                                                     <  1:2815999/150‑1 (MQ=255)
 ttGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatctca                                                                                                                                                     <  1:2365226/150‑1 (MQ=255)
  tGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatctcat                                                                                                                                                    <  2:230942/150‑1 (MQ=255)
  tGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcatctcat                                                                                                                                                    <  2:1298477/150‑1 (MQ=255)
   gTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCAtt                                                                                                                                                   <  2:4159181/150‑1 (MQ=255)
      gATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAAt                                                                                                                                                >  1:1023249/1‑150 (MQ=255)
      gATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAAt                                                                                                                                                >  2:825137/1‑150 (MQ=255)
        tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATgg                                                                                                                                              <  1:2101403/150‑1 (MQ=255)
        tGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATgg                                                                                                                                              <  1:1696961/150‑1 (MQ=255)
         gTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGa                                                                                                                                             >  1:4540590/1‑150 (MQ=255)
          tCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGaa                                                                                                                                            <  1:551567/150‑1 (MQ=255)
          tCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTAACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGaa                                                                                                                                            >  1:4173954/1‑150 (MQ=255)
           cAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGaaa                                                                                                                                           >  2:471076/1‑150 (MQ=255)
               aCCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCt                                                                                                                                       >  2:291091/1‑150 (MQ=255)
               aCCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCt                                                                                                                                       >  1:1485295/1‑150 (MQ=255)
                ccTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCtt                                                                                                                                      <  1:721943/150‑1 (MQ=255)
                ccTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCtt                                                                                                                                      <  2:1551413/150‑1 (MQ=255)
                   ctcATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTAATTAATGGAAATGCTTCCt                                                                                                                                   <  1:1143860/150‑1 (MQ=255)
                    tcATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTg                                                                                                                                  <  1:894790/150‑1 (MQ=255)
                      aTTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTAATTAATGGAAATGCTTCCTGcc                                                                                                                                >  1:3381108/1‑150 (MQ=255)
                       ttCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCg                                                                                                                               <  1:1287776/150‑1 (MQ=255)
                          tGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCAc                                                                                                                            >  1:4377148/1‑150 (MQ=255)
                              cATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagc                                                                                                                        <  2:1003019/150‑1 (MQ=255)
                              cATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagc                                                                                                                        <  2:2159926/150‑1 (MQ=255)
                                tACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGGTCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                                                      <  1:3984263/150‑1 (MQ=255)
                                tACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcag                                                                                                                      <  2:2229621/150‑1 (MQ=255)
                                 aCTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                     <  2:408919/150‑1 (MQ=255)
                                 aCTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                     >  2:3814788/1‑150 (MQ=255)
                                 aCTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                     >  2:3814709/1‑150 (MQ=255)
                                  cTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCt                                                                                                                    >  2:1997046/1‑150 (MQ=255)
                                  cTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCt                                                                                                                    <  1:1017802/150‑1 (MQ=255)
                                  cTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTACGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCt                                                                                                                    >  1:2937161/1‑150 (MQ=255)
                                   tGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTg                                                                                                                   >  1:1355403/1‑150 (MQ=255)
                                       gggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAgg                                                                                                               <  1:3140262/150‑1 (MQ=255)
                                       gggggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAgg                                                                                                               >  2:337781/1‑150 (MQ=255)
                                          ggCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAAt                                                                                                            >  1:1716780/1‑150 (MQ=255)
                                            cAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGt                                                                                                          >  1:835230/1‑150 (MQ=255)
                                             aGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGGATGTc                                                                                                         >  2:3698239/1‑150 (MQ=255)
                                              gCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCg                                                                                                        >  1:3014860/1‑150 (MQ=255)
                                                     ttATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTgtgt                                                                                                 >  2:3084905/1‑150 (MQ=255)
                                                            cgcgTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTa                                                                                          <  2:2102541/150‑1 (MQ=255)
                                                             gcgTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTAt                                                                                         >  1:1383347/1‑150 (MQ=255)
                                                              cgTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                                                        <  2:3889988/150‑1 (MQ=255)
                                                               gTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                                                       <  2:2595033/150‑1 (MQ=255)
                                                                    tcGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                                                  <  1:3844471/150‑1 (MQ=255)
                                                                    tcGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                                                  <  2:1832399/150‑1 (MQ=255)
                                                                     cGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                                                 >  1:2777111/1‑150 (MQ=255)
                                                                     cGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                                                 <  1:1552856/150‑1 (MQ=255)
                                                                      ggCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGACTGTCGCCTGTGTTGGGTTAGGGAATGttg                                                                                >  2:3018837/1‑149 (MQ=255)
                                                                      ggCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGttt                                                                                >  2:737387/1‑150 (MQ=255)
                                                                       gCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                                                               <  1:1997726/150‑1 (MQ=255)
                                                                           tctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGTAATGCTTACTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCtt                                                                           <  1:355680/150‑1 (MQ=255)
                                                                             tcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCg                                                                         <  1:3397958/150‑1 (MQ=255)
                                                                                ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTc                                                                      <  2:4508979/150‑1 (MQ=255)
                                                                                 gCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                                                     <  2:2493426/150‑1 (MQ=255)
                                                                                  cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTg                                                                    <  2:1108942/150‑1 (MQ=255)
                                                                                  cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTg                                                                    <  1:2367063/150‑1 (MQ=255)
                                                                                  cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTg                                                                    >  1:2387543/1‑150 (MQ=255)
                                                                                  cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTg                                                                    <  2:631249/150‑1 (MQ=255)
                                                                                  cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATATTCGGTCTg                                                                    >  1:1932581/1‑150 (MQ=255)
                                                                                     aaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTAATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCtt                                                                 <  2:1477573/150‑1 (MQ=255)
                                                                                     aaaCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGGTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCGTCGGTCTGCtt                                                                 >  1:1773855/1‑150 (MQ=255)
                                                                                        cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGc                                                              <  2:3029190/150‑1 (MQ=255)
                                                                                         ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCt                                                             >  2:2756573/1‑150 (MQ=255)
                                                                                             gAACAGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATATCATTAATGGAAATGATTCCTGCCGCACTAGCAGCTTAAGGAATGTCGCCTGTGTTTGGTTATGGAATTTTTATATTCGGTCTGCTTGGCTAttt                                                         >  2:3388884/1‑150 (MQ=255)
                                                                                               aCCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTg                                                       <  1:3124056/150‑1 (MQ=255)
                                                                                                   ctgtcTGGCGTTTTCGCTTGGTTGTGCGGCGGTGATCATGTTGCTCATCTCATTACTGGAACTGCTGCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                                                   <  2:4371978/146‑1 (MQ=255)
                                                                                                     tACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTgg                                                 <  2:4540840/150‑1 (MQ=255)
                                                                                                      aCTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTTGCCTGGa                                                >  1:2823081/1‑150 (MQ=255)
                                                                                                        tGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc                                              <  2:566080/150‑1 (MQ=255)
                                                                                                        tGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc                                              <  2:2385297/150‑1 (MQ=255)
                                                                                                            gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCCTCTCATTAATGGAAAAGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCAt                                          >  2:4285903/1‑150 (MQ=255)
                                                                                                            gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCAt                                          >  2:1494885/1‑150 (MQ=255)
                                                                                                            gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTAATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCAt                                          <  2:786036/150‑1 (MQ=255)
                                                                                                             ttttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGTCCTGGACCGCCTg                                         >  1:2672408/1‑150 (MQ=255)
                                                                                                             ttttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATg                                         <  2:1671087/150‑1 (MQ=255)
                                                                                                               ttCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGtt                                       <  1:3518826/150‑1 (MQ=255)
                                                                                                                  gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc                                    >  2:4059903/1‑150 (MQ=255)
                                                                                                                  gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc                                    >  2:831646/1‑150 (MQ=255)
                                                                                                                   cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca                                   >  1:251186/1‑150 (MQ=255)
                                                                                                                   cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca                                   <  2:1260107/150‑1 (MQ=255)
                                                                                                                     aGGTTTTGCGGCGGGGATCATGTTGCTCATCTAATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca                                 <  2:3606271/150‑1 (MQ=255)
                                                                                                                        ttttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGc                              >  1:1127825/1‑150 (MQ=255)
                                                                                                                         tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCt                             <  1:3720413/150‑1 (MQ=255)
                                                                                                                          ttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTc                            <  2:3317952/150‑1 (MQ=255)
                                                                                                                             cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATc                         >  1:1357023/1‑150 (MQ=255)
                                                                                                                              ggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATcc                        <  1:1091210/150‑1 (MQ=255)
                                                                                                                              ggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATcc                        <  2:1597672/150‑1 (MQ=255)
                                                                                                                              ggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATATTGGCCTGGACCGCATGTTGCCACATGCTCATcc                        <  2:399689/150‑1 (MQ=255)
                                                                                                                                        aTGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTa              >  1:4428790/1‑150 (MQ=255)
                                                                                                                                          gTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCGGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt            >  2:4231921/1‑150 (MQ=255)
                                                                                                                                            tGCTCATCTAATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc          <  2:3381108/150‑1 (MQ=255)
                                                                                                                                             gCTCATCTCATTAATGGACATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCa         <  2:2663062/150‑1 (MQ=255)
                                                                                                                                             gCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCa         >  1:2558712/1‑150 (MQ=255)
                                                                                                                                             gCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCa         >  1:1856342/1‑150 (MQ=255)
                                                                                                                                                catctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaaaa      <  2:4100872/150‑1 (MQ=255)
                                                                                                                                                 atctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaaaaa     >  1:3161206/1‑150 (MQ=255)
                                                                                                                                                   ctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATc   <  2:1802420/150‑1 (MQ=255)
                                                                                                                                                    tcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCg  <  2:1355403/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGGGGCAGCCACGTTTATTGGCGCGTTTCTCGGCGTTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCG  >  CP009273/3175902‑3176199

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: