Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,923,089 C→A 15.2% S52* (TCG→TAG)  sdaB → L‑serine dehydratase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,923,0890CA15.2% 701.5 / 82.9 222S52* (TCG→TAG) sdaBL‑serine dehydratase 2
Reads supporting (aligned to +/- strand):  ref base C (115/73);  new base A (15/19);  total (130/92)
Fisher's exact test for biased strand distribution p-value = 8.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

AGCGTATTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTC  >  CP009273/2922941‑2923238
                                                                                                                                                    |                                                                                                                                                     
  cGTATTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGgctcgct                                                                                                                                                    <  2:2962386/150‑1 (MQ=255)
    tATTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctct                                                                                                                                                  >  1:2275634/1‑150 (MQ=255)
    tATTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctct                                                                                                                                                  >  1:2021387/1‑150 (MQ=255)
     aTTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctctc                                                                                                                                                 >  2:3101294/1‑150 (MQ=255)
      ttCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctctct                                                                                                                                                >  1:670978/1‑150 (MQ=255)
      ttCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctctct                                                                                                                                                >  1:940137/1‑150 (MQ=255)
      ttCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGctctct                                                                                                                                                <  2:2101935/150‑1 (MQ=255)
          atatTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCAGACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGACCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGa                                                                                                                                            >  1:1536666/1‑150 (MQ=255)
           tatTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGAc                                                                                                                                           >  1:2159856/1‑150 (MQ=255)
            atTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGAcc                                                                                                                                          >  2:2414154/1‑150 (MQ=255)
             ttttCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCg                                                                                                                                         <  2:741832/150‑1 (MQ=255)
             ttttCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCg                                                                                                                                         <  2:2638367/150‑1 (MQ=255)
             ttttCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCg                                                                                                                                         >  1:2468422/1‑150 (MQ=255)
              tttCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCgg                                                                                                                                        >  2:1443081/1‑150 (MQ=255)
                tCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTa                                                                                                                                      <  1:444812/150‑1 (MQ=255)
                 cAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTaa                                                                                                                                     >  1:1665406/1‑150 (MQ=255)
                     aTCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGc                                                                                                                                 >  2:385659/1‑150 (MQ=255)
                     aTCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGc                                                                                                                                 >  1:1784323/1‑150 (MQ=255)
                     aTCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGc                                                                                                                                 <  2:1733721/150‑1 (MQ=255)
                     aTCGGCATTGGCCCTTCCAGGTCTCATACCGTTGGACCAAGGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGc                                                                                                                                 >  2:761673/1‑150 (MQ=255)
                      tCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGcc                                                                                                                                >  1:3803210/1‑150 (MQ=255)
                       cGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGcca                                                                                                                               >  1:2377651/1‑150 (MQ=255)
                         gCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGccacc                                                                                                                             <  2:2582400/150‑1 (MQ=255)
                          cATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGccacca                                                                                                                            >  1:2546600/1‑150 (MQ=255)
                           aTTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCcaccac                                                                                                                           >  2:1511324/1‑150 (MQ=255)
                               gCCCTTCCAGTTCTCATCCCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTg                                                                                                                       <  1:3349385/150‑1 (MQ=255)
                                 ccTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGat                                                                                                                     >  2:134746/1‑150 (MQ=255)
                                 ccTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGat                                                                                                                     >  2:2055335/1‑150 (MQ=255)
                                     ccAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATCTCg                                                                                                                 >  2:4238783/1‑150 (MQ=255)
                                     ccAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCg                                                                                                                 <  1:3083472/150‑1 (MQ=255)
                                        gTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCa                                                                                                              <  1:621477/150‑1 (MQ=255)
                                         ttCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCat                                                                                                             >  1:493072/1‑150 (MQ=255)
                                          tctcATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCatt                                                                                                            >  1:2647604/1‑150 (MQ=255)
                                            tcATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCattat                                                                                                          >  2:4467757/1‑150 (MQ=255)
                                             cATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCattatt                                                                                                         >  1:3844572/1‑150 (MQ=255)
                                                aCCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATg                                                                                                      >  2:3529063/1‑150 (MQ=255)
                                                   gTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGc                                                                                                   >  1:3059797/1‑150 (MQ=255)
                                                   gTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGc                                                                                                   <  2:3757701/150‑1 (MQ=255)
                                                     tGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCt                                                                                                 >  2:1618001/1‑150 (MQ=255)
                                                     tGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCt                                                                                                 <  1:1056768/150‑1 (MQ=255)
                                                      ggACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTg                                                                                                >  1:533457/1‑150 (MQ=255)
                                                       gACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTgg                                                                                               >  2:1731188/1‑150 (MQ=255)
                                                       gACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTgg                                                                                               <  2:22752/150‑1 (MQ=255)
                                                       gACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTgg                                                                                               >  1:1119813/1‑150 (MQ=255)
                                                          cAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCgg                                                                                            >  1:2346988/1‑150 (MQ=255)
                                                          cAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCgg                                                                                            >  1:3219246/1‑150 (MQ=255)
                                                            aTGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGt                                                                                          >  1:803313/1‑150 (MQ=255)
                                                              gAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTaa                                                                                        >  2:2485684/1‑150 (MQ=255)
                                                               aaaGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAAc                                                                                       >  2:3929021/1‑150 (MQ=255)
                                                                    gggTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGcc                                                                                  >  1:3047507/1‑150 (MQ=255)
                                                                        aaaCAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGAt                                                                              >  1:1587905/1‑150 (MQ=255)
                                                                         aaCAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATa                                                                             >  2:1448270/1‑150 (MQ=255)
                                                                           cAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATAcc                                                                           <  1:354120/150‑1 (MQ=255)
                                                                            aaTTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCg                                                                          <  2:1155265/150‑1 (MQ=255)
                                                                            aaTTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCg                                                                          <  2:1155220/150‑1 (MQ=255)
                                                                             aTTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGt                                                                         >  1:3649164/1‑150 (MQ=255)
                                                                             aTTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGt                                                                         >  1:1365242/1‑150 (MQ=255)
                                                                              tttACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTg                                                                        <  1:3138291/150‑1 (MQ=255)
                                                                              tttACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTg                                                                        >  1:1956628/1‑150 (MQ=255)
                                                                               ttACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTgg                                                                       >  1:1398827/1‑150 (MQ=255)
                                                                               ttACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTgg                                                                       >  1:4315797/1‑150 (MQ=255)
                                                                                  ccGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGata                                                                    <  1:394370/150‑1 (MQ=255)
                                                                                  ccGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGata                                                                    <  1:1196260/150‑1 (MQ=255)
                                                                                  ccGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGata                                                                    >  2:3415067/1‑150 (MQ=255)
                                                                                   cgacgaTCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                                                   >  2:2163096/1‑150 (MQ=255)
                                                                                   cgacgaTCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                                                   <  1:4273378/150‑1 (MQ=255)
                                                                                      cgaTCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                                                >  1:4492821/1‑150 (MQ=255)
                                                                                          cTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTcc                                                            <  1:2525712/150‑1 (MQ=255)
                                                                                          cTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTcc                                                            >  2:796484/1‑150 (MQ=255)
                                                                                           tGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCa                                                           >  2:4074567/1‑150 (MQ=255)
                                                                                             aTTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATc                                                         <  2:1245971/150‑1 (MQ=255)
                                                                                                 cccGTAACCTGCTTAAAGACGTGCCCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACTCTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGCTACCGTGGATATCGATTCCATccccc                                                     <  1:3899780/150‑2 (MQ=255)
                                                                                                   cGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGt                                                   >  2:1641593/1‑150 (MQ=255)
                                                                                                    gTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACAACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGtt                                                  >  2:488068/1‑150 (MQ=255)
                                                                                                      aaCCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGtttt                                                >  2:2782195/1‑150 (MQ=255)
                                                                                                        ccTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                                              <  2:860433/150‑1 (MQ=255)
                                                                                                         cTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAtt                                             <  2:4193313/150‑1 (MQ=255)
                                                                                                            cTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAg                                          >  1:991545/1‑150 (MQ=255)
                                                                                                            cTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAg                                          >  1:991587/1‑150 (MQ=255)
                                                                                                             ttAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAgg                                         >  2:1168041/1‑150 (MQ=255)
                                                                                                                aaGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGAtg                                      >  2:1870795/1‑150 (MQ=255)
                                                                                                                 aGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGAtgt                                     >  2:3669445/1‑150 (MQ=255)
                                                                                                                  gACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGAtgtg                                    <  1:3287452/150‑1 (MQ=255)
                                                                                                                   aCGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGa                                   <  2:1809200/150‑1 (MQ=255)
                                                                                                                       gACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATAc                               <  1:1703649/150‑1 (MQ=255)
                                                                                                                        aCCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACt                              <  2:3004242/150‑1 (MQ=255)
                                                                                                                          ccGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCa                            <  1:1442130/150‑1 (MQ=255)
                                                                                                                            gcgTGGTGGTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATg                          <  2:103141/150‑1 (MQ=255)
                                                                                                                              gtggtgGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGt                        >  1:3698262/1‑150 (MQ=255)
                                                                                                                               tggtggTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTc                       <  2:868249/150‑1 (MQ=255)
                                                                                                                               tggtggTTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTc                       <  1:558587/150‑1 (MQ=255)
                                                                                                                                ggtggtTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCg                      <  2:1927178/150‑1 (MQ=255)
                                                                                                                                 gtggtTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGc                     >  2:412559/1‑150 (MQ=255)
                                                                                                                                    gtTGACGTGTACGGCTAGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTg                  <  1:3906950/150‑1 (MQ=255)
                                                                                                                                      tGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGAt                >  2:454529/1‑150 (MQ=255)
                                                                                                                                              aCGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCaa        <  1:407526/150‑1 (MQ=255)
                                                                                                                                                 gctcgctcTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACg     >  2:3753925/1‑150 (MQ=255)
                                                                                                                                                 gctcgctcTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACg     <  2:3400999/150‑1 (MQ=255)
                                                                                                                                                 gctcgctcTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCACAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACg     <  1:3128989/150‑1 (MQ=255)
                                                                                                                                                  ctcgctcTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACgg    <  1:2887501/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGTATTCGATATTTTCAAAATCGGCATTGGCCCTTCCAGTTCTCATACCGTTGGACCAATGAAAGCGGGTAAACAATTTACCGACGATCTGATTGCCCGTAACCTGCTTAAAGACGTGACCCGCGTGGTGGTTGACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTCAGGATGTGAATACTCATGGTCGCCTGATGCTGGCAAACGGTC  >  CP009273/2922941‑2923238

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: