Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,474,993 G→T 44.5% 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,9930GT44.5% 221.8 / inf 261S259* (TCG→TAG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base G (74/71);  new base T (70/46);  total (144/117)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 264 total aligned reads displayed.

CCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTCAGCCAGC  >  CP009273/4474844‑4475142
                                                                                                                                                     |                                                                                                                                                     
ccGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCg                                                                                                                                                       <  1:1423376/150‑1 (MQ=255)
    gACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTAcgc                                                                                                                                                   >  1:2659095/1‑150 (MQ=255)
     aCCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGAcgcg                                                                                                                                                  <  2:3403740/150‑1 (MQ=255)
      ccATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGt                                                                                                                                                 <  1:3362271/150‑1 (MQ=255)
              aTCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTg                                                                                                                                         >  1:3488573/1‑150 (MQ=255)
              aTCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTg                                                                                                                                         >  1:3119297/1‑150 (MQ=255)
              aTCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTg                                                                                                                                         >  1:2432303/1‑150 (MQ=255)
               tCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGt                                                                                                                                        <  2:842416/150‑1 (MQ=255)
                  cGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACt                                                                                                                                     <  2:2717553/150‑1 (MQ=255)
                  cGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACt                                                                                                                                     <  2:2862059/150‑1 (MQ=255)
                  cGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACt                                                                                                                                     >  2:2682107/1‑150 (MQ=255)
                  cGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTACTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTAGTACt                                                                                                                                     >  2:2061575/1‑150 (MQ=255)
                     cGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCaa                                                                                                                                  >  2:477581/1‑150 (MQ=255)
                      gTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAAt                                                                                                                                 >  1:3064908/1‑150 (MQ=255)
                      gTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAAt                                                                                                                                 >  1:4005247/1‑150 (MQ=255)
                       tAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATc                                                                                                                                <  1:2110293/150‑1 (MQ=255)
                        aaaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCa                                                                                                                               <  2:547862/150‑1 (MQ=255)
                        aaaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGATGCCTTTGTACTCAATCa                                                                                                                               >  1:37886/1‑150 (MQ=255)
                         aaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCAc                                                                                                                              >  1:1783682/1‑150 (MQ=255)
                         aaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCAc                                                                                                                              >  1:2540532/1‑150 (MQ=255)
                          aCCGCTGCCGCACCGCACATATCGTACTTCATCTCTTCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCAcc                                                                                                                             <  2:4173043/150‑1 (MQ=255)
                          aCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCAcc                                                                                                                             <  2:3557806/150‑1 (MQ=255)
                               tGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACAt                                                                                                                        <  1:3012737/150‑1 (MQ=255)
                                gCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATc                                                                                                                       <  1:1453683/150‑1 (MQ=255)
                                   gcaccgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGc                                                                                                                    <  1:2699442/150‑1 (MQ=255)
                                       cgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGAtt                                                                                                                >  1:2326315/1‑150 (MQ=255)
                                         cacaTATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCg                                                                                                              >  2:830862/1‑150 (MQ=255)
                                           caTATCGTACTTCATCTCATCCATGCCTTCNGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtt                                                                                                            >  1:3858063/1‑150 (MQ=255)
                                            atatCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATAACGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGttt                                                                                                           >  2:1634987/1‑150 (MQ=255)
                                             tatCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTCCGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtttt                                                                                                          <  1:3769505/150‑1 (MQ=255)
                                                      tcatctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCt                                                                                                 >  1:2777550/1‑150 (MQ=255)
                                                       catctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTg                                                                                                <  2:1849580/150‑1 (MQ=255)
                                                          ctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgacc                                                                                             <  1:1311782/150‑1 (MQ=255)
                                                          ctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgacc                                                                                             <  1:1444210/150‑1 (MQ=255)
                                                             atcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgac                                                                                          <  1:969996/150‑1 (MQ=255)
                                                              tcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccgacc                                                                                         >  1:1889849/1‑150 (MQ=255)
                                                               ccATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                                                                                        <  2:1560924/150‑1 (MQ=255)
                                                                cATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGc                                                                                       <  2:3323871/150‑1 (MQ=255)
                                                                cATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGc                                                                                       >  2:2635792/1‑150 (MQ=255)
                                                                   gCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAg                                                                                    <  1:10284/150‑1 (MQ=255)
                                                                      ttCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATa                                                                                 <  1:1479115/150‑1 (MQ=255)
                                                                         tGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGa                                                                              >  1:2095790/1‑150 (MQ=255)
                                                                          gAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGaa                                                                             >  1:1151573/1‑150 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            <  1:2524141/150‑1 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            <  1:2581887/150‑1 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            <  1:418671/150‑1 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            <  2:3035292/150‑1 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGACAGATAGGAAt                                                                            <  1:2671041/150‑1 (MQ=255)
                                                                            aGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATg                                                                           <  1:1327594/150‑1 (MQ=255)
                                                                             ggCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGc                                                                          >  2:670869/1‑150 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  2:3666692/150‑1 (MQ=255)
                                                                                  gATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATccc                                                                     >  2:699660/1‑150 (MQ=255)
                                                                                   aTCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcccc                                                                    >  2:1494868/1‑150 (MQ=255)
                                                                                    tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCCGATAGGAATGCATCCCCa                                                                   >  2:4006591/1‑150 (MQ=255)
                                                                                     cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                  >  2:2527892/1‑150 (MQ=255)
                                                                                     cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                  >  2:2527815/1‑150 (MQ=255)
                                                                                     cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                  >  1:1391494/1‑150 (MQ=255)
                                                                                       agaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGct                                                                <  2:610207/150‑1 (MQ=255)
                                                                                        gaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctc                                                               <  1:2274302/150‑1 (MQ=255)
                                                                                        gaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCAACAGCACGATTGTGCGTGCATCTTCCTTCGCGTTGCCTTTGTACTCATTCACCGACATCAGCGATTCGTTTTTCGAAACCTGACCGTCAGCCAGATAGGAATGCATCCCAAGcta                                                               <  2:689274/150‑2 (MQ=255)
                                                                                         aTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctct                                                              >  2:1271213/1‑150 (MQ=255)
                                                                                            ccgccgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTc                                                           <  2:1100690/150‑1 (MQ=255)
                                                                                                   aGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctgct                                                    >  1:2206463/1‑150 (MQ=255)
                                                                                                   aGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctgct                                                    >  1:367123/1‑150 (MQ=255)
                                                                                                     tCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGACTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGt                                                  >  1:1768136/1‑150 (MQ=255)
                                                                                                        aaGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCg                                               <  1:156381/150‑1 (MQ=255)
                                                                                                          ggTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGcc                                             >  2:197437/1‑150 (MQ=255)
                                                                                                           gTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCg                                            >  2:1256199/1‑150 (MQ=255)
                                                                                                           gTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCg                                            >  1:698359/1‑150 (MQ=255)
                                                                                                            ttAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGa                                           >  1:3474571/1‑150 (MQ=255)
                                                                                                             tAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGAt                                          >  2:3373342/1‑150 (MQ=255)
                                                                                                             tAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGAt                                          >  1:553590/1‑150 (MQ=255)
                                                                                                              aaaCCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATa                                         >  2:3038810/1‑150 (MQ=255)
                                                                                                                aCCTTTACCCACCAGCACGATTGGGCGTGCTTCTTCAGACGCGTTGCCTTTGTAATCAATCACCGACATCAGCGATTCGTTTTGCGATCCCTGTCCGACCGCAAGATAGGAATGCATCCCAAGCTCTTTCATCTGCTGTTCGCAGATAAc                                       <  2:1028157/150‑1 (MQ=255)
                                                                                                                 ccTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcg                                      >  2:3216770/1‑150 (MQ=255)
                                                                                                                  cTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcgc                                     >  1:1945962/1‑150 (MQ=255)
                                                                                                                  cTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcgc                                     >  1:902120/1‑150 (MQ=255)
                                                                                                                   tttACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAAcgcg                                    <  1:4085861/150‑1 (MQ=255)
                                                                                                                       cccACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGtg                                <  2:201560/150‑1 (MQ=255)
                                                                                                                        ccaccaGCACGATTTGGCGTGAATCTTCCGACGCGTTGCAATTGTACACAATCACCGACATCAGCGATAAGTAATGCTAACACTGACCGACCGCCAGATAGGAATGCATCACCAGCTCTTTCATCTGCTGATCGCCCATAACGAAGGtga                               >  2:3871582/1‑150 (MQ=255)
                                                                                                                        ccaccaGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGtga                               >  2:773660/1‑150 (MQ=255)
                                                                                                                        ccaccaGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGtga                               >  2:1891166/1‑150 (MQ=255)
                                                                                                                            caGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGAc                           >  2:1037675/1‑150 (MQ=255)
                                                                                                                            caGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGAc                           >  2:1917559/1‑150 (MQ=255)
                                                                                                                             aGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACa                          >  1:1670050/1‑150 (MQ=255)
                                                                                                                               cACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACAtt                        >  1:2855049/1‑150 (MQ=255)
                                                                                                                                aCGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTc                       >  1:23472/1‑150 (MQ=255)
                                                                                                                                aCGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTc                       >  2:418518/1‑150 (MQ=255)
                                                                                                                                aCGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTc                       >  1:3260756/1‑150 (MQ=255)
                                                                                                                                 cGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCt                      >  2:3918175/1‑150 (MQ=255)
                                                                                                                                  gATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTGTCATCTGCTGTGCGCCGCTAACGCGGGTGATGACATTCtt                     >  2:2509130/1‑150 (MQ=255)
                                                                                                                                     tGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCt                  >  2:2747080/1‑150 (MQ=255)
                                                                                                                                         cGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAg              >  2:181482/1‑150 (MQ=255)
                                                                                                                                          gTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGc             >  2:3085693/1‑150 (MQ=255)
                                                                                                                                           tGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCt            >  2:87023/1‑150 (MQ=255)
                                                                                                                                           tGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCt            <  1:439635/150‑1 (MQ=255)
                                                                                                                                            gCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTg           <  2:3717455/150‑1 (MQ=255)
                                                                                                                                             cATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGt          >  2:2572346/1‑150 (MQ=255)
                                                                                                                                              aTCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTc         >  2:4129462/1‑150 (MQ=255)
                                                                                                                                                     gACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagccagc  >  2:2346669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTCAGCCAGC  >  CP009273/4474844‑4475142

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: