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,7140CA97.6% 700.9 / 4.3 220R325S (CGC→AGC) yfiNputative membrane‑anchored diguanylate cyclase with an N‑terminal periplasmic domain
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base A (110/104);  total (115/105)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 221 total aligned reads displayed.

CCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCATTG  >  CP009273/2736569‑2736861
                                                                                                                                                 |                                                                                                                                                   
ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGc                                                                                                                                                   <  2:3382799/148‑1 (MQ=255)
 cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGcc                                                                                                                                                  >  2:2478380/1‑148 (MQ=255)
  gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                                 >  2:1107242/1‑148 (MQ=255)
  gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                                 >  2:3022621/1‑148 (MQ=255)
  gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCGTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                                 >  2:2857960/1‑148 (MQ=255)
   aaaaaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTg                                                                                                                                                <  2:2854138/148‑1 (MQ=255)
     aaaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTggg                                                                                                                                              <  1:2485773/148‑1 (MQ=255)
        cGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGggcgg                                                                                                                                           <  2:2440675/148‑1 (MQ=255)
            ggCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGAt                                                                                                                                       <  2:151247/148‑1 (MQ=255)
               gTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGaa                                                                                                                                    <  2:171460/148‑1 (MQ=255)
                ttACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAAt                                                                                                                                   >  2:85194/1‑148 (MQ=255)
                ttACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAt                                                                                                                                   <  1:1061366/148‑1 (MQ=255)
                ttACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAt                                                                                                                                   <  2:1016499/148‑1 (MQ=255)
                ttACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAt                                                                                                                                   <  2:415082/148‑1 (MQ=255)
                 tACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGACtt                                                                                                                                  >  2:2037978/1‑148 (MQ=255)
                 tACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAtt                                                                                                                                  >  1:3407514/1‑148 (MQ=255)
                 tACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAtt                                                                                                                                  >  1:1175709/1‑148 (MQ=255)
                 tACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAAtt                                                                                                                                  >  1:3282730/1‑148 (MQ=255)
                  aCTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTc                                                                                                                                 >  1:668442/1‑148 (MQ=255)
                  aCTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTc                                                                                                                                 >  1:1181036/1‑148 (MQ=255)
                   cTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCg                                                                                                                                <  1:449160/148‑1 (MQ=255)
                     aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                              >  1:1857102/1‑148 (MQ=255)
                     aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                              >  1:1856457/1‑148 (MQ=255)
                      tttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTa                                                                                                                             <  2:582873/148‑1 (MQ=255)
                      tttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTa                                                                                                                             >  1:2827624/1‑148 (MQ=255)
                       ttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTAt                                                                                                                            >  1:393394/1‑148 (MQ=255)
                       ttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTAt                                                                                                                            >  2:1628550/1‑148 (MQ=255)
                        tCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATg                                                                                                                           <  2:1549365/148‑1 (MQ=255)
                         cTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATgg                                                                                                                          >  2:3135652/1‑148 (MQ=255)
                         cTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATgg                                                                                                                          >  2:495546/1‑148 (MQ=255)
                           tGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTg                                                                                                                        <  2:1265296/148‑1 (MQ=255)
                             aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGct                                                                                                                      >  1:3368055/1‑148 (MQ=255)
                              tGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGctc                                                                                                                     >  1:2564737/1‑148 (MQ=255)
                               ggCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGctct                                                                                                                    >  1:596083/1‑148 (MQ=255)
                               ggCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGctct                                                                                                                    >  1:596097/1‑148 (MQ=255)
                                     aaTTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTatgatg                                                                                                              >  2:2494189/1‑148 (MQ=255)
                                      aTTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGt                                                                                                             >  2:3233256/1‑148 (MQ=255)
                                       tttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTa                                                                                                            <  2:144650/148‑1 (MQ=255)
                                       tttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTa                                                                                                            >  2:3093561/1‑148 (MQ=255)
                                           aaaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGt                                                                                                        >  1:1092000/1‑148 (MQ=255)
                                           aaaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGt                                                                                                        <  2:1637964/148‑1 (MQ=255)
                                             aTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCg                                                                                                      <  2:750633/148‑1 (MQ=255)
                                               aCATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                    >  1:3280596/1‑148 (MQ=255)
                                               aCATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                    <  2:236485/148‑1 (MQ=255)
                                                 aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAAt                                                                                                  <  1:658528/148‑1 (MQ=255)
                                                 aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAAt                                                                                                  >  1:2892195/1‑148 (MQ=255)
                                                    aaTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTg                                                                                               >  1:158247/1‑148 (MQ=255)
                                                     aTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                              >  2:2312861/1‑148 (MQ=255)
                                                      tGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                             <  1:2847619/148‑1 (MQ=255)
                                                       gATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAg                                                                                            <  1:2569673/148‑1 (MQ=255)
                                                         tACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTCCAGTCGGAATCTGAAGTg                                                                                          >  2:1255587/1‑148 (MQ=255)
                                                         tACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTg                                                                                          >  2:1255695/1‑148 (MQ=255)
                                                            cTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgcag                                                                                       <  1:3154042/148‑1 (MQ=255)
                                                              ggggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgcagca                                                                                     >  2:2180787/1‑148 (MQ=255)
                                                               gggTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                    >  2:2576499/1‑148 (MQ=255)
                                                                 gTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGat                                                                                  <  2:2945886/148‑1 (MQ=255)
                                                                 gTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGat                                                                                  <  1:3477587/148‑1 (MQ=255)
                                                                     tGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGt                                                                              >  2:2918292/1‑148 (MQ=255)
                                                                      gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGtt                                                                             <  1:1240794/148‑1 (MQ=255)
                                                                           gggCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCa                                                                        <  1:3364124/148‑1 (MQ=255)
                                                                            ggCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCAc                                                                       <  2:307631/148‑1 (MQ=255)
                                                                              cGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTg                                                                     <  2:808590/148‑1 (MQ=255)
                                                                                 tAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGaca                                                                  <  1:3231465/148‑1 (MQ=255)
                                                                                   gagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGAcaca                                                                >  2:1772030/1‑148 (MQ=255)
                                                                                   gagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGAcaca                                                                >  1:2157887/1‑148 (MQ=255)
                                                                                       cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATc                                                            <  2:106216/148‑1 (MQ=255)
                                                                                         tgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                          >  1:2312220/1‑148 (MQ=255)
                                                                                         tgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                          <  1:3135652/148‑1 (MQ=255)
                                                                                            ttgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTaa                                                       <  1:1537897/148‑1 (MQ=255)
                                                                                             tgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                      <  1:822045/148‑1 (MQ=255)
                                                                                             tgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                      >  1:3040146/1‑148 (MQ=255)
                                                                                                 aTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTcc                                                  <  2:1675031/148‑1 (MQ=255)
                                                                                                 aTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTcc                                                  >  2:2653970/1‑148 (MQ=255)
                                                                                                   cGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCg                                                >  1:2772512/1‑148 (MQ=255)
                                                                                                      aaaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGttt                                             <  1:1918280/148‑1 (MQ=255)
                                                                                                      aaaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGttt                                             <  1:1918644/148‑1 (MQ=255)
                                                                                                         aCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGAt                                          >  2:200893/1‑148 (MQ=255)
                                                                                                         aCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGAt                                          >  1:1889268/1‑148 (MQ=255)
                                                                                                          cGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAATTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATc                                         >  2:3440168/1‑148 (MQ=255)
                                                                                                          cGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATc                                         <  1:2786460/148‑1 (MQ=255)
                                                                                                          cGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATc                                         <  1:3023249/148‑1 (MQ=255)
                                                                                                             ttAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                      >  2:634825/1‑148 (MQ=255)
                                                                                                              tAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCa                                     >  2:2256480/1‑148 (MQ=255)
                                                                                                               aGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCAt                                    <  2:451479/148‑1 (MQ=255)
                                                                                                                    aaTTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATgg                               >  2:1399670/1‑148 (MQ=255)
                                                                                                                       ttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCa                            >  2:1152919/1‑148 (MQ=255)
                                                                                                                       ttGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCa                            >  1:101164/1‑148 (MQ=255)
                                                                                                                        tGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCAt                           <  1:2857960/148‑1 (MQ=255)
                                                                                                                            gggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGa                       <  1:2520364/148‑1 (MQ=255)
                                                                                                                             ggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGac                      >  2:915676/1‑148 (MQ=255)
                                                                                                                             ggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGac                      >  2:316703/1‑148 (MQ=255)
                                                                                                                             ggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGac                      >  2:311200/1‑148 (MQ=255)
                                                                                                                                    aCATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATg               <  2:2845074/148‑1 (MQ=255)
                                                                                                                                    aCATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATg               <  1:3481013/148‑1 (MQ=255)
                                                                                                                                     cATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGa              >  2:3404564/1‑148 (MQ=255)
                                                                                                                                     cATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGa              <  2:1857102/148‑1 (MQ=255)
                                                                                                                                     cATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGa              <  2:1856457/148‑1 (MQ=255)
                                                                                                                                       tAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACa            >  1:2956984/1‑148 (MQ=255)
                                                                                                                                        aaaGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAt           >  1:2063583/1‑148 (MQ=255)
                                                                                                                                         aaGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAtt          <  2:3282730/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCATTG  >  CP009273/2736569‑2736861

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: