Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,030,073 C→A 13.5% E465* (GAA→TAA)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,030,0730CA13.5% 1031.1 / 97.4 311E465* (GAA→TAA) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (122/146);  new base A (18/25);  total (140/171)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

CGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCT  >  CP009273/3029926‑3030221
                                                                                                                                                   |                                                                                                                                                    
cGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGc                                                                                                                                                    >  1:3310464/1‑150 (MQ=255)
 gCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGcc                                                                                                                                                   <  2:1279400/150‑1 (MQ=255)
 gCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGcc                                                                                                                                                   <  2:2065924/150‑1 (MQ=255)
  cccACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCt                                                                                                                                                  <  1:2803311/150‑1 (MQ=255)
  cccACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTAGCCt                                                                                                                                                  <  1:2565187/150‑1 (MQ=255)
    caccaGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTg                                                                                                                                                >  1:2966047/1‑150 (MQ=255)
    caccaGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTg                                                                                                                                                >  2:2892736/1‑150 (MQ=255)
       caGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCaa                                                                                                                                             >  2:3353937/1‑150 (MQ=255)
        aGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAAt                                                                                                                                            >  2:2925707/1‑150 (MQ=255)
        aGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAAt                                                                                                                                            <  2:1848822/150‑1 (MQ=255)
           cGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGc                                                                                                                                         <  2:535433/150‑1 (MQ=255)
           cGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTCGCAATAGc                                                                                                                                         <  1:3813929/150‑1 (MQ=255)
            gCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCg                                                                                                                                        <  2:1697691/150‑1 (MQ=255)
              tGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGaa                                                                                                                                      >  2:2874867/1‑150 (MQ=255)
              tGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGaa                                                                                                                                      >  1:2264247/1‑150 (MQ=255)
              tGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGaa                                                                                                                                      <  2:3429956/150‑1 (MQ=255)
              tGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGaa                                                                                                                                      >  1:2498551/1‑150 (MQ=255)
               gTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAg                                                                                                                                     <  2:452090/150‑1 (MQ=255)
               gTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAg                                                                                                                                     >  1:2680120/1‑150 (MQ=255)
                 tGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAggg                                                                                                                                   >  1:4529944/1‑150 (MQ=255)
                     gcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCa                                                                                                                               <  2:3572108/150‑1 (MQ=255)
                     gcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCa                                                                                                                               <  1:669334/150‑1 (MQ=255)
                     gcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCa                                                                                                                               <  1:2913227/150‑1 (MQ=255)
                      cagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAg                                                                                                                              <  1:711559/150‑1 (MQ=255)
                      cagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAg                                                                                                                              <  1:1922900/150‑1 (MQ=255)
                      cagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAg                                                                                                                              >  2:657282/1‑150 (MQ=255)
                         aCGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCa                                                                                                                           <  1:3246975/150‑1 (MQ=255)
                           ggAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACt                                                                                                                         <  1:1798164/150‑1 (MQ=255)
                             aaaTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCa                                                                                                                       >  2:398020/1‑150 (MQ=255)
                             aaaTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCa                                                                                                                       <  1:3398978/150‑1 (MQ=255)
                             aaaTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCa                                                                                                                       >  1:167474/1‑150 (MQ=255)
                              aaTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCCATAGCGAAGGGTCCAGCCACTCAg                                                                                                                      >  2:3637810/1‑150 (MQ=255)
                              aaTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAg                                                                                                                      <  1:1320371/150‑1 (MQ=255)
                               aTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGt                                                                                                                     <  1:3993551/150‑1 (MQ=255)
                               aTGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGt                                                                                                                     <  1:4122028/150‑1 (MQ=255)
                                 gACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTaa                                                                                                                   <  2:2205062/150‑1 (MQ=255)
                                      tCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGt                                                                                                              >  1:101719/1‑150 (MQ=255)
                                       cAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGtt                                                                                                             <  1:4083697/150‑1 (MQ=255)
                                         aaCAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCg                                                                                                           >  1:533455/1‑150 (MQ=255)
                                         aaCAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCg                                                                                                           >  2:3399588/1‑150 (MQ=255)
                                           cAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGcc                                                                                                         >  1:3653034/1‑150 (MQ=255)
                                            aGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCa                                                                                                        <  1:3477129/150‑1 (MQ=255)
                                            aGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCa                                                                                                        >  1:2206182/1‑150 (MQ=255)
                                             gCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCaa                                                                                                       <  1:567214/150‑1 (MQ=255)
                                                    cGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCgttg                                                                                                <  1:3055814/150‑1 (MQ=255)
                                                     gggAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCgttgt                                                                                               <  2:4237245/150‑1 (MQ=255)
                                                       gAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGttgttg                                                                                             >  2:4481250/1‑150 (MQ=255)
                                                          cATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGaaa                                                                                          >  2:234446/1‑150 (MQ=255)
                                                           aTCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAag                                                                                         >  1:1502390/1‑150 (MQ=255)
                                                            tCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAaga                                                                                        <  2:1416631/150‑1 (MQ=255)
                                                              tGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGt                                                                                      <  1:1429968/150‑1 (MQ=255)
                                                                  cccACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGa                                                                                  <  1:3233860/150‑1 (MQ=255)
                                                                   ccACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGaa                                                                                 <  2:1870524/150‑1 (MQ=255)
                                                                    cACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAAt                                                                                <  2:1212064/150‑1 (MQ=255)
                                                                    cACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCATCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAAt                                                                                <  2:2278233/150‑1 (MQ=255)
                                                                      cGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAAttt                                                                              <  2:3251932/150‑1 (MQ=255)
                                                                       gggCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTa                                                                             >  1:905817/1‑150 (MQ=255)
                                                                           cAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCct                                                                         >  2:2785309/1‑150 (MQ=255)
                                                                           cAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCct                                                                         >  2:3603212/1‑150 (MQ=255)
                                                                           cAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCct                                                                         <  2:2527209/150‑1 (MQ=255)
                                                                            aGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCctc                                                                        <  2:1517849/150‑1 (MQ=255)
                                                                              cATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCtt                                                                      >  2:1178256/1‑150 (MQ=255)
                                                                                     cagcagCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGa                                                               >  1:385620/1‑150 (MQ=255)
                                                                                     cagcagCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGa                                                               >  1:3934932/1‑150 (MQ=255)
                                                                                       gcagcTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACa                                                             <  2:3755992/150‑1 (MQ=255)
                                                                                         agcTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACaaa                                                           >  2:1962327/1‑150 (MQ=255)
                                                                                            tGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAAccc                                                        <  1:3617215/150‑1 (MQ=255)
                                                                                               gcCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACccgcc                                                     >  2:1705582/1‑150 (MQ=255)
                                                                                                ccACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACccgccg                                                    >  2:2122393/1‑150 (MQ=255)
                                                                                                ccACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACccgccg                                                    >  2:1555345/1‑150 (MQ=255)
                                                                                                         tGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCAt                                           <  1:234455/150‑1 (MQ=255)
                                                                                                              aTTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCg                                      >  2:2400308/1‑150 (MQ=255)
                                                                                                                ttCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGcc                                    <  1:1930750/150‑1 (MQ=255)
                                                                                                                     ccGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACt                               <  1:2122393/150‑1 (MQ=255)
                                                                                                                      cGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACtt                              <  2:3224514/150‑1 (MQ=255)
                                                                                                                       gggCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTc                             >  1:4178276/1‑150 (MQ=255)
                                                                                                                       gggCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTc                             >  1:1185005/1‑150 (MQ=255)
                                                                                                                        ggCTTAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCa                            <  2:824142/150‑1 (MQ=255)
                                                                                                                         gCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAg                           >  2:2660747/1‑150 (MQ=255)
                                                                                                                         gCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAg                           <  2:3468768/150‑1 (MQ=255)
                                                                                                                            tAACGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGcat                        >  2:3145335/1‑150 (MQ=255)
                                                                                                                             aaCGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGcatc                       >  2:679106/1‑150 (MQ=255)
                                                                                                                             aaCGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGcatc                       <  2:4328988/150‑1 (MQ=255)
                                                                                                                               cGGACCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGcatcat                     <  2:2236403/150‑1 (MQ=255)
                                                                                                                                ggACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATg                    <  2:1547824/150‑1 (MQ=255)
                                                                                                                                ggACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATg                    >  1:1416538/1‑150 (MQ=255)
                                                                                                                                 gACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGc                   <  2:4001233/150‑1 (MQ=255)
                                                                                                                                  aCCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGcc                  <  2:4083982/150‑1 (MQ=255)
                                                                                                                                  aCCGTCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGcc                  <  2:1018775/150‑1 (MQ=255)
                                                                                                                                   ccGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCa                 >  1:638270/1‑150 (MQ=255)
                                                                                                                                     gTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAgg               >  1:521852/1‑150 (MQ=255)
                                                                                                                                      tCTGATACCACTTAGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAggg              <  2:2103397/150‑1 (MQ=255)
                                                                                                                                        tGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTa            >  1:4352980/1‑150 (MQ=255)
                                                                                                                                        tGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTa            >  1:374792/1‑150 (MQ=255)
                                                                                                                                          aTACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTaga          >  2:2547533/1‑150 (MQ=255)
                                                                                                                                          aTACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTaga          >  1:2853244/1‑150 (MQ=255)
                                                                                                                                             ccACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAgtg       >  1:1051838/1‑150 (MQ=255)
                                                                                                                                              cACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAgtgt      <  2:703945/150‑1 (MQ=255)
                                                                                                                                              cACTTCGCCGTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAgtgt      <  2:3310464/150‑1 (MQ=255)
                                                                                                                                                 ttCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTc   <  1:1094050/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCT  >  CP009273/3029926‑3030221

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: