Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,810,195 C→A 100% R54S (CGT→AGT)  yicC → UPF0701 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,810,1950CA99.5% 594.2 / ‑4.2 181R54S (CGT→AGT) yicCUPF0701 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (80/100);  total (81/100)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 181 total aligned reads displayed.

ACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTG  >  CP009273/3810051‑3810333
                                                                                                                                                |                                                                                                                                          
aCCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGtt                                                                                                                                         >  1:2123894/1‑148 (MQ=255)
aCCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGtt                                                                                                                                         <  2:830625/148‑1 (MQ=255)
 ccGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTtc                                                                                                                                        <  2:2395075/148‑1 (MQ=255)
  cGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTtct                                                                                                                                       <  1:2712915/148‑1 (MQ=255)
     cTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGc                                                                                                                                    <  2:2803506/148‑1 (MQ=255)
        cGCCCGGCGTGAAATCCAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTg                                                                                                                                 <  2:2679785/148‑1 (MQ=255)
        cGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTg                                                                                                                                 <  1:3029439/148‑1 (MQ=255)
           ccGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGAcg                                                                                                                              <  2:896719/148‑1 (MQ=255)
                 tGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTNTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGt                                                                                                                        <  1:1811056/148‑1 (MQ=255)
                  gAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTa                                                                                                                       >  2:2442352/1‑148 (MQ=255)
                  gAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTa                                                                                                                       >  2:1853292/1‑148 (MQ=255)
                  gAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTa                                                                                                                       <  2:910722/148‑1 (MQ=255)
                   aaaTCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGAcg                                                                                                                              >  1:2452419/1‑140 (MQ=255)
                     aTCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAg                                                                                                                    >  1:1295704/1‑148 (MQ=255)
                      tCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGt                                                                                                                   >  2:1198434/1‑148 (MQ=255)
                       cAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTg                                                                                                                  >  2:3395640/1‑148 (MQ=255)
                        aaGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTgg                                                                                                                 >  1:2595695/1‑148 (MQ=255)
                           ggTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAAt                                                                                                              >  1:1482385/1‑148 (MQ=255)
                           ggTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAAt                                                                                                              >  1:1482351/1‑148 (MQ=255)
                             tGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGt                                                                                                            <  2:2243179/148‑1 (MQ=255)
                             tGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGt                                                                                                            <  2:450096/148‑1 (MQ=255)
                              gAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTa                                                                                                           <  2:2965770/148‑1 (MQ=255)
                               aaTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTAc                                                                                                          <  2:3525686/148‑1 (MQ=255)
                               aaTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTAc                                                                                                          >  2:2247949/1‑148 (MQ=255)
                               aaTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTAc                                                                                                          <  2:3185859/148‑1 (MQ=255)
                                 tGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTAccc                                                                                                        <  1:3418468/148‑1 (MQ=255)
                                         cgcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTAt                                                                                                >  2:1979914/1‑148 (MQ=255)
                                         cgcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTAt                                                                                                >  2:179240/1‑148 (MQ=255)
                                          gcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATg                                                                                               <  2:2456075/148‑1 (MQ=255)
                                             aCCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGc                                                                                            >  1:1026082/1‑148 (MQ=255)
                                              ccTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGcc                                                                                           >  2:1802586/1‑148 (MQ=255)
                                              ccTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGcc                                                                                           <  2:1491899/148‑1 (MQ=255)
                                              ccTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGcc                                                                                           <  2:739859/148‑1 (MQ=255)
                                              ccTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCAGCGCTATGAGcc                                                                                           >  1:3289725/1‑148 (MQ=255)
                                               cTGTGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCa                                                                                          >  1:3281200/1‑148 (MQ=255)
                                                tGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAg                                                                                         <  1:802893/148‑1 (MQ=255)
                                                tGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAg                                                                                         <  1:1849637/148‑1 (MQ=255)
                                                 gggAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGa                                                                                        <  2:1635526/148‑1 (MQ=255)
                                                  ggAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGGGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGAt                                                                                       <  2:1937241/148‑1 (MQ=255)
                                                  ggAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGAt                                                                                       <  1:3070744/148‑1 (MQ=255)
                                                   gAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGAAATGTACCCTGCGCTATGAGCCAGATg                                                                                      <  2:655312/148‑1 (MQ=255)
                                                       tGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAg                                                                                  >  1:2119606/1‑148 (MQ=255)
                                                         cgcTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcg                                                                                <  1:2271209/148‑1 (MQ=255)
                                                            tCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCGGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcgcgc                                                                             >  1:1581048/1‑148 (MQ=255)
                                                             cGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGACAGTATTAGTTCTCGCCTGACGCTCGGTAAAGTGGAATGTACCCTCCGCTATGAGCCCGATGTTAGCGCGCa                                                                            >  1:3079264/1‑148 (MQ=255)
                                                               gTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAg                                                                          >  1:948699/1‑148 (MQ=255)
                                                               gTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAg                                                                          >  1:950021/1‑148 (MQ=255)
                                                                   aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGa                                                                      >  2:2151673/1‑148 (MQ=255)
                                                                   aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGa                                                                      >  2:926677/1‑148 (MQ=255)
                                                                   aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGa                                                                      >  2:926751/1‑148 (MQ=255)
                                                                    ccAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAg                                                                     <  2:597445/148‑1 (MQ=255)
                                                                    ccAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTATCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAg                                                                     <  1:350784/148‑1 (MQ=255)
                                                                    ccAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTATCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAg                                                                     <  1:648854/148‑1 (MQ=255)
                                                                    ccAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTATCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAg                                                                     <  1:351970/148‑1 (MQ=255)
                                                                      aGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCt                                                                   >  1:1869145/1‑148 (MQ=255)
                                                                        cGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGa                                                                 >  2:1393562/1‑148 (MQ=255)
                                                                        cGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGGTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGa                                                                 <  1:2847127/148‑1 (MQ=255)
                                                                         gTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGAt                                                                <  2:1295704/148‑1 (MQ=255)
                                                                           tATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATcc                                                              >  2:2597984/1‑148 (MQ=255)
                                                                            aTCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCt                                                             <  2:3246366/148‑1 (MQ=255)
                                                                            aTCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCt                                                             <  1:1198434/148‑1 (MQ=255)
                                                                             tCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGGAAAGTGGAATGTACCCTGCGCTCTGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTc                                                            >  2:715223/1‑148 (MQ=255)
                                                                                 gAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACg                                                        <  1:932700/148‑1 (MQ=255)
                                                                                  aaaCTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTCATATTTATCGCCTGACGCGCGGTAAAGAGGAAAGTACCCTGCGCTATGATCCAGATGTTATCGCGCAAGGAGATCTGAAACTCAACGa                                                       >  1:1739555/1‑148 (MQ=255)
                                                                                  aaaCTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGa                                                       >  2:2191368/1‑148 (MQ=255)
                                                                                    acttacttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGaaa                                                     <  2:2309165/148‑1 (MQ=255)
                                                                                      ttactgTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGaaaaa                                                   <  2:2465350/148‑1 (MQ=255)
                                                                                       tacttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAAc                                                  >  2:1027404/1‑148 (MQ=255)
                                                                                        acttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACt                                                 >  2:371164/1‑148 (MQ=255)
                                                                                           ttCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGc                                              <  1:1337913/148‑1 (MQ=255)
                                                                                           ttCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGc                                              <  2:2595695/148‑1 (MQ=255)
                                                                                               tCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTaaa                                          <  2:1949015/148‑1 (MQ=255)
                                                                                                  gCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAg                                       >  2:460007/1‑148 (MQ=255)
                                                                                                      gAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTgg                                   <  1:2514339/148‑1 (MQ=255)
                                                                                                      gAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTgg                                   <  2:530310/148‑1 (MQ=255)
                                                                                                      gAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTgg                                   <  1:2514322/148‑1 (MQ=255)
                                                                                                      gAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTgg                                   >  1:2156690/1‑148 (MQ=255)
                                                                                                            ttCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTg                             >  2:2593184/1‑148 (MQ=255)
                                                                                                               cGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcg                          <  1:1985687/148‑1 (MQ=255)
                                                                                                                gTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgc                         >  2:997357/1‑148 (MQ=255)
                                                                                                                gTAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgc                         >  2:3478309/1‑148 (MQ=255)
                                                                                                                 tAGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgcg                        >  1:1429597/1‑148 (MQ=255)
                                                                                                                  aGCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGa                       <  1:698302/148‑1 (MQ=255)
                                                                                                                   gCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGaa                      >  2:2537159/1‑148 (MQ=255)
                                                                                                                   gCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGaa                      <  2:1998853/148‑1 (MQ=255)
                                                                                                                   gCCTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGaa                      <  1:2316985/148‑1 (MQ=255)
                                                                                                                    ccTTGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAAc                     <  2:2637233/148‑1 (MQ=255)
                                                                                                                      ttGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTg                   <  1:2028713/148‑1 (MQ=255)
                                                                                                                      ttGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTg                   >  1:3135190/1‑148 (MQ=255)
                                                                                                                       tGAACCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCNAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTgg                  <  1:1807731/148‑1 (MQ=255)
                                                                                                                          aCCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTa               >  2:1960017/1‑148 (MQ=255)
                                                                                                                          aCCTGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTa               >  2:1960003/1‑148 (MQ=255)
                                                                                                                             tGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTcaaa            >  2:2560912/1‑148 (MQ=255)
                                                                                                                             tGTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaa            <  2:2898734/148‑1 (MQ=255)
                                                                                                                              gTCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaaa           <  1:3403339/148‑1 (MQ=255)
                                                                                                                               tCGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAAt          >  2:1556037/1‑148 (MQ=255)
                                                                                                                                cGTTCGCGAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATg         >  2:3382440/1‑148 (MQ=255)
                                                                                                                                    cgcgAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCaga     <  1:2966430/148‑1 (MQ=255)
                                                                                                                                    cgcgAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCaga     <  1:3339235/148‑1 (MQ=255)
                                                                                                                                       gAGCGTATTAGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTg  >  2:1604601/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                          
ACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTG  >  CP009273/3810051‑3810333

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: