Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,961,385 C→A 15.9% E137* (GAA→TAA)  ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,961,3850CA15.9% 1102.0 / 140.2 353E137* (GAA→TAA) ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (121/176);  new base A (20/36);  total (141/212)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 362 total aligned reads displayed.

CCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAA  >  CP009273/2961236‑2961534
                                                                                                                                                     |                                                                                                                                                     
ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTc                                                                                                                                                       <  1:3367596/150‑1 (MQ=255)
ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTACCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTc                                                                                                                                                       <  2:3512392/150‑1 (MQ=255)
 cAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                                                                                                                      >  1:2465991/1‑150 (MQ=255)
 cAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                                                                                                                      >  1:69665/1‑150 (MQ=255)
 cAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                                                                                                                      >  1:836940/1‑150 (MQ=255)
   gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTc                                                                                                                                                    <  2:135980/150‑1 (MQ=255)
   gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTc                                                                                                                                                    <  1:2785399/150‑1 (MQ=255)
      ttCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATa                                                                                                                                                 <  2:2924968/150‑1 (MQ=255)
        cGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACt                                                                                                                                               >  2:3129319/1‑150 (MQ=255)
        cGGCAATCGCCACACCAGGTGCGGCCGGGAATGCGCGGAGTCGCGTCTGGCGATATTGCCCAAACCCGGCAGTCAACTGCGACTGAGAAAGAATAGCGGCCCTCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACt                                                                                                                                               <  2:3537325/150‑1 (MQ=255)
         ggCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTg                                                                                                                                              <  2:4197796/150‑1 (MQ=255)
         ggCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTg                                                                                                                                              <  2:4493110/150‑1 (MQ=255)
           cAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGgc                                                                                                                                            >  1:3065464/1‑150 (MQ=255)
            aaTCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                                                                                                                           >  2:4359502/1‑150 (MQ=255)
              tCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCa                                                                                                                                         <  1:3002755/150‑1 (MQ=255)
                gCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAc                                                                                                                                       <  1:2228437/150‑1 (MQ=255)
                    cacCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACtctc                                                                                                                                   <  1:1448859/150‑1 (MQ=255)
                    cacCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACtctc                                                                                                                                   <  2:4255761/150‑1 (MQ=255)
                     acCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACACGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCg                                                                                                                                  <  1:4581408/150‑1 (MQ=255)
                     acCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCg                                                                                                                                  >  1:1352850/1‑150 (MQ=255)
                       cAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGtt                                                                                                                                <  1:3765844/150‑1 (MQ=255)
                        aGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAGCTCTCGTTg                                                                                                                               >  1:1072827/1‑150 (MQ=255)
                             ctgcCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTa                                                                                                                          <  1:2502321/150‑1 (MQ=255)
                              tgcCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                                                                                                                         >  1:1442772/1‑150 (MQ=255)
                              tgcCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTac                                                                                                                         >  1:4404905/1‑150 (MQ=255)
                               gcCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTacc                                                                                                                        <  2:4767595/150‑1 (MQ=255)
                                ccGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTacca                                                                                                                       <  1:4072118/150‑1 (MQ=255)
                                   gTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                                                                                                                    <  1:2835643/150‑1 (MQ=255)
                                      aTGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                                                                                                                 <  2:597328/150‑1 (MQ=255)
                                      aTGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                                                                                                                 >  1:203696/1‑150 (MQ=255)
                                        gcgcGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                                                                                                               >  2:21753/1‑150 (MQ=255)
                                         cgcgGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTacac                                                                                                              >  2:1123901/1‑150 (MQ=255)
                                              aTTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGc                                                                                                         >  1:3592711/1‑150 (MQ=255)
                                                tCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCaa                                                                                                       <  2:900808/150‑1 (MQ=255)
                                                tCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCACACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCaa                                                                                                       <  1:2429177/150‑1 (MQ=255)
                                                    gTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTgg                                                                                                   >  2:2828415/1‑150 (MQ=255)
                                                     tCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGc                                                                                                  <  1:666436/150‑1 (MQ=255)
                                                      cTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCg                                                                                                 >  2:4709367/1‑150 (MQ=255)
                                                        ggCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGAc                                                                                               <  1:209916/150‑1 (MQ=255)
                                                         gCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACg                                                                                              <  2:2979382/150‑1 (MQ=255)
                                                         gCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACg                                                                                              >  1:2386250/1‑150 (MQ=255)
                                                         gCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACg                                                                                              >  1:2628270/1‑150 (MQ=255)
                                                            atatTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTg                                                                                           >  2:2676765/1‑150 (MQ=255)
                                                              atTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGaa                                                                                         <  1:3286123/150‑1 (MQ=255)
                                                                tGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGaata                                                                                       >  1:3934303/1‑150 (MQ=255)
                                                                  cccAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGaataat                                                                                     >  1:1768471/1‑150 (MQ=255)
                                                                    cAAACAACGCACTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTg                                                                                   <  2:2084529/150‑1 (MQ=255)
                                                                     aaaCAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTgg                                                                                  >  1:4609573/1‑150 (MQ=255)
                                                                      aaCAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGt                                                                                 >  2:4619944/1‑150 (MQ=255)
                                                                         aaGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACg                                                                              >  1:3622263/1‑150 (MQ=255)
                                                                          aGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTAGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGc                                                                             >  1:1181164/1‑150 (MQ=255)
                                                                             cAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTa                                                                          <  2:3901510/150‑1 (MQ=255)
                                                                              aGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTaa                                                                         <  1:946077/150‑1 (MQ=255)
                                                                                  aaCTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAAc                                                                     <  2:65195/150‑1 (MQ=255)
                                                                                  aaCTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAAc                                                                     <  2:2461264/150‑1 (MQ=255)
                                                                                   aCTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAAcg                                                                    >  2:1757527/1‑150 (MQ=255)
                                                                                    cTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAAcgc                                                                   >  1:4004557/1‑150 (MQ=255)
                                                                                       cGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAAcgcgcg                                                                >  1:2132797/1‑150 (MQ=255)
                                                                                        gACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCgg                                                               >  1:2007401/1‑150 (MQ=255)
                                                                                         aCTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGa                                                              <  1:4741010/150‑1 (MQ=255)
                                                                                             agaaagaaTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACg                                                          <  2:1238926/150‑1 (MQ=255)
                                                                                               aaagaaTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGtt                                                        >  2:2587936/1‑150 (MQ=255)
                                                                                               aaagaaTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGtt                                                        <  2:2323625/150‑1 (MQ=255)
                                                                                                 agaaTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGttct                                                      >  2:919468/1‑150 (MQ=255)
                                                                                                    aTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGttcttcg                                                   >  1:1632309/1‑149 (MQ=255)
                                                                                                    aTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTtcttct                                                   >  1:4192479/1‑150 (MQ=255)
                                                                                                    aTAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTtcttct                                                   >  2:1692394/1‑150 (MQ=255)
                                                                                                     tAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTcttctt                                                  >  1:3865747/1‑150 (MQ=255)
                                                                                                       gCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCtttt                                                <  2:720746/150‑1 (MQ=255)
                                                                                                        cTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTa                                               <  2:491521/150‑1 (MQ=255)
                                                                                                         tGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTCAc                                              <  1:539975/150‑1 (MQ=255)
                                                                                                          gCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACg                                             >  1:1015417/1‑150 (MQ=255)
                                                                                                           ccATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACgg                                            >  1:977127/1‑150 (MQ=255)
                                                                                                              tCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAgg                                         <  1:877197/150‑1 (MQ=255)
                                                                                                                 gtgttgCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCATAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGa                                      <  2:553497/145‑1 (MQ=255)
                                                                                                                    tGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGt                                   >  1:796926/1‑150 (MQ=255)
                                                                                                                     ggCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTa                                  <  2:3657724/150‑1 (MQ=255)
                                                                                                                      gCAAGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTAt                                 <  1:2527914/150‑1 (MQ=255)
                                                                                                                        aaGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTAttt                               >  1:4492827/1‑150 (MQ=255)
                                                                                                                        aaGCGTCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTAttt                               <  1:3450212/150‑1 (MQ=255)
                                                                                                                         aGCGTCACCAGGACGGATTCTTCCCGTCCGTCATACTGGCGCAACTCACGTTGCTGTACCACCAGTACACCACGCACCTGGCGACGTTGACTAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTg                              <  2:3731141/150‑1 (MQ=255)
                                                                                                                           cGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGaa                            >  2:3799135/1‑150 (MQ=255)
                                                                                                                             tCACCAGGAAGGATTCTTCACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGc                          <  2:4404905/150‑1 (MQ=255)
                                                                                                                                ccAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTgg                       <  2:4099617/150‑1 (MQ=255)
                                                                                                                                   ggaaggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGt                    >  1:3558708/1‑150 (MQ=255)
                                                                                                                                   ggaaggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGt                    <  1:3129319/150‑1 (MQ=255)
                                                                                                                                    gaaggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTg                   <  2:4540675/150‑1 (MQ=255)
                                                                                                                                       ggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtt                <  1:1225538/150‑1 (MQ=255)
                                                                                                                                       ggaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtt                <  1:4372975/150‑1 (MQ=255)
                                                                                                                                        gaTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCttt               >  1:639271/1‑150 (MQ=255)
                                                                                                                                         aTTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCtttt              <  2:4505055/150‑1 (MQ=255)
                                                                                                                                            cttcACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTgcg           <  2:4136258/150‑1 (MQ=255)
                                                                                                                                              tcACTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCa         <  2:3611745/150‑1 (MQ=255)
                                                                                                                                                 cTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCt      >  2:404040/1‑150 (MQ=255)
                                                                                                                                                 cTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCt      <  2:18940/150‑1 (MQ=255)
                                                                                                                                                 cTTTAGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCt      <  2:117509/150‑1 (MQ=255)
                                                                                                                                                   ttCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGc    >  2:1794481/1‑150 (MQ=255)
                                                                                                                                                     cGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCaa  >  1:4018912/1‑150 (MQ=255)
                                                                                                                                                     cGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCaa  >  2:2136753/1‑150 (MQ=255)
                                                                                                                                                     cGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCGTTTGCGCATCTGCaa  >  1:2059914/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTTGAATAATTGGTACGCCTAAAAACGCGCGGAAACGTTCTTCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAA  >  CP009273/2961236‑2961534

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: