Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,265,415 G→T 15.9% R6L (CGA→CTA)  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,265,4150GT15.9% 763.4 / 84.2 238R6L (CGA→CTA) yjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (113/87);  new base T (16/22);  total (129/109)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 243 total aligned reads displayed.

TGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGG  >  CP009273/4265266‑4265564
                                                                                                                                                     |                                                                                                                                                     
tGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACg                                                                                                                                                       <  2:2467490/150‑1 (MQ=255)
  aCGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGac                                                                                                                                                     >  2:297100/1‑150 (MQ=255)
    gATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTacac                                                                                                                                                   >  1:1098725/1‑150 (MQ=255)
      ttAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCa                                                                                                                                                 >  1:1442184/1‑150 (MQ=255)
        aGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAAt                                                                                                                                               >  2:1572041/1‑150 (MQ=255)
         gAATGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAAtt                                                                                                                                              <  2:4322071/150‑1 (MQ=255)
          aaaGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTa                                                                                                                                             <  1:4068244/150‑1 (MQ=255)
          aaaGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGACATCGTGCACGACACCAATTa                                                                                                                                             >  2:2224387/1‑150 (MQ=255)
             ggTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTg                                                                                                                                          <  1:1859016/150‑1 (MQ=255)
             ggTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTg                                                                                                                                          <  2:1268294/150‑1 (MQ=255)
               tATCAATTTCAAATCAGGCAAAAGTGATATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGc                                                                                                                                        <  1:1166429/150‑1 (MQ=255)
                aTCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCg                                                                                                                                       >  2:1929453/1‑150 (MQ=255)
                   aaTTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAGTGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTg                                                                                                                                    >  2:3234848/1‑150 (MQ=255)
                      ttCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCg                                                                                                                                 >  2:2806397/1‑150 (MQ=255)
                       tCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCgg                                                                                                                                >  2:1831987/1‑150 (MQ=255)
                         aaaTCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGc                                                                                                                              <  1:3901698/150‑1 (MQ=255)
                          aaTCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCa                                                                                                                             >  1:2025855/1‑150 (MQ=255)
                          aaTCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCa                                                                                                                             >  1:1729505/1‑150 (MQ=255)
                           aTCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAt                                                                                                                            <  2:1484962/150‑1 (MQ=255)
                             cAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTa                                                                                                                          >  2:1027176/1‑150 (MQ=255)
                                  aaaaGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCttt                                                                                                                     >  2:521406/1‑150 (MQ=255)
                                  aaaaGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCttt                                                                                                                     >  2:1352155/1‑150 (MQ=255)
                                  aaaaGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCttt                                                                                                                     >  1:3461248/1‑150 (MQ=255)
                                    aaGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCttttt                                                                                                                   <  1:2329948/150‑1 (MQ=255)
                                       tGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGt                                                                                                                <  2:3531869/150‑1 (MQ=255)
                                        gCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGtt                                                                                                               <  2:3809008/150‑1 (MQ=255)
                                        gCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGtt                                                                                                               >  1:1493829/1‑150 (MQ=255)
                                        gCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGtt                                                                                                               >  1:160314/1‑150 (MQ=255)
                                           aTTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCCTTATCTTTTTAGTtctc                                                                                                            >  2:4027662/1‑150 (MQ=255)
                                            tttATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTctct                                                                                                           <  2:4347573/150‑1 (MQ=255)
                                            tttATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTctct                                                                                                           <  2:3731836/150‑1 (MQ=255)
                                             ttATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCtt                                                                                                          <  2:3619659/150‑1 (MQ=255)
                                                taCCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTcc                                                                                                       >  2:2477623/1‑150 (MQ=255)
                                                       aGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCAtt                                                                                                <  1:4107084/150‑1 (MQ=255)
                                                        gATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTc                                                                                               >  1:3301764/1‑150 (MQ=255)
                                                           ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCttt                                                                                            <  1:4177636/150‑1 (MQ=255)
                                                           ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCttt                                                                                            <  1:4619841/150‑1 (MQ=255)
                                                            tATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTc                                                                                           <  1:1925152/150‑1 (MQ=255)
                                                             aTCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCg                                                                                          >  2:2259663/1‑150 (MQ=255)
                                                                tAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTa                                                                                       >  1:2149686/1‑150 (MQ=255)
                                                                 aaaGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTAt                                                                                      >  2:4270952/1‑150 (MQ=255)
                                                                 aaaGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTAt                                                                                      >  2:26443/1‑150 (MQ=255)
                                                                  aaGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATg                                                                                     >  1:944703/1‑150 (MQ=255)
                                                                   aGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATgt                                                                                    >  2:3370984/1‑149 (MQ=255)
                                                                      cGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGAtt                                                                                 >  2:2296701/1‑150 (MQ=255)
                                                                      cGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGAtt                                                                                 >  2:2297062/1‑150 (MQ=255)
                                                                      cGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGAtt                                                                                 >  1:3239557/1‑150 (MQ=255)
                                                                       gTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTg                                                                                <  2:455413/150‑1 (MQ=255)
                                                                             aCCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGccttcc                                                                          <  2:3743448/150‑1 (MQ=255)
                                                                              cccAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTCCTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGccttcct                                                                         <  1:647545/150‑1 (MQ=255)
                                                                                cAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCttt                                                                       >  1:3515096/1‑150 (MQ=255)
                                                                                cAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCttt                                                                       <  1:3105977/150‑1 (MQ=255)
                                                                                cAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCttt                                                                       <  1:3679334/150‑1 (MQ=255)
                                                                                  gggTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTgg                                                                     <  1:1555946/150‑1 (MQ=255)
                                                                                  gggTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTgg                                                                     <  2:3286620/150‑1 (MQ=255)
                                                                                   ggTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTggg                                                                    <  2:621422/150‑1 (MQ=255)
                                                                                   ggTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTggg                                                                    >  1:3277067/1‑150 (MQ=255)
                                                                                         cACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGAGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGATATGGATTGCCTTCCTTTGGGCAAAAt                                                              >  2:2742108/1‑150 (MQ=255)
                                                                                           ccTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCa                                                            <  2:481091/150‑1 (MQ=255)
                                                                                              tGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGaa                                                         >  2:1418986/1‑150 (MQ=255)
                                                                                                 aaTGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTg                                                      >  2:2657394/1‑150 (MQ=255)
                                                                                                 aaTGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTg                                                      >  1:3944368/1‑150 (MQ=255)
                                                                                                  aTGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGa                                                     <  1:79797/150‑1 (MQ=255)
                                                                                                    ggCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGaat                                                   >  2:3787754/1‑150 (MQ=255)
                                                                                                    ggCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGaat                                                   <  1:3977799/150‑1 (MQ=255)
                                                                                                     gCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGaata                                                  <  1:2749836/150‑1 (MQ=255)
                                                                                                        gggTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATc                                               <  1:2251676/150‑1 (MQ=255)
                                                                                                        gggTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATc                                               >  2:449386/1‑150 (MQ=255)
                                                                                                          gTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAg                                             >  1:1212737/1‑150 (MQ=255)
                                                                                                          gTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAg                                             >  2:2910210/1‑150 (MQ=255)
                                                                                                           tataAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGc                                            <  1:2616038/150‑1 (MQ=255)
                                                                                                           tataAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGc                                            >  2:3880922/1‑150 (MQ=255)
                                                                                                           tataAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGc                                            <  1:4317461/150‑1 (MQ=255)
                                                                                                            ataAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCt                                           >  1:4228195/1‑150 (MQ=255)
                                                                                                                 caggcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGaa                                      <  1:490086/150‑1 (MQ=255)
                                                                                                                  aggcaggAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAAc                                     <  1:3471487/150‑1 (MQ=255)
                                                                                                                  aggcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAAc                                     >  2:1634141/1‑150 (MQ=255)
                                                                                                                    gcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCt                                   >  1:4218940/1‑150 (MQ=255)
                                                                                                                     caggAAATTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCtt                                  >  2:4090215/1‑150 (MQ=255)
                                                                                                                     caggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTGTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCtt                                  >  2:353916/1‑150 (MQ=255)
                                                                                                                       ggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTg                                >  2:4576467/1‑150 (MQ=255)
                                                                                                                          aaTTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTc                             <  2:2099858/150‑1 (MQ=255)
                                                                                                                          aaTTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTc                             <  2:2099748/150‑1 (MQ=255)
                                                                                                                           aTTGATAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCa                            <  2:1274183/150‑1 (MQ=255)
                                                                                                                           aTTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCa                            >  2:3971071/1‑150 (MQ=255)
                                                                                                                              gATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACt                         <  1:2997530/150‑1 (MQ=255)
                                                                                                                               aTAGCAATGAGTCATCGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTg                        <  2:2324728/150‑1 (MQ=255)
                                                                                                                                tAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTgg                       <  2:3991593/150‑1 (MQ=255)
                                                                                                                                tAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTgg                       >  2:805679/1‑150 (MQ=255)
                                                                                                                                         gTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATa              <  1:2353741/150‑1 (MQ=255)
                                                                                                                                          tCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATaa             >  2:1713128/1‑150 (MQ=255)
                                                                                                                                          tCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATaa             >  2:1138469/1‑150 (MQ=255)
                                                                                                                                           cATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATaaa            >  1:2769070/1‑150 (MQ=255)
                                                                                                                                           cATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATaaa            >  1:2435744/1‑150 (MQ=255)
                                                                                                                                              cGTGCACTACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATcc         <  1:20175/150‑1 (MQ=255)
                                                                                                                                              cGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATcc         >  2:441374/1‑150 (MQ=255)
                                                                                                                                                tGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGa       >  1:2603787/1‑150 (MQ=255)
                                                                                                                                                 gCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAg      <  1:3633501/150‑1 (MQ=255)
                                                                                                                                                   aCGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCt    >  1:4502974/1‑150 (MQ=255)
                                                                                                                                                     gACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTgg  >  1:935488/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGG  >  CP009273/4265266‑4265564

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: