Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,266,916 G→A 4.1% A506A (GCG→GCA yjcC → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,266,9160GA 4.1% 301.9 / 2.0 99A506A (GCG→GCAyjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (48/46);  new base A (2/2);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGA  >  CP009273/4266770‑4267064
                                                                                                                                                  |                                                                                                                                                    
aGTCATTTGATTGCGGAACACATCATCGGGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGAt                                                                                                                                                    <  2:27649/149‑1 (MQ=255)
  gcATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGGAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCCGGGATGGGTCTTTGCGAAGGCGATgc                                                                                                                                                  >  1:669453/2‑149 (MQ=255)
        gATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCAATGCcgccgc                                                                                                                                            >  1:1194580/1‑149 (MQ=255)
        gATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCAATGCcgccgc                                                                                                                                            >  1:251931/1‑149 (MQ=255)
           tGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAAAAGGTGATCAGGTTAACTGGCTGCTCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAgg                                                                                                                                         <  1:907990/149‑1 (MQ=255)
             cGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGGTCTTTGCGAAGGCGATGCCGCCGCAGgtg                                                                                                                                       >  1:101858/1‑149 (MQ=255)
               gAACCCATCATCGAGCCGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGtt                                                                                                                                     <  1:1141480/149‑1 (MQ=255)
               gAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCGAACGCGGCGTGCGCTATGGCCAGGGATGGTTCTTTGCGAAGGCGATGCGGCCGCAGgtgta                                                                                                                                     >  1:64258/1‑148 (MQ=255)
                  caccgcaccGAACTGACGCACAGCCTGGGGTAAAAAACGAACGCTGAAGCCGTCGAAACTGAGGAGCAGGTTACCTCGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTCTGCGAAGGCGATGCCGCCGCAGGTGTTTAt                                                                                                                                  <  1:436429/141‑1 (MQ=255)
                    catcatCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCCGGGATGATTCTTTGCGCAGGCGATGACGCGGCACGTGTTTATGc                                                                                                                                >  2:653550/1‑149 (MQ=255)
                     atcatcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCa                                                                                                                               >  1:228417/1‑149 (MQ=255)
                     atcatcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCa                                                                                                                               >  2:425074/1‑149 (MQ=255)
                     atcatcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCa                                                                                                                               >  1:507621/1‑149 (MQ=255)
                     atcatcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGGTTATGCa                                                                                                                               >  1:452838/1‑149 (MQ=255)
                      tcatcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCaa                                                                                                                              <  2:1077603/149‑1 (MQ=255)
                        ttcGGGCTGGCGCACAGCCTGGGGTTAAAAACGATCNCTGAAGGCGTCGAAACCGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTTCGCTATTGCNAGGGATGGTTCTTTGCGAAGGCGATGCCGCCTCAGGTGTTTATGCAatg                                                                                                                            <  2:916542/148‑1 (MQ=255)
                        atcGAGCTGGCGCAGAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAatg                                                                                                                            <  2:383306/149‑1 (MQ=255)
                        atcGAGCTGGCGCACTGCCTGGGGTTAACAACGATCGGTGCAGGCGTCACAACTAGGGAGCGGGACCACTCGCCGCCCAAACTCTGCGGGCGCTATTGCCAGGGATTGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTAGTCAatg                                                                                                                            <  1:555515/149‑1 (MQ=255)
                         tcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTGTGCCatgg                                                                                                                           >  2:693875/1‑149 (MQ=255)
                         tcGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAatgg                                                                                                                           >  1:74108/1‑149 (MQ=255)
                          cGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAATCGCTGCTTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAatgga                                                                                                                          <  1:1095098/149‑1 (MQ=255)
                            aGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGAGGCCGTCGCAGGTGTTTATGCAatggatg                                                                                                                        <  2:1194348/149‑1 (MQ=255)
                                ggCGCACAGCTTGGGGTTAAAAACGATTGCGGAAGGCGGCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGc                                                                                                                    <  1:128691/149‑1 (MQ=255)
                                  cgcACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCaa                                                                                                                  <  1:779838/149‑1 (MQ=255)
                                   gcACAGCCTGGGGTTAATAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTTGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATCGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAAt                                                                                                                 >  2:916045/1‑149 (MQ=255)
                                       aGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGATCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTAcc                                                                                                             >  1:952532/1‑149 (MQ=255)
                                         ccTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGGTGCCGCCGCAGGTGTTTATGCAATGGATGGAGCGATGAcctg                                                                                                           >  1:40012/1‑147 (MQ=255)
                                           tGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGCAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCcgcg                                                                                                         >  2:976457/1‑149 (MQ=255)
                                           tGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCcgcg                                                                                                         >  2:1168386/1‑149 (MQ=255)
                                           tGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCcgcg                                                                                                         >  2:541828/1‑149 (MQ=255)
                                           tGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCcgcg                                                                                                         >  1:631825/1‑149 (MQ=255)
                                              ggTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCgg                                                                                                      >  2:428965/1‑149 (MQ=255)
                                              ggTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGGGTTTATGCAATGGATGGAGCAATTACCCGCGCgg                                                                                                      >  2:429222/1‑149 (MQ=255)
                                                ttAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCgggg                                                                                                    >  1:504200/1‑148 (MQ=255)
                                                 tAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAg                                                                                                   <  1:295213/149‑1 (MQ=255)
                                                  aaataCGATCGCCGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGt                                                                                                  <  2:1019884/149‑1 (MQ=255)
                                                    aaaCGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGGTGGAGCAATTACCCGCGCGGGAGTTa                                                                                                >  2:1148562/1‑149 (MQ=255)
                                                    aaaCGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCGGCGCGGGAGTTa                                                                                                >  2:260082/1‑149 (MQ=255)
                                                        gATCGCTGAAGGCGTCGAAACTGAGTAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTCGCGAACGCGATGCTGCCGCAGGTGTTTATGCAATGGACGGAGCAATTACCCGCGCGGGAGTTAAcgc                                                                                            <  1:758391/149‑1 (MQ=255)
                                                        gATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCGGCAGGTGTTTATGCAATGGAGGGAGCAATTACCCGCGCGGGAGTTAAcgc                                                                                            >  1:419930/1‑149 (MQ=255)
                                                        gATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAAcgc                                                                                            >  1:54174/1‑149 (MQ=255)
                                                             cTGAAGGCGTCGGAACTGAGGAGGAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCggg                                                                                       <  2:474014/149‑1 (MQ=255)
                                                             cTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCggg                                                                                       <  1:465992/149‑1 (MQ=255)
                                                                  ggCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCTGGAGTTTACGCGCGGGCAATa                                                                                  <  2:1070159/149‑1 (MQ=255)
                                                                  gcgggcGAAACGGAGGAGCAGGTTAACTGGGTGCGCAAACGCGGCGTGCGGTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATa                                                                                  <  2:85272/144‑1 (MQ=255)
                                                                       cGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAAtt                                                                             >  1:120037/1‑149 (MQ=255)
                                                                          aaCTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACa                                                                          >  1:789888/1‑149 (MQ=255)
                                                                            cTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAgg                                                                        <  2:1207155/149‑1 (MQ=255)
                                                                                 gagCAGGTTAACTGGCTGCGCAAACGCCGCGTGCGCTATTGCCAGGGATGGTTCTTCGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACACGCGCGGGAGCTAAAGCGCGGGCAATAACATTACAGGCGGTg                                                                   >  1:562461/1‑149 (MQ=255)
                                                                                  gtccggTTAACGCGCTGCGCGAACTCCGCCTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATGACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTgg                                                                  <  2:1208524/145‑1 (MQ=255)
                                                                                  agCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCAATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTgg                                                                  <  2:251931/149‑1 (MQ=255)
                                                                                  agCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCAATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTgg                                                                  <  2:1194580/149‑1 (MQ=255)
                                                                                  agCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCATTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGTCAATAAAATTACAGGCGGTgg                                                                  <  1:698725/149‑1 (MQ=255)
                                                                                    cAGGTTAATTGGCTGCGGAAACGCGGAGTGGTGTTTTGTCATGGTTGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCg                                                                <  1:647438/149‑1 (MQ=255)
                                                                                           aCTGGCTGCGCAAACGCGGCGGGCGCGCTTGCCAGGGATCGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCGATAACATTGCGGGCGGTGGCGATAATCg                                                         <  2:146867/149‑1 (MQ=255)
                                                                                            cTGGCTGCGGAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGc                                                        <  2:127779/149‑1 (MQ=255)
                                                                                             tGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCt                                                       >  1:71692/1‑149 (MQ=255)
                                                                                               gccgcgcAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATAGAGCAATTAGCCGCGCGGGAGTTCACGCGCGCGCAATAAAATTACAGGCGGTGGCGATAATCGCTgg                                                     <  1:237160/146‑1 (MQ=255)
                                                                                                  gcgcAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAg                                                  >  1:675801/1‑149 (MQ=255)
                                                                                                      aaaCGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATGGCTGGGGGGGCg                                              >  2:554064/1‑149 (MQ=255)
                                                                                                       caCGCGGCGTGCGCTATTGCCAGGGATGGTTCGTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGa                                             <  2:505755/148‑1 (MQ=255)
                                                                                                       aaCGCGGCGTGCCCCATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGa                                             <  1:137597/149‑1 (MQ=255)
                                                                                                          gcggcgTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCa                                          <  1:774602/149‑1 (MQ=255)
                                                                                                           cggcgTGCGCCATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCaa                                         <  1:554881/149‑1 (MQ=255)
                                                                                                            ggcgTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCaaa                                        <  2:782014/149‑1 (MQ=255)
                                                                                                                tgcgcagctgcgatggAGGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGGTCAAACTGc                                    >  2:12630/15‑149 (MQ=255)
                                                                                                                tGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGc                                    >  1:1155181/1‑149 (MQ=255)
                                                                                                                tGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGc                                    >  1:401571/1‑149 (MQ=255)
                                                                                                                tGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGc                                    <  2:507621/149‑1 (MQ=255)
                                                                                                                tGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCGGGGAGTGCGATCAAACTGc                                    >  1:868924/1‑149 (MQ=255)
                                                                                                                ggcgcTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGc                                    <  2:228417/148‑1 (MQ=255)
                                                                                                                   gcTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGa                                 <  1:563935/149‑1 (MQ=255)
                                                                                                                   gcTATTGCCAGGGATGGTCCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGa                                 <  1:446824/149‑1 (MQ=255)
                                                                                                                    cTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGAc                                >  2:171440/1‑149 (MQ=255)
                                                                                                                      aTTGCCAGGGATGGTTCTTTGAGACGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACgg                              >  2:450822/1‑149 (MQ=255)
                                                                                                                          ccAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCCATAAAATTACAGGCGGTGGCGATAATCGGTGGGGGGGCGGACAAACGGCCGACGGaaaa                          >  2:595318/1‑149 (MQ=255)
                                                                                                                          ccAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGTGTGCGATCAAACTGCCGACGGaaaa                          <  1:621216/149‑1 (MQ=255)
                                                                                                                          ccAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGaaaa                          >  1:763538/1‑149 (MQ=255)
                                                                                                                          ccAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGAAGTGCGATCAAACTGCCGACGGaaaa                          >  1:224329/1‑149 (MQ=255)
                                                                                                                            aGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTCTGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAAcg                        >  2:100976/1‑149 (MQ=255)
                                                                                                                             gggATGGTTTCTTGCCAAGGCGATGCCGCCGCAGGTGTCTATCCAATGGATGGAGCAATTACCCCCGCAGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAAcgc                       <  1:290280/149‑1 (MQ=255)
                                                                                                                              ggATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAAcgcg                      >  2:1154981/1‑149 (MQ=255)
                                                                                                                              ggATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAAcgcg                      >  2:919046/1‑149 (MQ=255)
                                                                                                                               gATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCgg                     <  2:775571/149‑1 (MQ=255)
                                                                                                                                aTGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAAATACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCggg                    >  2:175958/1‑149 (MQ=255)
                                                                                                                                 tgcttCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGga                   <  2:1048041/146‑1 (MQ=255)
                                                                                                                                 tGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGga                   <  2:1048194/149‑1 (MQ=255)
                                                                                                                                 tGGTTCTTTGCGAAGGCGAGGCCGCCGTAGGTGTTTCGGGCATGGAGGGAGCGATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGga                   <  2:1182935/149‑1 (MQ=255)
                                                                                                                                   gTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGAGCGATCAAACTGCCGACGGAAAACGGGGgaga                 >  1:43937/1‑149 (MQ=255)
                                                                                                                                     tCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAg               >  2:449897/1‑149 (MQ=255)
                                                                                                                                      cTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAgg              <  1:215102/149‑1 (MQ=255)
                                                                                                                                      cTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAgg              <  2:74108/149‑1 (MQ=255)
                                                                                                                                       tttGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTCCAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGt             >  1:300466/1‑149 (MQ=255)
                                                                                                                                           cGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGAGGGACAGCGCGGGAGAAGGTctgc         >  1:1163357/1‑149 (MQ=255)
                                                                                                                                           cGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGACAACGCGGGAGAAGGGctgc         >  2:789444/1‑149 (MQ=255)
                                                                                                                                             aaGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTctgctg       <  2:452838/149‑1 (MQ=255)
                                                                                                                                             aaGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTctgctg       <  1:425074/149‑1 (MQ=255)
                                                                                                                                               ggCCATGCCGCCGCAGGGGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCg     <  1:86766/149‑1 (MQ=255)
                                                                                                                                                gCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCga    <  1:324401/149‑1 (MQ=255)
                                                                                                                                                  gATGCCGCCGCAGGTGTTTATGCAATGGATGGGGCCATGTCCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCgaga  <  2:355042/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCTATTGCCAGGGATGGTTCTTTGCGAAGGCGATGCCGCCGCAGGTGTTTATGCAATGGATGGAGCAATTACCCGCGCGGGAGTTAACGCGCGGGCAATAAAATTACAGGCGGTGGCGATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGA  >  CP009273/4266770‑4267064

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: