Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,474,993 G→T 19.1% 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,9930GT19.1% 809.1 / 118.0 269S259* (TCG→TAG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base G (117/100);  new base T (22/30);  total (139/130)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

CCGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTCAGCCAGC  >  CP009273/4474844‑4475142
                                                                                                                                                     |                                                                                                                                                     
ccGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCg                                                                                                                                                       >  2:3675098/1‑150 (MQ=255)
ccGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCg                                                                                                                                                       <  2:2870586/150‑1 (MQ=255)
ccGCGACCATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCg                                                                                                                                                       >  1:105041/1‑150 (MQ=255)
       cATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGtt                                                                                                                                                >  2:4117631/1‑150 (MQ=255)
       cATCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGtt                                                                                                                                                >  2:2790614/1‑150 (MQ=255)
        aTCCGCATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTg                                                                                                                                               >  1:1135091/1‑150 (MQ=255)
             cATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCttt                                                                                                                                          <  1:793586/150‑1 (MQ=255)
             cATCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCttt                                                                                                                                          <  1:65503/150‑1 (MQ=255)
               tCACGCCGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGt                                                                                                                                        >  1:515567/1‑150 (MQ=255)
                    ccGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCa                                                                                                                                   <  1:3658535/150‑1 (MQ=255)
                     cGTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCaa                                                                                                                                  <  1:4616291/150‑1 (MQ=255)
                      gTAAACCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAAt                                                                                                                                 >  1:1902746/1‑150 (MQ=255)
                         aaCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCAAAACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCAc                                                                                                                              >  1:3840528/1‑150 (MQ=255)
                          aCCGCTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCAcc                                                                                                                             <  1:2133619/150‑1 (MQ=255)
                              cTGCCGCACCGCACATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACa                                                                                                                         <  2:1290839/150‑1 (MQ=255)
                                   gcaccgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGc                                                                                                                    >  1:2328285/1‑150 (MQ=255)
                                   gcaccgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGc                                                                                                                    <  2:2920639/150‑1 (MQ=255)
                                     accgcacATATCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGa                                                                                                                  >  2:2361657/1‑150 (MQ=255)
                                            atatCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGttt                                                                                                           <  2:594262/150‑1 (MQ=255)
                                             tatCGTACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtttt                                                                                                          <  1:713090/150‑1 (MQ=255)
                                             tatCGTACGTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGtttt                                                                                                          <  1:2867655/150‑1 (MQ=255)
                                                  tACTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGaa                                                                                                     <  2:2260425/150‑1 (MQ=255)
                                                   aCTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAc                                                                                                    <  1:2508621/150‑1 (MQ=255)
                                                   aCTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAc                                                                                                    >  2:2475280/1‑150 (MQ=255)
                                                    cTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAcc                                                                                                   <  2:440934/150‑1 (MQ=255)
                                                    cTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAcc                                                                                                   <  2:3932417/150‑1 (MQ=255)
                                                    cTTCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAcc                                                                                                   >  2:2726445/1‑150 (MQ=255)
                                                     ttCATCTCATCCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAAccc                                                                                                  >  2:1205980/1‑150 (MQ=255)
                                                        atctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTga                                                                                               >  1:1185592/1‑150 (MQ=255)
                                                         tctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgac                                                                                              <  2:1634833/150‑1 (MQ=255)
                                                          ctcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgacc                                                                                             <  1:1358920/150‑1 (MQ=255)
                                                           tcatcCATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTgaccg                                                                                            <  2:4460195/150‑1 (MQ=255)
                                                               ccATGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGaccgaccg                                                                                        <  1:2283334/150‑1 (MQ=255)
                                                                 aTGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGcc                                                                                      >  2:3264454/1‑150 (MQ=255)
                                                                  tGCCTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCa                                                                                     >  1:3922663/1‑150 (MQ=255)
                                                                     cTTCTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGAt                                                                                  >  2:2076740/1‑150 (MQ=255)
                                                                        cTGAAGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAgg                                                                               <  2:3034454/150‑1 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            >  1:853210/1‑150 (MQ=255)
                                                                           aaGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAAt                                                                            >  1:3272185/1‑150 (MQ=255)
                                                                            aGGCTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATg                                                                           >  1:4431920/1‑150 (MQ=255)
                                                                               cTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCAt                                                                        >  1:2180876/1‑150 (MQ=255)
                                                                               cTTGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCAt                                                                        <  1:3008294/150‑1 (MQ=255)
                                                                                ttGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATc                                                                       >  2:2918086/1‑150 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  2:2206311/150‑1 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      >  2:2669164/1‑150 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  2:1972449/150‑1 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  1:1998351/150‑1 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  1:4019024/150‑1 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      <  1:4018865/150‑1 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      >  1:2734084/1‑150 (MQ=255)
                                                                                 tGATCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcc                                                                      >  2:3448931/1‑150 (MQ=255)
                                                                                   aTCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATcccc                                                                    >  2:1525780/1‑150 (MQ=255)
                                                                                    tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCa                                                                   <  2:2113732/150‑1 (MQ=255)
                                                                                    tCGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCa                                                                   <  1:2602519/150‑1 (MQ=255)
                                                                                     cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                  >  2:647678/1‑150 (MQ=255)
                                                                                     cGAGATACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAg                                                                  >  1:2026303/1‑150 (MQ=255)
                                                                                        gaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctc                                                               >  2:4355336/1‑150 (MQ=255)
                                                                                        gaTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctc                                                               >  2:422325/1‑150 (MQ=255)
                                                                                        gaTACCGCCGGAGTCGAAGGTCAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctc                                                               >  2:1918453/1‑150 (MQ=255)
                                                                                         aTACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGctct                                                              >  2:4122059/1‑150 (MQ=255)
                                                                                          tACCGCCGGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCtt                                                             <  2:1087752/150‑1 (MQ=255)
                                                                                               ccgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATc                                                        >  1:1450456/1‑150 (MQ=255)
                                                                                               ccgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATc                                                        <  2:3512526/150‑1 (MQ=255)
                                                                                                cgGAGTCGAAGGTTAAACCTTTACCCACCAGCTCGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATct                                                       <  2:1562966/150‑1 (MQ=255)
                                                                                                cgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATct                                                       >  1:2808972/1‑150 (MQ=255)
                                                                                                cgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATct                                                       >  1:2240390/1‑150 (MQ=255)
                                                                                                cgGAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATct                                                       >  2:381055/1‑150 (MQ=255)
                                                                                                 ggAGTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctg                                                      >  2:20457/1‑150 (MQ=255)
                                                                                                    gTCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATctgctg                                                   >  2:767619/1‑150 (MQ=255)
                                                                                                     tCGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGt                                                  >  1:755775/1‑150 (MQ=255)
                                                                                                      cGAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGtt                                                 <  1:4083502/150‑1 (MQ=255)
                                                                                                       gAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTc                                                >  1:1378252/1‑150 (MQ=255)
                                                                                                       gAAGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTc                                                <  2:741385/150‑1 (MQ=255)
                                                                                                        aaGGTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCg                                               <  2:4338006/150‑1 (MQ=255)
                                                                                                           gTTAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCg                                            <  1:4419592/150‑1 (MQ=255)
                                                                                                            ttAAACCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGa                                           >  2:2115773/1‑150 (MQ=255)
                                                                                                              aaaCCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATa                                         >  1:2951350/1‑150 (MQ=255)
                                                                                                              aaaCCTTTACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATa                                         <  2:1731313/150‑1 (MQ=255)
                                                                                                                    ttACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCgg                                   <  2:1488918/150‑1 (MQ=255)
                                                                                                                    ttACCCACCAGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCgg                                   >  1:1749426/1‑150 (MQ=255)
                                                                                                                         caccaGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGtgat                              <  1:3465114/150‑1 (MQ=255)
                                                                                                                            caGCACGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGAc                           <  2:2598639/150‑1 (MQ=255)
                                                                                                                             aGCACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACa                          >  1:2043994/1‑150 (MQ=255)
                                                                                                                               cACGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACAtt                        <  1:3800804/150‑1 (MQ=255)
                                                                                                                                aCGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTc                       <  2:2027071/150‑1 (MQ=255)
                                                                                                                                aCGATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTc                       <  1:3155368/150‑1 (MQ=255)
                                                                                                                                 cGATTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCt                      <  2:1827154/150‑1 (MQ=255)
                                                                                                                                  gATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCtt                     <  2:1619706/150‑1 (MQ=255)
                                                                                                                                  gATTGGGCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCtt                     <  1:705807/150‑1 (MQ=255)
                                                                                                                                   aTTGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTg                    >  2:2500351/1‑150 (MQ=255)
                                                                                                                                     tGGGCGTGCATCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCt                  <  1:788260/150‑1 (MQ=255)
                                                                                                                                      gggCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTg                 >  2:10675/1‑150 (MQ=255)
                                                                                                                                       ggCGTGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGt                <  1:4149287/150‑1 (MQ=255)
                                                                                                                                           tGCATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCt            <  1:3326553/150‑1 (MQ=255)
                                                                                                                                             cATCTTCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGt          <  2:515567/150‑1 (MQ=255)
                                                                                                                                              aTCTTCCTACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTc         <  1:2682845/150‑1 (MQ=255)
                                                                                                                                                 ttCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagc      <  1:811855/150‑1 (MQ=255)
                                                                                                                                                  tCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagcc     >  1:205713/1‑150 (MQ=255)
                                                                                                                                                  tCCGACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagcc     >  2:2296512/1‑150 (MQ=255)
                                                                                                                                                     gACGCGTTGCCTTTGTACTCAATCACCGACATCAGCGATTCGTTTTGCGAACCCTGACCGACCGCCAGATAGGAATGCATCCCCAGCTCTTTCATCTGCTGTTCGCCGATAACGCGGGTGATGACATTCTTGCTGTAGCTGTcagccagc  >  1:620334/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: