Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,596,499 G→A 6.0% G441G (GGG→GGA hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,596,4990GA 6.0% 704.1 / 22.7 215G441G (GGG→GGAhyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (106/96);  new base A (8/5);  total (114/101)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 225 total aligned reads displayed.

TGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAG  >  CP009273/2596350‑2596648
                                                                                                                                                     |                                                                                                                                                     
  ccTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCt                                                                                                                                                     <  2:2361657/150‑1 (MQ=255)
   cTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTg                                                                                                                                                    <  1:2888902/150‑1 (MQ=255)
     cgcgATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGc                                                                                                                                                  >  2:3172915/1‑150 (MQ=255)
          tatCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTaa                                                                                                                                             <  1:999215/150‑1 (MQ=255)
           atCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAAt                                                                                                                                            >  2:5178661/1‑150 (MQ=255)
           atCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAAt                                                                                                                                            >  2:4794415/1‑150 (MQ=255)
           atCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACGAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAAt                                                                                                                                            >  2:4794540/1‑150 (MQ=255)
            tCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAAtg                                                                                                                                           <  1:2032579/150‑1 (MQ=255)
             cAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAAtgt                                                                                                                                          >  2:2858039/1‑150 (MQ=255)
                ccATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCt                                                                                                                                       <  1:42178/150‑1 (MQ=255)
                 cATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAATGTGCtt                                                                                                                                      <  2:352084/150‑1 (MQ=255)
                  aTTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTc                                                                                                                                     >  1:338759/1‑150 (MQ=255)
                            tGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGt                                                                                                                           <  2:2951106/150‑1 (MQ=255)
                             gAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTa                                                                                                                          >  2:2478965/1‑150 (MQ=255)
                              aaTGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTAc                                                                                                                         >  2:4675058/1‑150 (MQ=255)
                               aTGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACg                                                                                                                        >  1:3015804/1‑150 (MQ=255)
                                 ggTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGt                                                                                                                      >  2:927459/1‑150 (MQ=255)
                                   ttttATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTAt                                                                                                                    <  2:1713903/150‑1 (MQ=255)
                                    tttATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTAtt                                                                                                                   <  1:2158121/150‑1 (MQ=255)
                                    tttATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTAtt                                                                                                                   >  2:4229152/1‑150 (MQ=255)
                                    tttATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTAAGGTAtt                                                                                                                   >  2:887373/1‑150 (MQ=255)
                                     ttATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTa                                                                                                                  >  2:3978259/1‑150 (MQ=255)
                                      tATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTAc                                                                                                                 <  1:1338615/150‑1 (MQ=255)
                                         cAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTCATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACttt                                                                                                              >  2:1204579/1‑150 (MQ=255)
                                         cAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACttt                                                                                                              >  2:1056738/1‑150 (MQ=255)
                                         cAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACttt                                                                                                              <  1:3253083/150‑1 (MQ=255)
                                         cAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACttt                                                                                                              <  1:4000132/150‑1 (MQ=255)
                                           gCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCGTCGTAAAAATGTACGGTATTACTTTCt                                                                                                            >  2:3131321/1‑150 (MQ=255)
                                            cGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCtg                                                                                                           >  1:3548994/1‑150 (MQ=255)
                                            cGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCtg                                                                                                           <  2:3230729/150‑1 (MQ=255)
                                                tGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTgtggt                                                                                                       <  1:3657770/150‑1 (MQ=255)
                                                          ggCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTAATTTCTGTGGTGCGCCGCGCa                                                                                             >  2:4892714/1‑150 (MQ=255)
                                                           gCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAg                                                                                            <  2:1331624/150‑1 (MQ=255)
                                                            cAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGt                                                                                           <  1:537977/150‑1 (MQ=255)
                                                                cGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTacac                                                                                       >  2:669718/1‑150 (MQ=255)
                                                                 gCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTCTGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTacaca                                                                                      >  2:59725/1‑150 (MQ=255)
                                                                    gTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTAGTGTCATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGc                                                                                   >  1:459677/1‑150 (MQ=255)
                                                                     ttCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCt                                                                                  <  1:4921404/150‑1 (MQ=255)
                                                                     ttCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCt                                                                                  <  1:1170331/150‑1 (MQ=255)
                                                                     ttCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCt                                                                                  <  1:291065/150‑1 (MQ=255)
                                                                        tcACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGaa                                                                               >  1:1985612/1‑150 (MQ=255)
                                                                           cTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAagag                                                                            >  1:258944/1‑150 (MQ=255)
                                                                           cTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAagag                                                                            <  1:3006042/150‑1 (MQ=255)
                                                                            tAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAgg                                                                           >  2:2829624/1‑150 (MQ=255)
                                                                             aaGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGc                                                                          >  2:56874/1‑150 (MQ=255)
                                                                                tCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGcaca                                                                       <  1:1516730/150‑1 (MQ=255)
                                                                                tCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGcaca                                                                       <  1:671914/150‑1 (MQ=255)
                                                                                tCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGcaca                                                                       <  1:3335933/150‑1 (MQ=255)
                                                                                  gtgtGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAgg                                                                     <  2:230751/150‑1 (MQ=255)
                                                                                   tgtgGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAATAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGa                                                                    <  1:456490/150‑1 (MQ=255)
                                                                                    gtgGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGaa                                                                   >  2:377371/1‑150 (MQ=255)
                                                                                     tgGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAg                                                                  <  2:2340669/150‑1 (MQ=255)
                                                                                      ggAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGt                                                                 >  2:2891011/1‑150 (MQ=255)
                                                                                      ggAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGt                                                                 <  2:3722043/150‑1 (MQ=255)
                                                                                       gAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTg                                                                <  1:4458170/150‑1 (MQ=255)
                                                                                        aaGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGc                                                               <  2:2992997/150‑1 (MQ=255)
                                                                                          gCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACAGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCa                                                             >  2:3268185/1‑150 (MQ=255)
                                                                                             gTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAAt                                                          >  2:5133478/1‑150 (MQ=255)
                                                                                               aGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATAc                                                        <  2:64496/150‑1 (MQ=255)
                                                                                                gcgcTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACg                                                       >  1:587019/1‑150 (MQ=255)
                                                                                                  gcTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgat                                                     <  1:2396394/150‑1 (MQ=255)
                                                                                                  gcTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgat                                                     >  1:2135218/1‑150 (MQ=255)
                                                                                                   cTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgatg                                                    >  1:2717788/1‑150 (MQ=255)
                                                                                                   cTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgatg                                                    <  1:3902247/150‑1 (MQ=255)
                                                                                                   cTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgatg                                                    <  2:1300312/150‑1 (MQ=255)
                                                                                                     aCAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACgatgat                                                  >  1:2987441/1‑150 (MQ=255)
                                                                                                        aCTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGAtcgt                                               >  1:4504067/1‑150 (MQ=255)
                                                                                                         cTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGAtcgtc                                              >  1:1148822/1‑150 (MQ=255)
                                                                                                         cTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGAtcgtc                                              <  2:77800/150‑1 (MQ=255)
                                                                                                         cTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGAtcgtc                                              >  2:2187601/1‑150 (MQ=255)
                                                                                                          ttGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGAtcgtcg                                             >  1:2416003/1‑150 (MQ=255)
                                                                                                               ggTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTCCACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATg                                        >  2:177710/1‑150 (MQ=255)
                                                                                                                gTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGc                                       <  2:4074984/150‑1 (MQ=255)
                                                                                                                gTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGACTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGc                                       <  2:2183912/150‑1 (MQ=255)
                                                                                                                 tCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCt                                      <  1:2683099/150‑1 (MQ=255)
                                                                                                                 tCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCt                                      <  1:4422875/150‑1 (MQ=255)
                                                                                                                 tCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCt                                      <  1:4575461/150‑1 (MQ=255)
                                                                                                                 tCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCt                                      <  1:4416754/150‑1 (MQ=255)
                                                                                                                    tATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACt                                   <  2:3926593/150‑1 (MQ=255)
                                                                                                                        gCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTc                               <  2:3318935/150‑1 (MQ=255)
                                                                                                                           aTGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTcgcg                            >  1:1062/1‑150 (MQ=255)
                                                                                                                            tGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCgg                           >  1:1002767/1‑150 (MQ=255)
                                                                                                                             ggTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGACGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGc                          >  2:1229268/1‑150 (MQ=255)
                                                                                                                               tAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCAc                        <  2:95394/150‑1 (MQ=255)
                                                                                                                               tAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCAc                        >  1:1608838/1‑150 (MQ=255)
                                                                                                                                 aTGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCActc                      >  2:1793838/1‑150 (MQ=255)
                                                                                                                                         aGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTAtt              >  1:1685666/1‑150 (MQ=255)
                                                                                                                                         aGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTAtt              >  2:1592819/1‑150 (MQ=255)
                                                                                                                                         aGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTAtt              <  1:3090786/150‑1 (MQ=255)
                                                                                                                                           tCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTaa            >  2:5140378/1‑150 (MQ=255)
                                                                                                                                           tCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTaa            <  2:3666575/150‑1 (MQ=255)
                                                                                                                                           tCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTaa            >  2:2562442/1‑150 (MQ=255)
                                                                                                                                             aCTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAAtt          <  1:3986017/150‑1 (MQ=255)
                                                                                                                                               tggtggGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTgc        >  2:1158152/1‑150 (MQ=255)
                                                                                                                                               tggtggGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTgc        <  2:3548994/150‑1 (MQ=255)
                                                                                                                                                   gggCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCtt    <  2:201143/150‑1 (MQ=255)
                                                                                                                                                    ggCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTCTTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTa   >  2:1076593/1‑150 (MQ=255)
                                                                                                                                                    ggCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTa   <  2:5133501/150‑1 (MQ=255)
                                                                                                                                                    ggCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTa   <  1:257933/150‑1 (MQ=255)
                                                                                                                                                     gCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAg  <  2:2987308/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGGCAGTCGCTGTTCTCACTAAGTCGTGTGGAAGCCGTAGCGCTACAACTTGCGGGTCCTATTGCTATGGTAATGCTGGCAGTCACTGGTGGGCTGGCAGTAATGTGCTTCGTAAAAATGTACGGTATTACTTTCTGTGGTGCGCCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAG  >  CP009273/2596350‑2596648

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: