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,7140CA98.1% 827.1 / 2.2 255R325S (CGC→AGC) yfiNputative membrane‑anchored diguanylate cyclase with an N‑terminal periplasmic domain
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (126/124);  total (129/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 255 total aligned reads displayed.

CCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCAT  >  CP009273/2736569‑2736859
                                                                                                                                                 |                                                                                                                                                 
ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCt                                                                                                                                               <  2:923664/150‑1 (MQ=255)
ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                               <  2:4144459/150‑1 (MQ=255)
ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                               <  2:3682972/150‑1 (MQ=255)
ccGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCt                                                                                                                                               <  2:4854679/150‑1 (MQ=255)
 cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTg                                                                                                                                              >  1:1100176/1‑150 (MQ=255)
 cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTg                                                                                                                                              <  1:4742124/150‑1 (MQ=255)
 cGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTg                                                                                                                                              <  2:148900/150‑1 (MQ=255)
  gAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTgg                                                                                                                                             >  2:1325720/1‑150 (MQ=255)
     aaaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGggcg                                                                                                                                          <  2:2718404/150‑1 (MQ=255)
     aaaCGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGggcg                                                                                                                                          <  2:2410046/150‑1 (MQ=255)
         gTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGa                                                                                                                                      <  1:2947415/150‑1 (MQ=255)
         gTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGa                                                                                                                                      <  1:513908/150‑1 (MQ=255)
          tCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGAt                                                                                                                                     >  1:4089641/1‑150 (MQ=255)
            ggCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGa                                                                                                                                   <  1:1006899/150‑1 (MQ=255)
                  aCTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGc                                                                                                                             >  1:1576751/1‑150 (MQ=255)
                   cTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCt                                                                                                                            >  1:3679147/1‑150 (MQ=255)
                     aTTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTAt                                                                                                                          <  2:2581488/150‑1 (MQ=255)
                       ttCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATgg                                                                                                                        <  1:2606977/150‑1 (MQ=255)
                        tCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGt                                                                                                                       <  2:4049123/150‑1 (MQ=255)
                             aTGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGctct                                                                                                                  >  2:3724706/1‑150 (MQ=255)
                               ggCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTat                                                                                                                >  2:2616398/1‑150 (MQ=255)
                                 cGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTatga                                                                                                              <  1:553206/150‑1 (MQ=255)
                                        ttCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAg                                                                                                       >  1:2231295/1‑150 (MQ=255)
                                          cAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTc                                                                                                     <  1:1475952/150‑1 (MQ=255)
                                           aaaTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCg                                                                                                    >  2:3253598/1‑150 (MQ=255)
                                             aTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                  >  1:4870401/1‑150 (MQ=255)
                                             aTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                  <  1:1630649/150‑1 (MQ=255)
                                             aTACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGa                                                                                                  <  1:936810/150‑1 (MQ=255)
                                              tACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGaa                                                                                                 <  2:3598479/150‑1 (MQ=255)
                                               aCATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAAt                                                                                                >  1:5038889/1‑150 (MQ=255)
                                                cATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATc                                                                                               <  2:244024/150‑1 (MQ=255)
                                                 aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCt                                                                                              <  1:2444088/150‑1 (MQ=255)
                                                 aTCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCt                                                                                              <  2:1636753/150‑1 (MQ=255)
                                                   cAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                            <  1:3477121/150‑1 (MQ=255)
                                                   cAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGa                                                                                            <  1:2722673/150‑1 (MQ=255)
                                                    aaTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                           >  2:3196595/1‑150 (MQ=255)
                                                    aaTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                           <  2:5014793/150‑1 (MQ=255)
                                                    aaTGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATATGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGACTCTGaa                                                                                           >  1:4060481/1‑150 (MQ=255)
                                                    aaTGAGACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGaa                                                                                           <  2:1959676/150‑1 (MQ=255)
                                                        aTACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                       <  1:4619211/150‑1 (MQ=255)
                                                        aTACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                       >  2:1598614/1‑150 (MQ=255)
                                                        aTACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                       >  2:4604346/1‑150 (MQ=255)
                                                        aTACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgc                                                                                       <  1:2809466/150‑1 (MQ=255)
                                                            cTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTgcagca                                                                                   >  1:970188/1‑150 (MQ=255)
                                                             tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                  >  2:4616472/1‑150 (MQ=255)
                                                             tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                  >  1:4468519/1‑150 (MQ=255)
                                                             tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                  >  2:1238996/1‑150 (MQ=255)
                                                             tGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGcagcag                                                                                  >  1:1563262/1‑150 (MQ=255)
                                                                      gCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCa                                                                         <  1:5059839/150‑1 (MQ=255)
                                                                       cGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAg                                                                        >  2:2767066/1‑150 (MQ=255)
                                                                       cGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAg                                                                        >  2:745318/1‑150 (MQ=255)
                                                                           gggCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACt                                                                    >  1:2361139/1‑150 (MQ=255)
                                                                           gggCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACt                                                                    >  2:4369648/1‑150 (MQ=255)
                                                                             gCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGa                                                                  >  2:3976180/1‑150 (MQ=255)
                                                                              cGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCCGCAGATATGTTCAGCACTGac                                                                 >  1:689305/1‑150 (MQ=255)
                                                                                aTAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGacac                                                               >  1:1097663/1‑150 (MQ=255)
                                                                                 tAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGAcaca                                                              <  1:3905878/150‑1 (MQ=255)
                                                                                 tAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGAcaca                                                              >  2:4291047/1‑150 (MQ=255)
                                                                                  agagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCCGCAGATATGTTCAGCACTGACACaa                                                             >  2:2060808/1‑150 (MQ=255)
                                                                                  agagTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACaa                                                             >  2:3900656/1‑150 (MQ=255)
                                                                                    agTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAAt                                                           <  2:1824647/150‑1 (MQ=255)
                                                                                       cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                        >  1:4226525/1‑150 (MQ=255)
                                                                                       cTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCtt                                                        >  2:693632/1‑150 (MQ=255)
                                                                                        ttgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCttt                                                       <  2:130811/150‑1 (MQ=255)
                                                                                         tgattgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTa                                                      >  1:3433317/1‑150 (MQ=255)
                                                                                           attgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                    >  2:3916663/1‑150 (MQ=255)
                                                                                           attgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                    >  2:4457233/1‑150 (MQ=255)
                                                                                           attgaAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGGTCAGCACTGACACAAATCTTTAAt                                                    <  2:1305943/150‑1 (MQ=255)
                                                                                           attgaAATAGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAAt                                                    >  2:225222/1‑150 (MQ=255)
                                                                                                aaTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTccc                                               >  1:4895720/1‑150 (MQ=255)
                                                                                                 aTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCg                                              >  1:1015530/1‑150 (MQ=255)
                                                                                                 aTCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCg                                              >  1:1155348/1‑150 (MQ=255)
                                                                                                   cGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGtt                                            <  1:3550732/150‑1 (MQ=255)
                                                                                                    gCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGttt                                           <  2:4667309/150‑1 (MQ=255)
                                                                                                     cAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                          <  2:2708992/150‑1 (MQ=255)
                                                                                                     cAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                          <  2:4519977/150‑1 (MQ=255)
                                                                                                     cAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                          >  1:4922452/1‑150 (MQ=255)
                                                                                                     cAAAACGGTGAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTg                                          <  2:2832299/150‑1 (MQ=255)
                                                                                                      aaaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGa                                         <  2:4768606/150‑1 (MQ=255)
                                                                                                      aaaaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGa                                         >  2:4083569/1‑150 (MQ=255)
                                                                                                        aaCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATc                                       >  2:3699692/1‑150 (MQ=255)
                                                                                                         aCGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCt                                      <  2:3063430/150‑1 (MQ=255)
                                                                                                          cGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCtt                                     <  2:1949309/150‑1 (MQ=255)
                                                                                                           ggTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                    >  2:1312513/1‑150 (MQ=255)
                                                                                                           ggTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTc                                    <  2:1737605/150‑1 (MQ=255)
                                                                                                                gCTGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAAt                               <  2:3800259/150‑1 (MQ=255)
                                                                                                                  tGAATTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATgg                             >  2:4200530/1‑150 (MQ=255)
                                                                                                                    aaTTTGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGcc                           >  1:4037998/1‑150 (MQ=255)
                                                                                                                        tGGCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCa                       <  2:5143641/150‑1 (MQ=255)
                                                                                                                         ggCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAg                      <  2:4452954/150‑1 (MQ=255)
                                                                                                                         ggCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAg                      >  1:1185661/1‑150 (MQ=255)
                                                                                                                          gCGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGa                     <  2:789818/150‑1 (MQ=255)
                                                                                                                           cGGGCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGac                    >  1:3268164/1‑150 (MQ=255)
                                                                                                                            gggCTGCGACATAAAGCATACAGCCTGGGCGGCGATNAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGacc                   >  1:2680649/1‑150 (MQ=255)
                                                                                                                            gggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGacc                   >  1:4274933/1‑150 (MQ=255)
                                                                                                                             ggCTGCGACATAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGacca                  >  2:3162687/1‑150 (MQ=255)
                                                                                                                                    aCATAAAGCAGACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTTTGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCTTAATGGCCATCAGACCACCATGAc           <  1:4795727/150‑1 (MQ=255)
                                                                                                                                       tAAAGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACAtt        >  1:2959834/1‑150 (MQ=255)
                                                                                                                                        aaaGCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTa       >  2:2389779/1‑150 (MQ=255)
                                                                                                                                           gCATACAGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGc    <  2:1576751/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGAAAAACGTCGGCGTTACTATTTCTTGATGGCGATAATTTCAAATACATCAATGATACCTGGGGTCATGCGACGGGCGATAGAGTCTTGATTGAAATCGCAAAACGGTTAGCTGAATTTGGCGGGCTGCGACATAAAGCATACCGCCTGGGCGGCGATGAATTCGCTATGGTGCTCTATGATGTACAGTCGGAATCTGAAGTGCAGCAGATATGTTCAGCACTGACACAAATCTTTAATCTCCCGTTTGATCTTCATAATGGCCATCAGACCACCATGACATTAAGCAT  >  CP009273/2736569‑2736859

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: