Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,736,714 C→A 100% R325S (CGC→AGC)  yfiN → putative membrane‑anchored diguanylate cyclase with an N‑terminal periplasmic domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,736,7140CA96.8% 675.4 / 11.7 218R325S (CGC→AGC) yfiNputative membrane‑anchored diguanylate cyclase with an N‑terminal periplasmic domain
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (128/83);  total (130/88)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 218 total aligned reads displayed.

CCCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCA  >  CP009273/2736568‑2736858
                                                                                                                                                  |                                                                                                                                                
cccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCg                                                                                                                                                 <  2:3566969/148‑1 (MQ=255)
 ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGc                                                                                                                                                >  2:2091372/1‑148 (MQ=255)
 ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGc                                                                                                                                                <  1:2014388/148‑1 (MQ=255)
  cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGcc                                                                                                                                               <  2:2642554/148‑1 (MQ=255)
  cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGcc                                                                                                                                               <  2:3449242/148‑1 (MQ=255)
  cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATAAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGcc                                                                                                                                               >  1:3551647/1‑148 (MQ=255)
   gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                              <  2:689230/148‑1 (MQ=255)
   gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                              <  2:248868/148‑1 (MQ=255)
       aaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGggc                                                                                                                                          >  2:885391/1‑148 (MQ=255)
         cGTCGGCGTTACTATTTCTTTATGGCGAAAATTTCAAATAAATCAAAGATACCTTGTGTCATGCGACGGGCGATAGAGTCTTTAATGAAATAGCAAAACGGTTAGCTCAATTTGGCGGGCTGCGACATAAAGCATACAGCCAGgtcgg                                                                                                                                        >  2:815477/1‑148 (MQ=255)
         cGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGggcgg                                                                                                                                        >  2:2784212/1‑148 (MQ=255)
              gCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATg                                                                                                                                   >  2:519590/1‑148 (MQ=255)
                 ttACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAt                                                                                                                                <  1:582130/148‑1 (MQ=255)
                    cTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCg                                                                                                                             >  2:2616289/1‑148 (MQ=255)
                     tATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGc                                                                                                                            >  1:633588/1‑148 (MQ=255)
                      aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                           >  1:1948104/1‑148 (MQ=255)
                      aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                           <  1:1008563/148‑1 (MQ=255)
                        ttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTAt                                                                                                                         >  1:3391460/1‑148 (MQ=255)
                            tggtGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTg                                                                                                                     >  1:2845014/4‑148 (MQ=255)
                              aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGct                                                                                                                   <  2:225929/148‑1 (MQ=255)
                              aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGct                                                                                                                   >  1:975199/1‑148 (MQ=255)
                              aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGct                                                                                                                   >  2:1358840/1‑148 (MQ=255)
                                 gCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTa                                                                                                                >  1:1477861/1‑148 (MQ=255)
                                  cGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTat                                                                                                               >  1:1167662/1‑148 (MQ=255)
                                  cGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTat                                                                                                               >  2:3107280/1‑148 (MQ=255)
                                     tAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTatgat                                                                                                            >  1:1876166/1‑148 (MQ=255)
                                        tttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTa                                                                                                         <  1:1324577/148‑1 (MQ=255)
                                         ttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTAc                                                                                                        >  2:290008/1‑148 (MQ=255)
                                         ttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTAc                                                                                                        >  2:1403250/1‑148 (MQ=255)
                                           cAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAg                                                                                                      <  2:1698654/148‑1 (MQ=255)
                                            aaaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGt                                                                                                     <  1:3139507/148‑1 (MQ=255)
                                             aaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTc                                                                                                    >  2:2832561/1‑148 (MQ=255)
                                              aTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCg                                                                                                   <  2:668613/148‑1 (MQ=255)
                                               tACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCgg                                                                                                  >  1:482310/1‑148 (MQ=255)
                                               tACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCgg                                                                                                  >  1:482786/1‑148 (MQ=255)
                                                 cATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGaa                                                                                                >  1:2510015/1‑148 (MQ=255)
                                                   tCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATc                                                                                              >  1:2964145/1‑148 (MQ=255)
                                                    cAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCt                                                                                             >  2:1652102/1‑148 (MQ=255)
                                                      aTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                           <  1:289663/148‑1 (MQ=255)
                                                       tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  2:923508/1‑148 (MQ=255)
                                                       tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  2:2980696/1‑148 (MQ=255)
                                                           aCCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                      >  2:1300747/1‑148 (MQ=255)
                                                              tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgcagc                                                                                   >  2:819849/1‑148 (MQ=255)
                                                                gggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                 <  2:632232/148‑1 (MQ=255)
                                                                gggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                 <  2:635231/148‑1 (MQ=255)
                                                                gggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                 >  2:3487696/1‑148 (MQ=255)
                                                                  gTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGat                                                                               <  1:318397/148‑1 (MQ=255)
                                                                   tCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGata                                                                              >  1:2056296/1‑148 (MQ=255)
                                                                   tCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGata                                                                              >  2:2234311/1‑148 (MQ=255)
                                                                   tCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGata                                                                              >  1:2388714/1‑148 (MQ=255)
                                                                     aTGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTCTGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATg                                                                            >  2:3264124/1‑148 (MQ=255)
                                                                     aTGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATg                                                                            >  2:2520288/1‑148 (MQ=255)
                                                                      tGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGt                                                                           >  2:1209685/1‑148 (MQ=255)
                                                                      tGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGt                                                                           <  1:2625215/148‑1 (MQ=255)
                                                                       gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGtt                                                                          >  1:2488816/1‑148 (MQ=255)
                                                                           cGGGCGATAGAGTCTTGATTGAAATCGCAAACCGGTTAGCTGAATTTGGCGGGCTGCGTCATTAAGCATACAGCCTGGGCGTCGATTTATTCGATATGGTGGTCTATGTTGTACTGTCGGAATCTGAAGTGCAGCAGTTATGTTCAGc                                                                      <  1:3052268/148‑1 (MQ=255)
                                                                            gggCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCa                                                                     <  1:2524364/148‑1 (MQ=255)
                                                                                   agagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGGGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGacac                                                              >  2:899540/1‑148 (MQ=255)
                                                                                    gagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGAcaca                                                             >  1:2652954/1‑148 (MQ=255)
                                                                                     agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:382571/148‑1 (MQ=255)
                                                                                     agTCTTGATTGAAATCGCAAAACGGTTACCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:1487537/148‑1 (MQ=255)
                                                                                     agTCTTGATTGAAATCGCAAAACGGTGAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  1:431469/148‑1 (MQ=255)
                                                                                           gattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCttt                                                      <  2:2277070/148‑1 (MQ=255)
                                                                                              tgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                   >  1:2588412/1‑148 (MQ=255)
                                                                                               gaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAtc                                                  >  2:895626/1‑148 (MQ=255)
                                                                                                  aTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTcc                                               >  2:2684989/1‑148 (MQ=255)
                                                                                                   tCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTccc                                              >  2:1150430/1‑148 (MQ=255)
                                                                                                   tCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTccc                                              >  1:1799590/1‑148 (MQ=255)
                                                                                                   tCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTccc                                              >  2:1723357/1‑148 (MQ=255)
                                                                                                   tCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTccc                                              >  2:1150319/1‑148 (MQ=255)
                                                                                                      cAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGtt                                           >  1:1110669/1‑148 (MQ=255)
                                                                                                        aaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                         >  1:1617739/1‑148 (MQ=255)
                                                                                                         aaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGa                                        <  1:52177/148‑1 (MQ=255)
                                                                                                             gTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCtt                                    >  1:3450750/1‑148 (MQ=255)
                                                                                                              ttAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                   >  2:2484859/1‑148 (MQ=255)
                                                                                                              ttAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                   >  2:34960/1‑148 (MQ=255)
                                                                                                              ttAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                   <  1:836921/148‑1 (MQ=255)
                                                                                                               tAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCa                                  >  2:1895252/1‑148 (MQ=255)
                                                                                                               tAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCa                                  >  2:2593169/1‑148 (MQ=255)
                                                                                                                aGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCAt                                 <  1:916005/148‑1 (MQ=255)
                                                                                                                 gCTGAATTTGGCGGGCTGCGACATAAAGCATACATCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATa                                <  2:1095011/148‑1 (MQ=255)
                                                                                                                  cTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATaa                               <  2:1948104/148‑1 (MQ=255)
                                                                                                                     aaTTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATgg                            >  2:111597/1‑148 (MQ=255)
                                                                                                                     aaTTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATgg                            <  1:1535623/148‑1 (MQ=255)
                                                                                                                       tttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGcc                          >  1:484134/1‑148 (MQ=255)
                                                                                                                        ttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCa                         >  2:3378526/1‑148 (MQ=255)
                                                                                                                        ttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCa                         >  2:1610625/1‑148 (MQ=255)
                                                                                                                        ttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCa                         >  1:779304/1‑148 (MQ=255)
                                                                                                                         tGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCAt                        <  2:3341286/148‑1 (MQ=255)
                                                                                                                             gggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGa                    <  1:2613195/148‑1 (MQ=255)
                                                                                                                                cTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGacca                 <  1:1358840/148‑1 (MQ=255)
                                                                                                                                     aCATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATg            >  1:3329367/1‑148 (MQ=255)
                                                                                                                                       aTAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGAc          >  2:1265313/1‑148 (MQ=255)
                                                                                                                                       aTAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGAc          <  2:1516242/148‑1 (MQ=255)
                                                                                                                                         aaaGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAt        >  1:1202982/1‑148 (MQ=255)
                                                                                                                                         aaaGCATACATCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAt        <  2:975199/148‑1 (MQ=255)
                                                                                                                                             cATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAg    >  1:2153036/1‑148 (MQ=255)
                                                                                                                                             cATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAg    >  1:2119281/1‑148 (MQ=255)
                                                                                                                                               tACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCa  >  1:2765328/1‑148 (MQ=255)
                                                                                                                                               tACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCa  <  2:1759395/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCA  >  CP009273/2736568‑2736858

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: