Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 453,656 T→A 31.2% S259T (TCC→ACC)  clpX → ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273453,6560TA31.2% 405.1 / 269.9 269S259T (TCC→ACC) clpXATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base T (81/103);  new base A (42/42);  total (123/146)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 275 total aligned reads displayed.

AGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTG  >  CP009273/453507‑453805
                                                                                                                                                     |                                                                                                                                                     
agGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAttt                                                                                                                                                       >  1:553273/1‑150 (MQ=255)
 ggCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTc                                                                                                                                                      >  1:3661881/1‑150 (MQ=255)
   cACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTccc                                                                                                                                                    <  2:1331039/150‑1 (MQ=255)
   cACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTAccc                                                                                                                                                    <  1:4095001/150‑1 (MQ=255)
     cTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCAc                                                                                                                                                  >  1:798469/1‑150 (MQ=255)
         tGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCgtg                                                                                                                                              <  1:2542550/150‑1 (MQ=255)
              cTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaa                                                                                                                                         >  2:30034/1‑150 (MQ=255)
              cTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaa                                                                                                                                         >  1:2292101/1‑150 (MQ=255)
               tGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaaa                                                                                                                                        <  2:1197135/150‑1 (MQ=255)
                 aTCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAAcc                                                                                                                                      >  1:3355143/1‑150 (MQ=255)
                 aTCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAAcc                                                                                                                                      <  1:1222102/150‑1 (MQ=255)
                  tCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTATGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCg                                                                                                                                     <  1:273928/150‑1 (MQ=255)
                   cGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCgg                                                                                                                                    >  1:2008318/1‑150 (MQ=255)
                   cGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCgg                                                                                                                                    >  1:926080/1‑150 (MQ=255)
                    gAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGc                                                                                                                                   >  1:1756058/1‑150 (MQ=255)
                      aGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTc                                                                                                                                 >  2:3725407/1‑150 (MQ=255)
                      aGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTc                                                                                                                                 >  2:4561379/1‑150 (MQ=255)
                       ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGATTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTcc                                                                                                                                <  1:3413329/150‑1 (MQ=255)
                       ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTTTAGAAACCGGCTcc                                                                                                                                <  1:2545731/150‑1 (MQ=255)
                       ggTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTcc                                                                                                                                <  2:3387873/150‑1 (MQ=255)
                        gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCg                                                                                                                               <  2:1516211/150‑1 (MQ=255)
                        gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCg                                                                                                                               <  1:3530562/150‑1 (MQ=255)
                           cGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCa                                                                                                                            <  1:4243145/150‑1 (MQ=255)
                              tAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTg                                                                                                                         >  1:1560131/1‑150 (MQ=255)
                              tAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTg                                                                                                                         <  1:4685010/150‑1 (MQ=255)
                               aGCTGCTGTTCCACCGCATGGTGGGCGTAAACATCCGCAGCAGGTATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTgg                                                                                                                        <  2:4168317/150‑1 (MQ=255)
                               aGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTgg                                                                                                                        <  2:2714387/150‑1 (MQ=255)
                               aGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTgg                                                                                                                        >  1:1086937/1‑150 (MQ=255)
                               aGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTgg                                                                                                                        >  1:3471992/1‑150 (MQ=255)
                                 ctgctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGtt                                                                                                                      <  2:1471235/150‑1 (MQ=255)
                                  tgctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGttt                                                                                                                     <  1:916558/150‑1 (MQ=255)
                                     tgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTgg                                                                                                                  >  2:1828437/1‑150 (MQ=255)
                                         ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGTAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGcgcg                                                                                                              <  1:2072251/150‑1 (MQ=255)
                                          cACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGa                                                                                                             <  2:4283151/150‑1 (MQ=255)
                                                aaGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCTCGACGGTaa                                                                                                       >  2:3877417/1‑150 (MQ=255)
                                                aaGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTaa                                                                                                       >  1:2605957/1‑150 (MQ=255)
                                                aaGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCAGTGGTTTTGGCTCGACGGTaa                                                                                                       >  2:4163892/1‑150 (MQ=255)
                                                  ggTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCAGTGGTTTTGGCGCGACGGTaaaa                                                                                                     >  1:151303/1‑150 (MQ=255)
                                                    tGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGc                                                                                                   >  2:2632781/1‑150 (MQ=255)
                                                     gggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCg                                                                                                  >  2:805219/1‑150 (MQ=255)
                                                     gggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCg                                                                                                  >  1:3072458/1‑150 (MQ=255)
                                                        cttAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAg                                                                                               <  1:4453292/148‑1 (MQ=255)
                                                           aaaCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTcc                                                                                            >  2:3118056/1‑150 (MQ=255)
                                                            aaCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCg                                                                                           >  1:4614869/1‑150 (MQ=255)
                                                             aCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGa                                                                                          <  2:2216147/150‑1 (MQ=255)
                                                               aTCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACa                                                                                        >  1:2342345/1‑150 (MQ=255)
                                                               aTCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACa                                                                                        >  1:4356488/1‑150 (MQ=255)
                                                                 ccGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACaaa                                                                                      >  1:312920/1‑150 (MQ=255)
                                                                  cGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaag                                                                                     <  1:1377722/150‑1 (MQ=255)
                                                                  cGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaag                                                                                     <  2:2530307/150‑1 (MQ=255)
                                                                   gcagcaGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaagc                                                                                    <  1:2467832/150‑1 (MQ=255)
                                                                     agcagGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaagcaa                                                                                  <  2:4063506/150‑1 (MQ=255)
                                                                      gcagGAATTCTTGCAGGTTGAGACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaagcaag                                                                                 <  2:3277085/150‑1 (MQ=255)
                                                                            aTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAgg                                                                           >  1:3681482/1‑150 (MQ=255)
                                                                             ttCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGc                                                                          <  2:3471992/150‑1 (MQ=255)
                                                                               cTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGa                                                                        >  2:488458/1‑150 (MQ=255)
                                                                                  gCAGGTTGATACCTCTACGAGCCTGTTTATTTGTGGGGGTGCGTTTGCCGGTCGGTTTAATGTGATTGCCCACCGTGTAGAAACCGGCTCCGGCTTTGGTTTTGGCGCGACGGTAAATGGGAAGTCCGACAAAGCATGCGAAGGCGAgct                                                                     <  2:645398/150‑1 (MQ=255)
                                                                                       ttGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTgg                                                                <  1:4122797/150‑1 (MQ=255)
                                                                                        tGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgc                                                               >  2:746385/1‑150 (MQ=255)
                                                                                        tGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgc                                                               <  2:2008318/150‑1 (MQ=255)
                                                                                         gATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgcg                                                              <  2:3983705/150‑1 (MQ=255)
                                                                                           tACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCa                                                            <  2:888395/150‑1 (MQ=255)
                                                                                           tACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCa                                                            <  1:758686/150‑1 (MQ=255)
                                                                                             ccTCGAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAgg                                                          <  2:1343681/150‑1 (MQ=255)
                                                                                              ctctGAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGt                                                         >  1:2849839/1‑150 (MQ=255)
                                                                                                   aGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAAc                                                    >  1:393004/1‑150 (MQ=255)
                                                                                                       ccTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGa                                                >  1:1827495/1‑150 (MQ=255)
                                                                                                        cTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGaa                                               <  2:3427900/150‑1 (MQ=255)
                                                                                                         tGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                              >  2:2719896/1‑150 (MQ=255)
                                                                                                           tttATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                            <  2:2911838/150‑1 (MQ=255)
                                                                                                           tttATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGAt                                            <  1:1390281/150‑1 (MQ=255)
                                                                                                             tATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCt                                          <  2:553273/150‑1 (MQ=255)
                                                                                                             tATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCt                                          >  1:1651475/1‑150 (MQ=255)
                                                                                                                  gtgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCa                                     <  1:488458/150‑1 (MQ=255)
                                                                                                                   tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCaa                                    >  2:2430882/1‑150 (MQ=255)
                                                                                                                   tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCaa                                    >  1:984116/1‑150 (MQ=255)
                                                                                                                    gtcGGTGCTTTTGCGGTTCTGTATAATGTGATTACCCTCTGTGTAGAAACCGGCTCCGGCATTGTTTTTGGCGAGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAg                                   <  1:855707/148‑1 (MQ=255)
                                                                                                                    ggCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAg                                   <  2:3559335/150‑1 (MQ=255)
                                                                                                                      cGGTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGtt                                 <  2:4271712/150‑1 (MQ=255)
                                                                                                                       ggTGCGTTTGCCGGTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGttt                                <  2:1098573/150‑1 (MQ=255)
                                                                                                                         tGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTgg                              <  1:124134/150‑1 (MQ=255)
                                                                                                                          gCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGt                             >  1:376053/1‑150 (MQ=255)
                                                                                                                             tttGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCtt                          <  2:831899/150‑1 (MQ=255)
                                                                                                                               tGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTAt                        >  1:1083334/1‑150 (MQ=255)
                                                                                                                                   ggTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                    >  2:1268536/1‑150 (MQ=255)
                                                                                                                                   ggTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                    <  2:1772468/150‑1 (MQ=255)
                                                                                                                                   ggTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCACGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGACGATCTGATCAAGTTTGGTCTTATCCCt                    >  2:1930953/1‑150 (MQ=255)
                                                                                                                                   ggTCTGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCt                    >  2:4330662/1‑150 (MQ=255)
                                                                                                                                     tCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGa                  <  2:926080/150‑1 (MQ=255)
                                                                                                                                       tGGATAAAGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                >  2:3041487/1‑150 (MQ=255)
                                                                                                                                          aTAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTa             <  2:2849631/150‑1 (MQ=255)
                                                                                                                                          aTAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTa             <  2:3619800/150‑1 (MQ=255)
                                                                                                                                             aaGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTg          <  2:2304427/150‑1 (MQ=255)
                                                                                                                                             aaGTGATTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTg          >  2:738931/1‑150 (MQ=255)
                                                                                                                                               gTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgt        <  2:1051156/150‑1 (MQ=255)
                                                                                                                                                tGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtc       <  2:3355143/150‑1 (MQ=255)
                                                                                                                                                 gATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtcg      <  1:2715919/150‑1 (MQ=255)
                                                                                                                                                 gATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtcg      >  1:626624/1‑150 (MQ=255)
                                                                                                                                                  aTTACCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtcgt     <  2:4403850/150‑1 (MQ=255)
                                                                                                                                                   tttCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtcgtc    >  2:2800776/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTG  >  CP009273/453507‑453805

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: