Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,495,516 C→A 17.4% S297* (TCG→TAG)  frlB → fructoselysine‑6‑P‑deglycase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,495,5160CA17.4% 810.2 / 99.6 259S297* (TCG→TAG) frlBfructoselysine‑6‑P‑deglycase
Reads supporting (aligned to +/- strand):  ref base C (121/92);  new base A (21/25);  total (142/117)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

CCGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTCCCGGCCCG  >  CP009273/3495367‑3495665
                                                                                                                                                     |                                                                                                                                                     
   ccATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCa                                                                                                                                                    <  1:4415647/150‑1 (MQ=255)
    cATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCaa                                                                                                                                                   <  1:1838569/150‑1 (MQ=255)
     aTGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAATTCGCCACACCACCGAACAAGCCATAAACTTTGTTAAACAACCTACTGACAACGTGATCGTAATCGATTAAGACGAAATTTAGCAAg                                                                                                                                                  >  1:401656/1‑150 (MQ=255)
      tGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAgg                                                                                                                                                 <  1:3722854/150‑1 (MQ=255)
      tGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAgg                                                                                                                                                 <  2:2718153/150‑1 (MQ=255)
         cccGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCt                                                                                                                                              <  1:2894229/150‑1 (MQ=255)
                     tGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGtggct                                                                                                                                  <  2:4363014/150‑1 (MQ=255)
                     tGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGtggct                                                                                                                                  <  2:2031189/150‑1 (MQ=255)
                      gTCGAACCGGGCGTTCCGTTCCTGTTCCAGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGtggctg                                                                                                                                 >  2:2440576/1‑150 (MQ=255)
                        cGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGtggctggc                                                                                                                               >  2:3895072/1‑150 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              >  2:2873333/1‑150 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              >  1:3113766/1‑150 (MQ=255)
                               ggCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTc                                                                                                                        >  2:3617273/1‑150 (MQ=255)
                                 cgttccgttcCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCt                                                                                                                      >  2:2074326/1‑150 (MQ=255)
                                 cgttccgttcCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCt                                                                                                                      >  1:3581950/1‑150 (MQ=255)
                                 cgttccgttcCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCt                                                                                                                      >  1:3955869/1‑150 (MQ=255)
                                 cgttccgttcCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCt                                                                                                                      >  2:1069020/1‑150 (MQ=255)
                                   ttccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGa                                                                                                                    <  1:3308973/150‑1 (MQ=255)
                                     ccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATg                                                                                                                  <  1:1077383/150‑1 (MQ=255)
                                     ccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATg                                                                                                                  <  2:3369737/150‑1 (MQ=255)
                                     ccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATg                                                                                                                  >  2:1116201/1‑150 (MQ=255)
                                      cgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGt                                                                                                                 >  2:1628072/1‑150 (MQ=255)
                                       gttcctgttcctGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGtt                                                                                                                >  1:3809072/1‑150 (MQ=255)
                                        ttcctgttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTc                                                                                                               <  2:3754558/150‑1 (MQ=255)
                                        ttcctgttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTc                                                                                                               >  1:4470784/1‑150 (MQ=255)
                                           ctgttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTg                                                                                                            <  2:2937936/150‑1 (MQ=255)
                                             gttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGcc                                                                                                          >  2:2366632/1‑150 (MQ=255)
                                              ttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCa                                                                                                         <  2:1196210/150‑1 (MQ=255)
                                               tcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCaa                                                                                                        <  1:1108013/150‑1 (MQ=255)
                                                  tgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATgg                                                                                                     <  2:2888045/150‑1 (MQ=255)
                                                     tCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGt                                                                                                  >  2:521840/1‑150 (MQ=255)
                                                       ggCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTgg                                                                                                >  1:1340111/1‑150 (MQ=255)
                                                        gCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGc                                                                                               >  1:245485/1‑150 (MQ=255)
                                                        gCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGc                                                                                               >  2:3240442/1‑150 (MQ=255)
                                                         cAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGct                                                                                              <  1:1197429/150‑1 (MQ=255)
                                                           atgatgAAAGTCGCCTCACCACCGACCGCGCCATTAACTTTTTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGATGGCACCTTTCCTGTTTTTCGTTCCTATGGAGTGGatct                                                                                            <  2:1852983/150‑1 (MQ=255)
                                                           atgatgAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGctct                                                                                            <  1:915612/150‑1 (MQ=255)
                                                              atgaAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGct                                                                                         >  1:4251951/1‑150 (MQ=255)
                                                                 aaaGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGctact                                                                                      >  1:4374207/1‑150 (MQ=255)
                                                                     tCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCt                                                                                  >  1:92346/1‑150 (MQ=255)
                                                                       gCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGt                                                                                >  2:3952247/1‑150 (MQ=255)
                                                                       gCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGt                                                                                >  2:1026546/1‑150 (MQ=255)
                                                                        ccACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTc                                                                               <  2:182395/150‑1 (MQ=255)
                                                                         cacaCCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCt                                                                              <  1:2777207/150‑1 (MQ=255)
                                                                          acacCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTa                                                                             <  2:1315959/150‑1 (MQ=255)
                                                                            accaccGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTAtt                                                                           <  2:3960389/150‑1 (MQ=255)
                                                                             ccaccGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTAttt                                                                          <  2:235896/150‑1 (MQ=255)
                                                                                ccGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACa                                                                       >  1:1558735/1‑150 (MQ=255)
                                                                                   aaCGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAg                                                                    >  1:3866137/1‑150 (MQ=255)
                                                                                   aaCGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAg                                                                    >  1:581098/1‑150 (MQ=255)
                                                                                    aCGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGa                                                                   >  1:1524247/1‑150 (MQ=255)
                                                                                     cgcgCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGAt                                                                  <  1:2184713/150‑1 (MQ=255)
                                                                                       cgcCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATca                                                                >  2:2975240/1‑150 (MQ=255)
                                                                                        gcCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATcac                                                               >  1:154229/1‑150 (MQ=255)
                                                                                         ccATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATATTTTACGCCGAAATTTAGCTAGGGCTGATCCCGTGTCTGGCACCGTTCCCGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATcaca                                                              <  2:1511783/150‑1 (MQ=255)
                                                                                           aTTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAc                                                            >  1:458825/1‑150 (MQ=255)
                                                                                             tAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAccc                                                          <  2:186927/150‑1 (MQ=255)
                                                                                               aCTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCgg                                                        <  2:663055/150‑1 (MQ=255)
                                                                                                cTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGa                                                       <  1:4120825/150‑1 (MQ=255)
                                                                                                 tttGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGAt                                                      <  2:1903/150‑1 (MQ=255)
                                                                                                 tttGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGAt                                                      <  1:2984313/150‑1 (MQ=255)
                                                                                                      tAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcg                                                 >  2:4114189/1‑150 (MQ=255)
                                                                                                       aaaCAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgc                                                >  2:3396861/1‑150 (MQ=255)
                                                                                                       aaaCAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgc                                                >  2:1030856/1‑150 (MQ=255)
                                                                                                        aaCAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgcc                                               >  2:593475/1‑150 (MQ=255)
                                                                                                        aaCAGCGTACTGACAACGCGATCGTCCTCGATTACGCCGAAATCTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGCTGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTATACTTACAAAGATCACAACCCGGATGAAcgcc                                               <  2:4413521/150‑1 (MQ=255)
                                                                                                          cAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgccgc                                             <  1:352610/150‑1 (MQ=255)
                                                                                                             cGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTAt                                          >  2:2626683/1‑150 (MQ=255)
                                                                                                                aCTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTAc                                       >  2:3937661/1‑150 (MQ=255)
                                                                                                                aCTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTAc                                       <  2:1493061/150‑1 (MQ=255)
                                                                                                                 cTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACg                                      >  1:3340140/1‑150 (MQ=255)
                                                                                                                    aCAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGGCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACggtg                                   <  1:2612696/150‑1 (MQ=255)
                                                                                                                    aCAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACggtg                                   <  2:3696652/150‑1 (MQ=255)
                                                                                                                       aCGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTc                                <  2:4274190/150‑1 (MQ=255)
                                                                                                                        cGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCt                               <  2:847530/150‑1 (MQ=255)
                                                                                                                        cGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCt                               >  2:252232/1‑150 (MQ=255)
                                                                                                                         gTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtg                              >  2:4112788/1‑150 (MQ=255)
                                                                                                                          tGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtgg                             >  2:3828067/1‑150 (MQ=255)
                                                                                                                            aTCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtggtg                           <  2:771205/150‑1 (MQ=255)
                                                                                                                            aTCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtggtg                           <  1:3416978/150‑1 (MQ=255)
                                                                                                                             tCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtggtgg                          >  2:825438/1‑150 (MQ=255)
                                                                                                                               gTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGaa                        >  1:3681992/1‑150 (MQ=255)
                                                                                                                               gTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGaa                        >  1:3488614/1‑150 (MQ=255)
                                                                                                                                tCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAat                       >  1:518035/1‑150 (MQ=255)
                                                                                                                                 cATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAata                      >  2:2494109/1‑150 (MQ=255)
                                                                                                                                  aTCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAatat                     >  1:1357764/1‑150 (MQ=255)
                                                                                                                                   tCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATAtt                    <  1:2558503/150‑1 (MQ=255)
                                                                                                                                     gATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTaa                  >  2:4542959/1‑150 (MQ=255)
                                                                                                                                       ttACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATc                <  1:2012144/150‑1 (MQ=255)
                                                                                                                                        tACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATcc               >  2:2193265/1‑150 (MQ=255)
                                                                                                                                        tACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCCGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATcc               >  1:4406549/1‑150 (MQ=255)
                                                                                                                                         aCGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATccc              <  1:2180579/150‑1 (MQ=255)
                                                                                                                                          cGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCct             <  2:3113766/150‑1 (MQ=255)
                                                                                                                                           gCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCctc            >  1:2335290/1‑150 (MQ=255)
                                                                                                                                               aaaTTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTccc        >  2:1937420/1‑150 (MQ=255)
                                                                                                                                               aaaTTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTccc        <  2:2370544/150‑1 (MQ=255)
                                                                                                                                               aaaTTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTccc        <  1:2641000/150‑1 (MQ=255)
                                                                                                                                                 aTTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccgg      >  2:2373760/1‑150 (MQ=255)
                                                                                                                                                  tttAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccggc     >  1:4324189/1‑150 (MQ=255)
                                                                                                                                                     cGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccggcccg  >  2:2754549/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTCCCGGCCCG  >  CP009273/3495367‑3495665

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: