Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,474,993 G→T 41.7% S259* (TCG→TAG)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,474,9930GT41.7% 284.0 / 304.9 256S259* (TCG→TAG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base G (79/70);  new base T (53/54);  total (132/124)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

CGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTCAGCCAG  >  CP009273/4474845‑4475141
                                                                                                                                                    |                                                                                                                                                    
cgcgACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGa                                                                                                                                                     <  1:3222009/150‑1 (MQ=255)
       aTCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTg                                                                                                                                              >  1:1518881/1‑150 (MQ=255)
       aTCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTg                                                                                                                                              <  1:164256/150‑1 (MQ=255)
       aTCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTg                                                                                                                                              >  2:124969/1‑150 (MQ=255)
         ccGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGcc                                                                                                                                            >  1:4624733/1‑150 (MQ=255)
          cGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCt                                                                                                                                           >  2:2922555/1‑150 (MQ=255)
          cGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCt                                                                                                                                           >  2:384123/1‑150 (MQ=255)
            cATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCttt                                                                                                                                         <  1:4153550/150‑1 (MQ=255)
             aTCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTg                                                                                                                                        <  1:2923588/150‑1 (MQ=255)
                     gTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAAt                                                                                                                                >  1:1119505/1‑150 (MQ=255)
                      tAAACCGCTGCCGCACCGCACATATCGTACTTCATTTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATc                                                                                                                               <  2:42061/150‑1 (MQ=255)
                      tAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATc                                                                                                                               >  1:1226131/1‑150 (MQ=255)
                      tAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATc                                                                                                                               <  2:3038778/150‑1 (MQ=255)
                       aaaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCa                                                                                                                              >  2:1115567/1‑150 (MQ=255)
                            gCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGAc                                                                                                                         <  2:2534367/150‑1 (MQ=255)
                             cTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACa                                                                                                                        >  2:677799/1‑150 (MQ=255)
                             cTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACa                                                                                                                        >  2:3859997/1‑150 (MQ=255)
                             cTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACa                                                                                                                        >  2:1650728/1‑150 (MQ=255)
                              tGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACAt                                                                                                                       <  1:3476575/150‑1 (MQ=255)
                                 cgcaccgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAg                                                                                                                    <  1:843725/150‑1 (MQ=255)
                                 cgcaccgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAg                                                                                                                    <  2:1942604/150‑1 (MQ=255)
                                      cgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGAtt                                                                                                               <  2:2996953/150‑1 (MQ=255)
                                         acaTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGt                                                                                                            >  1:2057451/1‑150 (MQ=255)
                                          caTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGATATACCGCCGGAGTCGAAGGTTAAACCTATACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtt                                                                                                           >  2:2256948/1‑150 (MQ=255)
                                          caTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtt                                                                                                           <  2:1970525/150‑1 (MQ=255)
                                          caTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtt                                                                                                           <  1:3521470/150‑1 (MQ=255)
                                          caTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtt                                                                                                           >  2:1494011/1‑150 (MQ=255)
                                            tatCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtttt                                                                                                         <  2:4291622/150‑1 (MQ=255)
                                             atCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTg                                                                                                        <  2:2371733/150‑1 (MQ=255)
                                              tCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGc                                                                                                       >  2:3081325/1‑150 (MQ=255)
                                               cGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCg                                                                                                      >  2:2458785/1‑150 (MQ=255)
                                                gTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGa                                                                                                     >  2:301622/1‑150 (MQ=255)
                                                  aCTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAc                                                                                                   >  1:1263632/1‑150 (MQ=255)
                                                   cTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAcc                                                                                                  <  2:3144866/150‑1 (MQ=255)
                                                    ttCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAccc                                                                                                 >  2:4419641/1‑150 (MQ=255)
                                                    ttCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAccc                                                                                                 >  1:2768145/1‑150 (MQ=255)
                                                      catctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTg                                                                                               <  2:1722766/150‑1 (MQ=255)
                                                      catctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTg                                                                                               <  1:1699817/150‑1 (MQ=255)
                                                       atctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTga                                                                                              <  2:16377/150‑1 (MQ=255)
                                                        tctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgac                                                                                             <  2:4431827/150‑1 (MQ=255)
                                                            atcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                                                                                         >  1:2184623/1‑150 (MQ=255)
                                                            atcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                                                                                         >  2:2662086/1‑150 (MQ=255)
                                                             tcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgacc                                                                                        <  2:715795/150‑1 (MQ=255)
                                                               cATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGc                                                                                      >  2:3181197/1‑150 (MQ=255)
                                                                 tGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCa                                                                                    <  2:3220468/150‑1 (MQ=255)
                                                                  gCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAg                                                                                   >  1:2700453/1‑150 (MQ=255)
                                                                   ccTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGa                                                                                  >  2:3425386/1‑150 (MQ=255)
                                                                   ccTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGa                                                                                  >  2:1062579/1‑150 (MQ=255)
                                                                    cTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGAt                                                                                 <  1:3561332/150‑1 (MQ=255)
                                                                    cTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGAt                                                                                 <  1:3561828/150‑1 (MQ=255)
                                                                     ttCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATa                                                                                <  1:3249522/150‑1 (MQ=255)
                                                                      tCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAg                                                                               >  2:3768128/1‑150 (MQ=255)
                                                                       cTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAgg                                                                              >  1:4433951/1‑150 (MQ=255)
                                                                          aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                           >  1:651774/1‑150 (MQ=255)
                                                                              cTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCAt                                                                       <  1:4506059/150‑1 (MQ=255)
                                                                                 gATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc                                                                    >  1:2843934/1‑150 (MQ=255)
                                                                                 gATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc                                                                    >  2:2069640/1‑150 (MQ=255)
                                                                                   tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCa                                                                  <  2:3524249/150‑1 (MQ=255)
                                                                                   tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCa                                                                  <  1:412191/150‑1 (MQ=255)
                                                                                   tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCa                                                                  <  1:2923962/150‑1 (MQ=255)
                                                                                    cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                 >  1:4179545/1‑150 (MQ=255)
                                                                                     gagaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGc                                                                >  1:1045516/1‑150 (MQ=255)
                                                                                      agaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGct                                                               >  2:622819/1‑150 (MQ=255)
                                                                                          aCCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCttt                                                           >  2:897072/1‑150 (MQ=255)
                                                                                                ggAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctg                                                     >  2:3697971/1‑150 (MQ=255)
                                                                                                 gAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctgc                                                    >  1:2396982/1‑150 (MQ=255)
                                                                                                   gTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctgctg                                                  <  2:803260/150‑1 (MQ=255)
                                                                                                        aGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGc                                             >  1:2568754/1‑150 (MQ=255)
                                                                                                         ggTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGcc                                            >  1:4580320/1‑150 (MQ=255)
                                                                                                               aCCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAc                                      <  1:3061259/150‑1 (MQ=255)
                                                                                                                ccTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcg                                     <  2:4062272/150‑1 (MQ=255)
                                                                                                                ccTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcg                                     <  2:2193882/150‑1 (MQ=255)
                                                                                                                 cTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcgc                                    >  2:2692732/1‑150 (MQ=255)
                                                                                                                 cTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcgc                                    >  2:3459995/1‑150 (MQ=255)
                                                                                                                   ttACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCgg                                  <  1:970813/150‑1 (MQ=255)
                                                                                                                    tACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCggg                                 <  1:2363502/150‑1 (MQ=255)
                                                                                                                    tACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCggg                                 <  1:220591/150‑1 (MQ=255)
                                                                                                                     aCCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGt                                >  1:3761571/1‑150 (MQ=255)
                                                                                                                      cccACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGtg                               <  1:1402835/150‑1 (MQ=255)
                                                                                                                            aGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACa                         >  1:1013233/1‑150 (MQ=255)
                                                                                                                             gCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACAt                        >  2:2963174/1‑150 (MQ=255)
                                                                                                                              cACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACAtt                       <  1:3650231/150‑1 (MQ=255)
                                                                                                                              cACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACAtt                       <  1:1464454/150‑1 (MQ=255)
                                                                                                                                   ttGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGc                  <  1:4087428/150‑1 (MQ=255)
                                                                                                                                    tGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCt                 >  2:4512290/1‑150 (MQ=255)
                                                                                                                                      ggCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGt               <  1:1696377/150‑1 (MQ=255)
                                                                                                                                        cGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAg             <  2:2123992/150‑1 (MQ=255)
                                                                                                                                        cGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAg             >  2:480298/1‑150 (MQ=255)
                                                                                                                                         gTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGc            >  2:2761139/1‑150 (MQ=255)
                                                                                                                                         gTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGc            >  1:2052749/1‑150 (MQ=255)
                                                                                                                                         gTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGc            >  2:609519/1‑150 (MQ=255)
                                                                                                                                         gTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGAAATTCTTGCTGTAGc            >  2:1354009/1‑150 (MQ=255)
                                                                                                                                          tGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCt           <  1:2700669/150‑1 (MQ=255)
                                                                                                                                            cATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGt         <  1:497589/150‑1 (MQ=255)
                                                                                                                                             aTCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTc        >  2:1020534/1‑150 (MQ=255)
                                                                                                                                             aTCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTc        >  2:3745653/1‑150 (MQ=255)
                                                                                                                                               cTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcag      >  1:392618/1‑150 (MQ=255)
                                                                                                                                                ttCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagc     >  2:3338737/1‑150 (MQ=255)
                                                                                                                                                ttCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagc     >  2:143674/1‑150 (MQ=255)
                                                                                                                                                  ccGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagcca   <  2:3673508/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTCAGCCAG  >  CP009273/4474845‑4475141

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: