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,7140CA99.0% 648.3 / ‑1.7 200R325S (CGC→AGC) yfiNputative membrane‑anchored diguanylate cyclase with an N‑terminal periplasmic domain
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (116/82);  total (118/82)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

GAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCA  >  CP009273/2736571‑2736858
                                                                                                                                               |                                                                                                                                                
gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                              >  2:2708167/1‑148 (MQ=255)
   aaaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTggg                                                                                                                                           >  2:3173344/1‑148 (MQ=255)
   aaaCGTCGGCGTTAATATTTATTGATGGCGATAATTTCAAATACAACAATGATACCGGGGGTCATGCGACGGGCGATAGAGTCTTGATTTAAATAGCAAAACGGTTAGCTGAATTTTGCGGGCTGCGACATAAAGAATAAAACCTggg                                                                                                                                           >  1:1767885/1‑148 (MQ=255)
          ggCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGAt                                                                                                                                    >  1:3047600/1‑148 (MQ=255)
            cGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGa                                                                                                                                  <  2:1116955/148‑1 (MQ=255)
               tACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAtt                                                                                                                               >  1:380396/1‑148 (MQ=255)
                aCTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTc                                                                                                                              >  2:580467/1‑148 (MQ=255)
                 cTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCg                                                                                                                             >  2:1379095/1‑148 (MQ=255)
                 cTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCg                                                                                                                             <  1:1698418/148‑1 (MQ=255)
                   aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                           <  2:2007599/148‑1 (MQ=255)
                    tttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTa                                                                                                                          <  2:2354471/148‑1 (MQ=255)
                          gATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGc                                                                                                                    >  1:3381739/1‑148 (MQ=255)
                           aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGct                                                                                                                   <  2:2124836/148‑1 (MQ=255)
                             ggCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGctct                                                                                                                 >  1:3700514/1‑148 (MQ=255)
                             ggCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGAATACAGCCTGGGCGGCGATGAATTCTCTATGGTGctct                                                                                                                 >  1:3001725/1‑148 (MQ=255)
                                    aTTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGt                                                                                                          >  1:1553820/1‑148 (MQ=255)
                                    aTTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGt                                                                                                          >  1:313905/1‑148 (MQ=255)
                                     tttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTa                                                                                                         >  2:3286166/1‑148 (MQ=255)
                                     tttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTa                                                                                                         >  2:713822/1‑148 (MQ=255)
                                      ttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTAc                                                                                                        >  2:1756928/1‑148 (MQ=255)
                                          aaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTc                                                                                                    <  1:580467/148‑1 (MQ=255)
                                            tACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCgg                                                                                                  >  1:1261826/1‑148 (MQ=255)
                                             aCATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                 <  2:1035314/148‑1 (MQ=255)
                                             aCATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                 >  2:798730/1‑148 (MQ=255)
                                              cATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGaa                                                                                                <  2:2025462/148‑1 (MQ=255)
                                               aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAAt                                                                                               >  2:2708289/1‑148 (MQ=255)
                                               aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAAt                                                                                               <  2:313905/148‑1 (MQ=255)
                                                tCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATc                                                                                              >  1:2693542/1‑148 (MQ=255)
                                                tCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATc                                                                                              >  1:3432678/1‑148 (MQ=255)
                                                tCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATc                                                                                              >  1:549780/1‑148 (MQ=255)
                                                   aTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                           <  1:2819049/148‑1 (MQ=255)
                                                   aTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                           >  1:1794105/1‑148 (MQ=255)
                                                    tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  1:326792/1‑148 (MQ=255)
                                                    tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  2:1640961/1‑148 (MQ=255)
                                                    tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  2:2247060/1‑148 (MQ=255)
                                                    tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                          >  2:3727322/1‑148 (MQ=255)
                                                     gATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAg                                                                                         >  1:2232067/1‑148 (MQ=255)
                                                      aTACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGt                                                                                        >  1:41643/1‑148 (MQ=255)
                                                       tACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTg                                                                                       >  2:1304782/1‑148 (MQ=255)
                                                        aCCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                      >  1:1254921/1‑148 (MQ=255)
                                                        aCCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                      <  1:1484991/148‑1 (MQ=255)
                                                           tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgcagc                                                                                   >  1:2566456/1‑148 (MQ=255)
                                                             gggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                 >  1:2587761/1‑148 (MQ=255)
                                                              ggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGa                                                                                <  1:1904622/148‑1 (MQ=255)
                                                               gTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGat                                                                               <  1:3466611/148‑1 (MQ=255)
                                                                tCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGata                                                                              <  2:2439224/148‑1 (MQ=255)
                                                                    gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGtt                                                                          <  1:1358001/148‑1 (MQ=255)
                                                                    gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGtt                                                                          >  2:2907869/1‑148 (MQ=255)
                                                                    gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGtt                                                                          >  2:1200002/1‑148 (MQ=255)
                                                                       aCGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAGTTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAg                                                                       >  2:3692191/1‑148 (MQ=255)
                                                                           gCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACt                                                                   <  2:1550463/148‑1 (MQ=255)
                                                                            cGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTg                                                                  <  2:1143355/148‑1 (MQ=255)
                                                                            cGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTg                                                                  <  1:1993435/148‑1 (MQ=255)
                                                                            cGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTg                                                                  <  2:170499/148‑1 (MQ=255)
                                                                               tAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGaca                                                               >  1:13240/1‑148 (MQ=255)
                                                                               tAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGaca                                                               >  1:14198/1‑148 (MQ=255)
                                                                                  agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:765655/148‑1 (MQ=255)
                                                                                  agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:255272/148‑1 (MQ=255)
                                                                                  agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:765748/148‑1 (MQ=255)
                                                                                  agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            <  2:253987/148‑1 (MQ=255)
                                                                                  agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                            >  1:971038/1‑148 (MQ=255)
                                                                                    tcgtgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAAt                                                          <  2:2676834/145‑1 (MQ=255)
                                                                                    tCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAAt                                                          >  2:1873543/1‑148 (MQ=255)
                                                                                    tCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCGCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAAt                                                          <  2:2261969/148‑1 (MQ=255)
                                                                                     cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATc                                                         >  2:3283236/1‑148 (MQ=255)
                                                                                     cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATc                                                         >  2:2173238/1‑148 (MQ=255)
                                                                                     cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATc                                                         >  2:3105093/1‑148 (MQ=255)
                                                                                     cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATc                                                         >  2:324642/1‑148 (MQ=255)
                                                                                       tgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                       >  1:165879/1‑148 (MQ=255)
                                                                                       tgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                       >  1:3287926/1‑148 (MQ=255)
                                                                                         attgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTa                                                     >  2:652914/1‑148 (MQ=255)
                                                                                         attgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTa                                                     >  2:2725037/1‑148 (MQ=255)
                                                                                           tgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                   <  2:3215869/148‑1 (MQ=255)
                                                                                                     aaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                         <  2:1803840/148‑1 (MQ=255)
                                                                                                     aaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGACGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                         >  1:657010/1‑148 (MQ=255)
                                                                                                       aCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGAt                                       >  2:1323863/1‑148 (MQ=255)
                                                                                                         ggTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCt                                     <  2:1191277/148‑1 (MQ=255)
                                                                                                           ttAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                   <  2:2569580/148‑1 (MQ=255)
                                                                                                              gCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATa                                <  2:2282935/148‑1 (MQ=255)
                                                                                                              gCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATa                                <  1:3632157/148‑1 (MQ=255)
                                                                                                               cTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATaa                               <  2:522657/148‑1 (MQ=255)
                                                                                                               cTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATaa                               >  1:3340538/1‑148 (MQ=255)
                                                                                                                 gAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATg                             <  1:928605/148‑1 (MQ=255)
                                                                                                                    tttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGcc                          >  1:2847924/1‑148 (MQ=255)
                                                                                                                       ggCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATc                       <  1:1229243/148‑1 (MQ=255)
                                                                                                                         cGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAg                     <  1:2380446/148‑1 (MQ=255)
                                                                                                                             cTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGacca                 <  1:50708/148‑1 (MQ=255)
                                                                                                                              tGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGaccac                >  2:1666038/1‑148 (MQ=255)
                                                                                                                               gCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGaccacc               <  1:1685893/148‑1 (MQ=255)
                                                                                                                               gCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGaccacc               >  2:2981869/1‑148 (MQ=255)
                                                                                                                               gCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGaccacc               >  2:1308747/1‑148 (MQ=255)
                                                                                                                                cGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGAccacca              >  1:2548440/1‑148 (MQ=255)
                                                                                                                                  aCATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATg            <  2:2063285/148‑1 (MQ=255)
                                                                                                                                     tAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACa         <  1:733440/148‑1 (MQ=255)
                                                                                                                                      aaaGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAt        >  1:368741/1‑148 (MQ=255)
                                                                                                                                        aGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTa      <  2:483565/148‑1 (MQ=255)
                                                                                                                                        aGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTa      >  2:1986225/1‑148 (MQ=255)
                                                                                                                                         gCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTaa     >  2:3075838/1‑148 (MQ=255)
                                                                                                                                         gCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTaa     >  2:3075800/1‑148 (MQ=255)
                                                                                                                                         gCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTaa     <  2:3381791/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCA  >  CP009273/2736571‑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: