Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,310,690 C→A 7.2% R889L (CGC→CTC)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,310,6900CA 7.2% 1322.8 / 51.7 345R889L (CGC→CTC) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (136/184);  new base A (6/19);  total (142/203)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 368 total aligned reads displayed.

TAACGTCAGCGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAGC  >  CP009273/2310542‑2310839
                                                                                                                                                    |                                                                                                                                                     
taACGTCAGCGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAAtgcg                                                                                                                                                      <  1:790486/150‑1 (MQ=255)
  aCGTCAGCGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAAtgcgtt                                                                                                                                                    <  2:2579074/150‑1 (MQ=255)
    gTCAGCGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAAtgagttgc                                                                                                                                                  >  1:2770463/1‑150 (MQ=255)
         cGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCaa                                                                                                                                             <  2:4454854/150‑1 (MQ=255)
         cGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCaa                                                                                                                                             >  1:4291266/1‑150 (MQ=255)
          gTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAg                                                                                                                                            >  1:3522688/1‑150 (MQ=255)
            cTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCg                                                                                                                                          >  1:1059823/1‑150 (MQ=255)
             tGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCgg                                                                                                                                         <  2:1559337/150‑1 (MQ=255)
                  tATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAcc                                                                                                                                    <  1:3597495/150‑1 (MQ=255)
                     cacaTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccctc                                                                                                                                 >  1:2570711/1‑150 (MQ=255)
                     cacaTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccatc                                                                                                                                 >  2:787387/1‑150 (MQ=255)
                     cacaTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccatc                                                                                                                                 <  2:4505499/150‑1 (MQ=255)
                      acaTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccatcc                                                                                                                                <  1:3681958/150‑1 (MQ=255)
                       caTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccatcca                                                                                                                               <  1:4559507/150‑1 (MQ=255)
                        aTCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAAccatccat                                                                                                                              <  1:2430926/150‑1 (MQ=255)
                          ccAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCatat                                                                                                                            <  1:3342498/150‑1 (MQ=255)
                            aGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATAttt                                                                                                                          <  2:3854597/150‑1 (MQ=255)
                              cGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTgg                                                                                                                        <  1:2926570/150‑1 (MQ=255)
                              cGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTgg                                                                                                                        >  2:2450716/1‑150 (MQ=255)
                               gTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGc                                                                                                                       >  2:840577/1‑150 (MQ=255)
                                 tACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCAt                                                                                                                     >  1:1483659/1‑150 (MQ=255)
                                    cGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGtt                                                                                                                  <  2:4697495/150‑1 (MQ=255)
                                     ggCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTg                                                                                                                 <  2:2762321/150‑1 (MQ=255)
                                     ggCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTg                                                                                                                 >  2:2832310/1‑150 (MQ=255)
                                       cTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGAc                                                                                                               <  1:1159021/150‑1 (MQ=255)
                                        ttCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGCTGACg                                                                                                              >  1:4509954/1‑150 (MQ=255)
                                          cGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTc                                                                                                            <  2:4247464/150‑1 (MQ=255)
                                            aCAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGc                                                                                                          >  2:1018422/1‑150 (MQ=255)
                                                gCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAg                                                                                                      <  2:1114655/150‑1 (MQ=255)
                                                gCAGCTGTCCATACCGAACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAg                                                                                                      <  2:1156122/150‑1 (MQ=255)
                                                   gCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCAc                                                                                                   <  2:1627251/150‑1 (MQ=255)
                                                        ccATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGatat                                                                                              <  1:1003564/150‑1 (MQ=255)
                                                         cATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATc                                                                                             <  1:3115460/150‑1 (MQ=255)
                                                          aTACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCa                                                                                            >  2:4049206/1‑150 (MQ=255)
                                                          aTACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCa                                                                                            >  1:4779291/1‑150 (MQ=255)
                                                             ccGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAata                                                                                         <  2:3265324/150‑1 (MQ=255)
                                                              cGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAatat                                                                                        <  1:249657/150‑1 (MQ=255)
                                                               ggACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATg                                                                                       <  2:489158/150‑1 (MQ=255)
                                                               ggACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATg                                                                                       >  2:3074262/1‑150 (MQ=255)
                                                                   tCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTc                                                                                   >  2:1603081/1‑150 (MQ=255)
                                                                   tCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTc                                                                                   <  2:3023806/150‑1 (MQ=255)
                                                                    ccAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCt                                                                                  <  1:1849030/150‑1 (MQ=255)
                                                                    ccAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCt                                                                                  <  2:2277132/150‑1 (MQ=255)
                                                                     cAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCtt                                                                                 <  2:1294823/150‑1 (MQ=255)
                                                                      aGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTg                                                                                >  1:4463686/1‑150 (MQ=255)
                                                                      aGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATAGTTGGCATGTTGACGTCGCTAAGCACGATATCACTATGATTCTTg                                                                                >  1:237184/1‑150 (MQ=255)
                                                                         cacCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTa                                                                             <  1:595869/150‑1 (MQ=255)
                                                                           ccGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAg                                                                           <  1:1249483/150‑1 (MQ=255)
                                                                           ccGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAg                                                                           <  2:1740967/150‑1 (MQ=255)
                                                                             gctgctTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTa                                                                         <  1:1419388/150‑1 (MQ=255)
                                                                             gctgctTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTa                                                                         >  2:3486426/1‑150 (MQ=255)
                                                                              ctgctTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTAc                                                                        <  2:1323584/150‑1 (MQ=255)
                                                                              ctgctTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTAc                                                                        >  1:3126838/1‑150 (MQ=255)
                                                                               tgctTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACa                                                                       >  2:679291/1‑150 (MQ=255)
                                                                                      ctTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcg                                                                >  2:455853/1‑150 (MQ=255)
                                                                                      ctTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAgcg                                                                <  1:2463001/150‑1 (MQ=255)
                                                                                        tCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGATGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCa                                                              >  1:1147481/1‑150 (MQ=255)
                                                                                        tCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGATGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCa                                                              >  1:1046280/1‑150 (MQ=255)
                                                                                         cAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCAt                                                             >  1:3775677/1‑150 (MQ=255)
                                                                                          aGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATc                                                            <  2:2369165/150‑1 (MQ=255)
                                                                                           gCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCg                                                           >  1:4221143/1‑150 (MQ=255)
                                                                                               aCGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGc                                                       >  2:1923526/1‑150 (MQ=255)
                                                                                                cGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGcc                                                      >  1:395408/1‑150 (MQ=255)
                                                                                                 gCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCca                                                     <  2:2109480/150‑1 (MQ=255)
                                                                                                 gCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCca                                                     <  2:463849/150‑1 (MQ=255)
                                                                                                  cATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCcat                                                    <  2:355773/150‑1 (MQ=255)
                                                                                                    ttAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCcatca                                                  >  1:1710723/1‑150 (MQ=255)
                                                                                                      aGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCAtt                                                <  2:1208440/150‑1 (MQ=255)
                                                                                                      aGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCAtt                                                <  2:4146138/150‑1 (MQ=255)
                                                                                                      aGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCAtt                                                >  2:3682306/1‑150 (MQ=255)
                                                                                                      aGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCAtt                                                <  2:2711456/150‑1 (MQ=255)
                                                                                                      aGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCAtt                                                >  2:220347/1‑150 (MQ=255)
                                                                                                             aCTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGtt                                         <  1:1174376/150‑1 (MQ=255)
                                                                                                              cTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGttt                                        >  1:4734550/1‑150 (MQ=255)
                                                                                                              cTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGttt                                        <  1:1528865/150‑1 (MQ=255)
                                                                                                                    ttACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACAtt                                  >  1:3145165/1‑150 (MQ=255)
                                                                                                                     tACCGGCAACGTCAGTCCCAACGGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTg                                 <  2:4184707/150‑1 (MQ=255)
                                                                                                                      aCCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTCACCATCCATATTTGGCATGTGGACGTCGCTACGCACGATCTCAATATGATGCTTGCTCAGTGCATTAAGCGCCTCGACGCCCTCCTTCGCGGTTGTACATTGa                                >  2:2938401/1‑150 (MQ=255)
                                                                                                                      aCCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGa                                <  2:4443402/150‑1 (MQ=255)
                                                                                                                      aCCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGa                                >  2:3849268/1‑150 (MQ=255)
                                                                                                                       ccGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGATGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGAt                               <  2:1046280/150‑1 (MQ=255)
                                                                                                                         ggCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAg                             >  1:915158/1‑150 (MQ=255)
                                                                                                                           cAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGcc                           >  1:19239/1‑150 (MQ=255)
                                                                                                                            aaCGTCAGTCCCAACTGACGAATGAGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGccc                          <  2:3126853/150‑1 (MQ=255)
                                                                                                                               gTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAAc                       >  1:1984116/1‑150 (MQ=255)
                                                                                                                                 cAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGa                     <  1:679291/150‑1 (MQ=255)
                                                                                                                                  aGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGAt                    <  1:1275316/150‑1 (MQ=255)
                                                                                                                                       cAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCaa               <  1:2914059/150‑1 (MQ=255)
                                                                                                                                       cAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCaa               <  1:810174/150‑1 (MQ=255)
                                                                                                                                         aCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACt             <  1:3464661/150‑1 (MQ=255)
                                                                                                                                         aCTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCACGGTTTTACATTGATAGCCCAACGATCCCAACt             <  1:412551/150‑1 (MQ=255)
                                                                                                                                           tGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGa           <  1:3330726/150‑1 (MQ=255)
                                                                                                                                             aCGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATc         <  2:1414750/150‑1 (MQ=255)
                                                                                                                                               gAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTg       >  1:3281157/1‑150 (MQ=255)
                                                                                                                                               gAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTg       >  1:4668478/1‑150 (MQ=255)
                                                                                                                                                  tgcgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCa    <  2:1941258/150‑1 (MQ=255)
                                                                                                                                                   gcgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAg   <  2:3818661/150‑1 (MQ=255)
                                                                                                                                                   gcgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAg   <  2:3778144/150‑1 (MQ=255)
                                                                                                                                                    cgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAGc  <  1:2206886/150‑1 (MQ=255)
                                                                                                                                                    cgttgcgtCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAGc  >  2:2481476/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAACGTCAGCGTCTGTTTTATCACATCCAGCGTTACCGGCTTCGACAGGCAGCTGTCCATACCGGACTCCAGACACCGCTGCTTCTCTTCAGCCAACGCATTAGCAGTTACTCCGATTACCGGCAACGTCAGTCCCAACTGACGAATGCGTTGCGTCAAGCGGTAACCATCCATATTTGGCATGTTGACGTCGCTAAGCACGATATCAATATGATTCTTGCTAAGTACATTAAGCGCATCGACGCCATCATTCGCGGTTTTACATTGATAGCCCAACGATCCCAACTGATCTGCCAGC  >  CP009273/2310542‑2310839

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: